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CONSERVANCY DISTRICTS AS FLOOD CONTROL ORGANIZATIONS# 


By Cloyde C. Chambers,1 F. ASCE 


SYNOPSIS 


A well designed conservancy district is a valuable organization for de- 
loping, constructing and administering a comprehensive water management 
ogram for an intrastate watershed. 

The evidence, as provided by many operating districts, indicates that this 
pe organization is well adapted to the solution of flood problems, especially 
1en approached on a drainage area basis. It is most applicable in the field 
at lies between those problems that can be solved best be existing local 
vernmental agencies and the problems of the major rivers, which are inter- 
ate in character and in which a federal interest is predominant. 

Some of the most successful flood control projects have been achieved by 
mservancy districts which have kept a large degree of control of the plan- 
ng and operation at the local level and which have assumed a fair share of 
e project costs. Local project control and maintenance can be provided best 

there is a conservancy district or similar organization with the legal 
thority and financial ability to participate in the project, even though there 
e large state and federal interests involved. 

A conservancy district, to operate successfully, must be a political sub- 
vision of the state, organized on a problem area basis. It must be organized 

the request of and have the support of the local people. It must have 
thority to plan and to execute plans. It must have financial responsibility 
d the power to levy taxes and assessments. It must have authority to ac- 
ire property, including the right of eminent domain, and power to establish 


Note.—Discussion open until September 1, 1960. To extend the closing date one 

onth, a written request must be filed with the Executive Secretary, ASCE. This paper 
art of the copyrighted Journal of the Hydraulics Division, Proceedings of the Amer- 
1 Society of Civil Engineers, Vol. 86, No. HY 4, April, 1960. 

a Presented at the May, 1959 ASCE Convention in Cleveland, Ohio. 

1 Cons. Engr., Dayton, Ohio. 


i 1 


2 April, 1960 HY 


and enforce rules and regulations. It must, also, have authority to cooperat 
with other levels of government and their agencies. 


INTRODUCTION 


The engineer engaged in public works must accomplish his goals throug 
the medium of diverse types of public agencies and organizations which ma 
have a variety of legal powers, administrative authorizations and objective: 
The organizational means and available legal powers for solving his probler 
and constructing his works are, therefore, as important to him as the mé 
terials, tools and equipment he selects or designs. He should have a know: 
edge of the advantages and disadvantages of the various types of public orgam 
zations available for his use. This is particularly true in dealing with floc 
control problems, each one of which is generally unique in character. H. 
recommendations, as to the provisions to be contained in legislation creatir 
the powers, authority and organization with which he will have to work, shou 
be given much weight, because of the effect the legislation will have on tl 
quality of the work to be produced. 

In drafting legislation for flood control purposes, it is important that prc 
visions for the administrative organization be designed to be adaptable to 
variety of flood control problems. The enactment of the first enabling legis 
lation, making possible the creation of conservancy districts, required ; 
much pioneering as the technical problems involved in the design of projec 
to be carried out under the legislation. Negligence in either case would hav 
increased the chance of failure. 


USE OF TERMS 


The term ‘conservancy district,” for the purposes of this paper, is a] 
plied to a variety of special local organizations created by legislation : 
special taxing and assessment districts with designated powers and some | 
cal financial responsibility. 

The term “flood control” is interpreted to cover the broad field of mea: 
ures that can be used to prevent or mitigate losses and other ill effects, whi: 
normally occur asa result of heavy precipitation. These measures can | 
grouped, roughly, into three main classifications (the first has a tendency 
merge with the second): 


1. Flood prevention measures, which include, among others, those mea 
ures for preventing or retarding the runoff of excess water from the land a 
for retarding the flow of small streams or brooks. 

2. Flood control measures, consisting principally of works that are d 
signed to control the flow of flood waters by the temporary storage of exce 
water in retarding basins or reservoirs, the structures of which can have 
wide range in size. 
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3. Flood protection measures include local works at critical or important 
Qints in the flood plains, such as channel improvements, protective levees, 
alls, etc., either as a sole remedy or for Supplementing one or both of the 
ther types of measures. 


CONSERVANCY DISTRICTS 


Conservancy districts vary in organization, powers and objectives, in the 
everal states where such districts are found, because the wording of the 
espective enabling legislations, and the authorities given to the districts are 
uite different. In some cases, this variation is the result of differences in 
le state constitutions and the existing laws to which the districts must con- 
rm. In other instances, the conditions and problems to be met are not the 
ame. States with semi-arid areas may give irrigation a prominent place 
long with flood control, as, for example, the New Mexico law under which the 
fiddle Rio Grande Conservancy District has operated. Power may be given 
nportance in others, as in the New York law, under which the Hudson River 
egulating District operates. 

Some of the provisions, such as method of organizing, assessing for bene- 
ts, etc., that are now incorporated in conservancy acts have been used for 
any years in other enabling legislation for the solution of flood problems. 
aws for the creation of levee and drainage districts date back to the days of 
e early settlements of the country. These earlier laws provided for dis- 
‘icts with a limited scope of activities and frequently authorized project con- 
Truction without adequate provisions for maintenance. The boundaries of 
ich districts were often limited by lines dividing other existing political 
ibdivisions. The enabling legislation was frequently drawn with a particular 
‘oblem in mind, consequently it did not fit other problems, or it lacked some 
ements that would have provided for good flexible operation. Eventually, the 
ed for legislation with a broader viewpoint, and with provisions for coordi- 
iting a variety of related problems, became apparent. 

Flood problems and their solutions are generally unique. It is not often 
at blueprints, successfully used at one location, can be moved bodily to 
other problem area and achieve the same degree of success. 

A conservancy district law should be broad enough to permit a complete 
d honest inventory of the flood problems in the various parts of the district 
d to give freedom to ferret out and explore the various possibilities for the 
lution of the problems. The best solution for major flood problems is 
ually found in a variety of measures rather than in a single one. 

The need for an organization that can fill the gap between those problems 
at can be solved by existing local agencies and the interstate problems on 
r larger river systems, in which the federal government has an interest, 
Ss been emphasized by the interest shown in conservancy districts in Ohio 
d other states. In addition, it has been shown that conservancy districts are 
metimes useful in the solution of the local flood problems and the problems 
river systems, such as the Muskingum watershed, which are tributary to 
2 larger systems, such as the Ohio and Mississippi Rivers. 
Jt is rare that existing political subdivisions have the facilities, authority 
d flexibility of action necessary to cope even satisfactorily with local prob- 
ms. It is not often that the boundaries of flood problem areas are common 
the boundaries of such subdivisions. Flood control and water use problems 
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tend to divide themselves into stream drainage areas. The Miami Cc 
servancy District and the Muskingum Watershed Conservancy District ¢ 
examples of situations in which a group of several communities in one drai 
age area have been given an opportunity to pool their problems for unifi 
solution. Satisfactory solutions for the problems in those districts would ha 
been impossible, had such an opportunity not been created. 

There are examples of communities electing to stay out of conservar 
districts covering large drainage areas in order to organize their own sm 
districts to provide for their own local protection. Such protection has ge 
erally been less adequate than that which would have been provided by 1 
larger district. As a result, these communities, as, for example, Mt. Vern 
Ohio, in the Muskingum area, and Springfield, Ohio, in the Miami area, he 
suffered repeated, serious flood losses because of inadequate protection. 


OHIO EXPERIENCE 


Interest in enabling legislation for the creation of conservancy distri 
first became active in Ohio, immediately following the floods of March, 19 
The enactment of the Ohio conservancy act was prompted by a realization th 
in some cases, it would be desirable to work out flood problems on a ri' 
drainage-wide basis, rather than solely by individual local projects, and 
permit a group of communities in one valley to pool their problems for unif 
solution. Without such enabling legislation it seemed impossible to solve so 
of the flood problems of the state. The persistent and well-guided determil 
tion of the people in the Miami Valley to do something about the flood probl 
in that area gave strong impetus to the enactment of the first Ohio cons 
vancy act. It was a case of a new invention being developed to meet a new 
urgent need. 

The acceptance of the conservancy type organization for flood control ; 
allied problems has been more pronounced in Ohio than in some other sta 
because of the early recognition of the need for such an organization. 
greater use in that state has been influenced to some extent by the success 
the Miami Conservancy District. The Miami project was financed with fu 
derived from assessments on appraised benefits to property in the distr: 
The original cost of the project has now been paid in full. The Upper Sci 
Drainage and Conservancy District, the first to be organized in Ohio, has a 
been financed solely with assessments levied on benefited property. Both d 
tricts have remained active and have provided a high standard of maintena 
for their projects with funds raised through maintenance assessments 
property receiving project benefits. 

A total of twenty-two districts and four sub-districts have been organi 
under the Ohio act. They range in size from a district involving a small gr 
of farmers to one that includes about one-fifth of the state. Fourteen of 
twenty-two districts have included flood control along stream valleys as” 
of their principal purposes. Two others include flood control and beach e 
sion prevention along the shore of Lake Erie, as their prime objectives. 

Of the fourteen districts that included flood control along stream valleys 
a principal purpose, seven have been completed or are in the process of ¢ 
structing their initial programs; two are in the development stage; one is 
active; and four are either in the process of being dissolved or have beenc 
solved. 7 


Ty 4 CONSERVANCY DISTRICTS 4) 


Some districts have become inactive or disorganized because it was found 
hat the estimated costs of the proposed work either exceeded the estimated 
enefits to be derived therefrom, or were higher than the people would accept. 
nm some instances, it was found that the proposed plans did not offer a suita- 
le or acceptable solution to the problems, and in others, the districts were 
wromoted without conducting a thorough and honest educational program among 
he local people most vitally affected by the proposed program. 


COMPONENTS OF GOOD CONSERVANCY LAWS 


Current thinking on the development, management and control of the na- 
ion’s water resources requires that the enabling legislation for the creation 
f conservancy districts provide for a broad view of the complicated and in- 
errelated water, soil and sanitary problems, and for considerable latitude in 
uthority for solving such problems. 

As conservancy districts generally do not coincide with the usual political 
Oundaries, it is necessary to provide for their establishment through special 
egislation for creating a separate political subdivision of a state. Such a type 
f organization is also necessary because of taxing and other powers that dis- 
ricts must have in order to operate. 

The enabling acts must provide for the free expression of the desires of the 
eople and for an opportunity for them to accept or reject any proposed dis- 
rict. Objections to conservancy districts are very common and arise for 
many reasons. 

Some conservancy district enabling laws have been criticized on the 
rounds that they are too long or complicated (the Ohio law, for example) and 
hat the steps required for development and operation of a program are too 
umerous and time consuming, especially for small or local projects. On the 
ther hand, attempts at brevity generally omit some one or more essential 
lements of a workable law. In addition, after a proposed development passes 
hrough all the checks and balances of a well-drafted conservancy district law, 
ts soundness is fairly well established from a practical as well as a techni- 
al viewpoint. 

Sometimes, objections to conservancy acts are made because the acts are 
idged by the provisions of individual sections rather than as an act composed 
f many interrelated working parts. Objections are also raised over the use 
fa conservancy district for the solution of local flood problems because the 
ork can be done through existing governmental agencies, thus avoiding the 
reation of a new organization. 

There are some advantages in the use of conservancy districts for some 
cal problems. As previously mentioned, boundaries of existing agencies do 
ot often fit desirable boundaries for flood problem solutions. In addition, the | 
lanaging directors of a district usually serve out of a civic interest or duty 
ather than for any financial or political gain. Their attention is directed to- 
ard a single objective and its allied features, and they are not burdened with 
ther day-to-day routine activities. 

The degree of influence a conservancy district may have on the solution of 

ood control and allied problems may depend, in some cases, upon the poli- 

ies of the state and federal governments at the time the project is undertaken. 

At the time when the conservancy district idea began to be developed, the 

ctive interest of the federal government in flood control was confined largely — 
= 
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to the main rivers that had been or were being developed for general navigz 
tion. Then the field was open for conservancy districts to operate under tt 
original concept of a self-financing organization. 

Changes in federal and state policies, with respect to flood control sin« 
the conservancy district idea was first brought into prominent use, especial. 
during the past quarter of a century, have tended to change the function of suc 
districts, but they have not entirely eliminated the need for them. Many of tt 
newer districts are being established to cooperate on federal projects rathe 
than to operate independently. 

Persons who believe that the entire development, planning, design, cor 
struction, maintenance and operation of all flood control projects, regardles 
of type and size, are the sole responsibility of federal and state government: 
may find little use for conservancy districts as flood control organization 
The writer, who is an advocate of retaining a high degree of local respons: 
bility and control, does not subscribe to that belief. 

In some river drainage areas, it is unwise to divorce the flood problen 
in the area from the needs for the development and wise use of the water re 
sources of that area. In other areas, the water problems are tied closely | 
the problems of good management of other natural resources. A conservant 
district created for purposes that will include all related and associated prol 
lems is in a better position to develop or guide the development of a compre 
hensive plan that will have a broad scope and balance between functions thz 
is an organization that has a predominant interest in only one phase of tl 
over-all problem. 

Many valleys, where conservancy districts are organized for flood contr 
purposes, are experiencing an increasingly rapid and complex econom 
growth. These developments create water, sanitation, highway and othe 
problems, many of which have important relationships with the work of tl 
districts. In some states, very moderate amendments to their conservan 
acts would enable conservancy districts to better serve those communitii 
having such closely related problems. At least, the conservancy distri 
might be the organization for dealing with all water problems, rather than wi 
a single purpose such as flood control. 

There have been alterations, amendments, and, in some cases, redraftil 
of conservancy acts, where experience has dictated it to be necessary, in 0: 
der to meet changing or overlooked conditions. The original concept for co: 
servancy districts was for flood control, with only rather general provisio 
for other conservation and water problems. The tendency in amendments h: 
been to expand and detail those provisions. Experience with conservancy di: 
tricts has not suggested the need for radical changes in the original concefy 
However, in some cases, there is a need for new powers or a new approach 
implement that original concept, as for example, how to assess or charge f 
water use in a ground-water replenishment project. 

Conservancy districts should have some sort of zoning and regulating ar 
thority, in order that serious flood problems can be prevented before th 
arise, and in order that encroachments and objectionable obstructions to flo 
flow are not built in the flood plains. This authority may be exercised, eith 
alone or in cooperation with some other proper agency, to restrict develo 
ments in flood plains until it is economically and otherwise feasible to provi 
adequate flood protection. ag 

Most conservancy districts have the authority to move and reestabli 
communities that will be disturbed by the construction of district project 
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This authority should be broadened to permit moving communities out of flood 
plains, where it is economically and otherwise practical to do so, and where it 
is more feasible than protecting them against flooding. 

In some cases, through the authority of conservancy districts, the pattern 
of railroad tracks has been changed to permit the construction of the most 
feasible plan for flood protection and, with a little forethought, at the same 
time, to provide a more efficient plan viewed from railroad and general wel- 
fare viewpoints. Examples are the work done by conservancy districts at 
Pueblo, Colorado and Massillon, Ohio. 

In order to carry out its work, a conservancy district must be able to raise 
funds. This is usually accomplished either through an ad valorum tax on the 
taxable property in the district, by assessments made on the property in pro- 
portion to the benefits received, or both. These methods are used to pay for 
the preliminary expense of organizing the district, the preparation of the plan 
of works, the construction of the works and their operation and maintenance. 

Some states have had difficulty in incorporating within their constitutional 
framework legislation that will provide an equitable and legal means of dis- 
tributing the preliminary expense of organization of the district, preparing 
plans, making appraisals and other expenses prior to the receipt of funds from 
assessments on appraised benefits. Some suggest that these expenses be met 
with appropriations by the state. This is equitable only in case there is defi- 
nite state interest in developing the possibilities of the proposed project and 
should be done in a manner that does not diminish local control over the af- 
fairs of the district. Preliminary funds should be made available promptly 

after the district is organized in order that planning can go forward while ini- 
tial enthusiasm is high. 
: The problem of preliminary funds has been solved successfully in several 
cases by advances of funds or loans made by other agencies, such as counties, 
“Organizations, such as Chambers of Commerce, or individuals, with the under - 
standing that such funds will be repaid when and if funds are made available 
from assessments on benefits or from other sources. The State of Indiana, 
for example, has made loans to conservancy districts for preliminary ex- 
-penses from its Flood Control Revolving Fund. Financing in this manner is a 
test of the local appraisal of the merits of the proposed program and the sin- 
-cerity of the local sponsors. 
: A conservancy district has an advantage over some other types of districts 
‘and some local governmental agencies in the provisions for maintenance of 
the works of the district. The right of the district to make annual mainten- 
ance assessments onbenefited properties provides funds to keep the improve- 
ments well maintained. This helps to keep the works ina high state of effi- 
ciency at all times, and at less expense through the correction of maintenance 
‘needs while they are small instead of by large maintenance or rebuilding jobs 
at rare intervals. The contrast between the types of maintenance usually pro- 
ided by a conservancy district and that of some existing local governmental 
agencies is often very striking. 
i The function of conservancy districts as organizations to raise funds for 
the construction, maintenance and operation of projects may have been re- 
duced as federal and state appropriations have become increasingly available 
for such purposes. However, the need for such organizations to provide the 
local cooperation for many projects still exists. 

An equitable subdivision of costs among the federal government, state and 

local interests, in proportion to their respective benefits, may be appropriate 
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in some cases. The appropriations by the federal government have tended to 
exceed their fair share of the cost. The tendency for states to make appropri- 
ations to reduce or eliminate local flood assessments is becoming more 
prevalent. The need for an alert local agency is as great now as it ever has 
been, not only as an agency to provide those funds required locally, but also 
to protect the local interests. There may seem to be little need for a conser- 
vancy district in the case of asystem of federal flood control reservoirs to be 
constructed on a main tributary river system in the area but designed almost 
exclusively for benefits to be provided outside the area. However, a local or- 
ganization, with the broad powers that a conservancy district should have, 
should be able to ferret out any alternate or amended plan which would give 
comparable outside benefits and, at the same time, be more desirable from a 
local viewpoint. 

Federal participation in the cost of channel improvement, levee and other 
local projects is generally limited, under present federal policy, to the first 
construction cost. Some local agency must assume the responsibility for ac- 
quiring the necessary lands, rights-of-way and easements for constructing the 
project and must agree to maintain and operate it after it is built. Similar re- 
quirements must be met for projects in which the United States Department of 
Agriculture participates, under Federal Statute Public Law 566. These are 
proper functions of a conservancy district, and if administered by such an or- 
ganization will insure the greatest value to the respective areas from the 
federal contribution. 

A conservancy district is in a good position, also, to coordinate the activi- 
ties of various agencies and units of government interested in flood contro! 
and other water problems in its area. It can provide the assurance and co- 
operation required on federal projects through one organization, and eliminate 
the need for the federal governments to enter into many agreements. 

The broad powers generally vested in a conservancy district give oppor- 
tunities to provide other benefits when local flood protection works are buil 
through an urban area. With proper planning, unpleasant river fronts can be 
changed into the most attractive parts of a community. In some cases, the 
authority of conservancy districts should be broadened to permit more active 
cooperation in other community activities closely allied to flood control. 


SUMMARY 


A conservancy district, to serve present day needs, must be adaptable to ; 
broad range of purposes. It is not enough that they provide for flood contro 
alone; they must be able to cope with many other water problems and relate 
activities. 

Enabling legislation for the creation of a conservancy district should pro 
vide for at least the following features and powers: 


1. It should be a political subdivision of the state, be created under stat 
law, and should have perpetual existence except as provisions for its dissolu 
tion are provided in the enabling legislation. 

2. The initiative for the creation of a conservancy district should com 
from local interested persons, corporations or political subdivisions by peti 
tion, vote or otherwise, and those opposed to it should have an opportunity t 
register and argue their opposition. 
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3. The purposes for which a district may be organized should include the 
prevention or control of floods and the development and wise use of the water 
resources of the district, in accordance withneeds when such development can 
be accomplished to good advantage along with flood control. 

4. It should have the authority to (1) develop, or cause to be developed, a 
plan for the improvements necessary to meet the purposes for which the dis- 


trict is created, (2) cause such a plan to be officially adopted after passing 
through the democratic processes necessary to give proper checks and 
balances, and (3) amend the plan from time to time, as needed. 


5. It should have the authority, alone or in cooperation with other agencies, 
to construct, maintain and operate its works and improvements. 

6. It should have the authority to establish and enforce the rules and regu- 
lations necessary to accomplish its purposes and protect its works. 

7. It should have the authority to control, hold, and acquire by purchase, 
condemnation, donation or otherwise, all lands and rights necessary to con- 
struct, use and maintain its works. 

8. It should have the right of eminent domain, including a dominant right 


over utilities and subordinate public corporations. 


9. It should have the authority to secure by special assessments the funds 
necessary to carry out its purposes, and the authority to issue bonds in antici- 
pation of their collection. 

10. It should have the authority to enter into contracts or other arrange- 
ments with the United States and state governments, or any agencies thereof, 
and with persons and corporations, for cooperation and assistance in carrying 
out any of its purposes, objectives and authorities. 


Like many other undertakings, the degree of success that can be expected 
from a conservancy district will depend largely upon how well it is adminis - 
tered. A conservancy district created under an ideal law can be only moder- 
ately satisfactory if it gets into incompetent or insincere hands. On the other 
hand, one created under a law that may lack some desirable features can be 
successful when used by a community with vision and administered by compe- 


tent and sincere leadership. Such leadership exists in nearly every com- 


munity. In any case, success will require sacrifices and services from people 
who will not receive their full reward in financial compensation. 

The final worth of a water control and use project will depend upon how 
well society is able to adapt itself to the benefits the project offers. There 
have been projects of sound economics and engineering which have had dis- 
appointing results because the people in the area have been unable to use the 
opportunities offered to good advantage. Less attractive projects have be- 
‘come worthy because the people have adapted every possibility to good use. 
The real satisfaction of the project builders does not come at the comple- 
tion of the project, but rather when they go back after many years and find that 
‘the people have been able to take full advantage of the opportunities that were 
offered, and when they see that the test of time has proved their efforts to be 
‘worthwhile. The more the most directly benefited individuals have contributed 
‘toward the initiation, development, construction, maintenance and operation of 
aconservancy district project, the better chance there is for the realization 
of such satisfaction. A conservancy district can be a worthwhile tool for a 


ommunity with foresight. 
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A COMPARISON OF STREAM VELOCITY METERS 


By F. Wayne Townsend! and F. A. Blust2 


SYNOPSIS 


A description of the U. S. Lake Survey method of making stream flow 
Measurements is given, and comparisons of current velocities simultaneously 
: measured by cup and screw type current meters in the lower Niagara River 
are presented. The conclusion is reached that the two types of meters give 
identical results in Lake Survey flow measurements. 


INTRODUCTION 


Among the responsibilities of the United States Lake Survey is the meas- 
‘urement of the flows in the outflow rivers of the Great Lakes, which include 
‘the St. Marys, St. Clair, Detroit, Niagara, and St. Lawrence rivers. Such 
“measurements were first made by the Lake Survey in the 1860’s and have been 
‘made many times since. Fig. 1 shows the Great Lakes and their outflow 
‘Tivers. In recent years, vastly increasing hydro-electric power development 
and the agreements incident thereto between the United States and Canada for 
the sharing of the Great Lakes outflows have rendered the accurate deter- 
mination of these flows more important than ever before. 
ce es es eS ee 
Note.—Discussion open until September 1, 1960. To extend the closing date one 
month, a written request must be filed with the Executive Secretary, ASCE. This paper 
is part of the copyrighted Journal of the Hydraulics Division, Proceedings of the Ameri- 
can Society of Civil Engineers, Vol. 86, No. HY 4, April, 1960. 


1 U.S, Lake Survey, Corps of Engrs., U.S. Army, Detroit, Mich. 
2 U.S, Lake Survey, Corps of Engrs., U.S. Army, Detroit, Mich. 


11 


Pe on man, 


12 April, 1960 HY 4 


Recent flow measuring practice by the Lake Survey has involved the use of 
cup-type current meters which respond to the water movement by means of a 
number of cups revolving around a common vertical axis. There is considera- 
ble evidence3-9 that a revolving cup meter will over-register, that is, give a 
velocity higher than actual in turbulent water or as a result of movement of 
the meter. There is also evidence3,4,5,6,8,9 that a screw-type current meter 
which responds to the water movement by means of a number of vanes re- 
volving around a common horizontal axis does not have this tendency to over- 
register. Because of the present need for greatest possible accuracy of flow 
measurements, a comparison was made of current velocities obtained simul- 
taneously by cup meters in present use and screw meters of latest design. 
The data on which this comparison was based were obtained under conditions 
typical of those encountered by the Lake Survey in measuring flows of the 
outflow rivers of the Great Lakes, and the comparison is intended to show only 
whether or not the two types of meters will give significantly different re- 
sults under these conditions. 
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Lake Superior 


St. Marys River 


Straits of Mackinac 


St. Clair 
River 


“fr . Niagara River 
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FIG. 1.—THE GREAT LAKES AND THEIR OUTFLOW RIVERS 
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This paper presents first a brief description of the Lake Survey method of 
flow measurement to help the reader understand the comparison that was 
made, and then a description and analysis of the comparison. 


UNITED STATES LAKE SURVEY METHOD OF FLOW MEASUREMENT 


In selecting a vertical plane in the river, called a hydraulic section, past 

which the flow is to be measured, the flow characteristics looked for are lack 
of turbulence, uniform direction parallel to the banks, and velocities ranging 
between 1.5 and 5.0 ft per sec. Channel configurations sought are parallel 
banks, uniform cross sections upstream and downstream from the section, and 
relatively steep slopes at the water’s edge. A section having these flow and 
‘channel characteristics can be selected tentatively for a given reach of the 
iriver by means of field reconnaissance and inspection of available hydrogra- 
phic and topographic surveys. Before a section is finally selected, however, 
exploratory soundings and current velocity observations are made at intervals 
of 100 ft above and below the proposed section to verify the previous findings. 
The nature of the Great Lakes’ connecting rivers is such that an acceptable 
hydraulic section usually can be found in any desired reach. 
_ In the Lake Survey method, the hydraulic section is divided by vertical 
lines into a number (usually 10) of areas of nearly equal width, called panels. 
Current velocity measurements are made at a single specific location in each 
successive panel starting at one side or the other of the river. This location 
is usually at the horizontal mid-point of the panel and at a depth of 0.4 of the 
total depth at that point. A set of such velocity measurements for all panels in 
the hydraulic section, requiring two to three hours to obtain, is considered a 
single flow measurement, and the total flow computed from the data obtained 
is considered to be the average flow existing during the flow measurement. 

The flow in cubic feet per second through each panel is computed by multi- 
plying the measured velocity in feet per second by the area of the panel in 
square feet and by three coefficients which are called the directional coeffi- 
cient, the vertical coefficient, and the horizontal coefficient. It is the purpose 
of the coefficients to convert the measured velocities into the mean velocities 
in each panel. The total river flow is computed by adding the panel flows. 
The direction of flow through each panel is determined by transit tracking 
of floats, and the sine of the angle between the flow direction and hydraulic 
section is the directional coefficient. Normally, the directional coefficients 
are determined before the first flow measurements at a new hydraulic section 
are made, and they are not re-determined as long as the physical character - 
istics of the river in the vicinity of the hydraulic section remain the same. 
A vertical coefficient for each panel is determined from a number of si- 
multaneous observations of velocity at the measuring point of the panel at 0.4 

lepth and at each of the fractional depths 0.1, 0.2, 0.3, 0.5, 0.6, 0.7, 0.8, and 
‘4 The average of the ratios of the velocities at the various depths to the 
simultaneous velocities at 0.4 depth is the vertical coefficient for the panel. 
Data for determination of vertical coefficients are usually obtained for a new 
| ae section before the first flow measurements are made, and the de - 
ived coefficients, like the directional coefficients, are not re-deter mined as 
ong asthe physical characteristics of the river in the vicinity of the hydraulic 
ipetion remain the same. Under some circumstances, measurements of flow 
we made by metering the current velocities simultaneously at 0.2 and 0.8 
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depth in each panel instead of 0.4 depth. When this is done, vertical coeffi- 
cients are not used; the flow through each panel is computed by multiplying the 
mean of the velocities simultaneously measured at 0.2 and 0.8 depth by the 
panel area and by the directional and horizontal coefficients only. 

Horizontal coefficients are determined for each series of flow measure- 
ments at a given hydraulic section. Flow measurements are usually made in 
groups of about 20 over a period of about 10 days, and a new set of horizontal 
coefficients is derived for each such group of measurements from the 0.4 
depth velocities observed during the measurements. The average of the ve- 
locities observed at 0.4 depth in each panel is plotted as an ordinate against 
the horizontal location of the measuring point as an abscissa. A smooth curve 
is drawn through the points so plotted for all the panels. This curve is con- 
sidered to show the horizontal distribution of current velocities at 0.4 depth 
across the hydraulic section, and the horizontal coefficient for each panel is 
the ratio of the average ordinate of the curve over the width of the panel to the 
average of the observed 0.4 depth velocities. 

Most Lake Survey flow measurements are made from a catamaran. Some- 
times special conditions require that measurements be made from a bridge or 
a cableway. In measuring froma catamaran, velocity observations at each 
panel at 0.4 depth are made simultaneously with three current meters. The 
three meters are suspended on the line of the hydraulic section by cables, with 
a separation of 3 ft between adjacent meters. 


COMPARATIVE TEST OF CUP AND SCREW METERS 


A group of flow measurements was made by the Lake Survey in the fall of 
1958 at the Stella Niagara hydraulic section in the lower Niagara River. The 
metering was done from a catamaran, simultaneous velocities for 5-min peri- 
ods being observed with three current meters suspended at 0.4 depth in each 
panel. This provided a good opportunity to compare velocities obtained by dif- 
ferent types of meters, and in all of the 20 flow measurements made over a 
period of 16 days, one cup and two screw meters were used simultaneously. 

The cup meter used was the Price meter. Screw meters used were the 
Neyrpic, which is manufactured in France, and the Ott, which is manufactured 
in the Federal Republic of Germany. These three meters are shown in Fig. 2. 

Three meters of each make were used for the measurements. Throughout 
a single flow measurement, the same three meters, one of each make, were 
operated in the same relative positions. In the successive measurements, a 
different meter of each make was used then in the preceding measurement and 
the three makes were placed in different relative positions so that about the 
same amount of data was obtained for each make in a given position. 

Physical characteristics of the Stella Niagara hydraulic section are showr 
by the bottom profile in Fig. 3. The section is seen to be broad, shallow, anc 
uniform. Also given in Fig. 3 is a curve showing the horizontal distribution o! 
current velocities at 0.4 depth. In measuring the flow past the section, 
velocities were metered successively at the locations in Fig. 3 numbered J 
through 10 and labeled Panel Points. _ ; 


ANALYSIS OF DATA oh ae 3 


Table 1 shows the average velocities obtained at each metering locatio1 
for each make of meter for the entire group of 20 flow measurements, and the 
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average velocities for each make of meter for all metering locations and all 


measurements. ; 
The horizontal distribution of velocity curve in Fig. 3 indicates that simul- 
taneous velocities recorded by the three meters may not be the same because 


(a) HORIZONTAL DISTRIBUTION OF VELOCITY 


Velocity, in feet per second 


I 
Panel point 


United States shore 


2 
fo} 
= 
“” 
inq 
Ss 
Bo) 
oO 
& 
oO 
oO 


750 1000 1250 1500 
Distance from Canadian shore, in feet | 


Depth, in feet 


FIG. 3.—STELLA NIAGARA HYDRAULIC SECTION 


of the slightly different positions of the meters on the line of the hydraulic 
section, particularly at panel points 1 and 10. Velocities at each panel point, 
compared with respect to position instead of type of meter, should reflect the 


TABLE 1.—AVERAGE VELOCITIES IN fps OBTAINED FROM 
TWENTY 5-MIN METERINGS 


Panel Point 


CABANOULWNE 


shape of the horizontal distribution of velocity curve. Average velocities ob-_ 
tained from the twenty 5-min meterings at each of the three meter positions — 
left, center, and right—looking upstream at panel points 1 and 10 were com- 
puted, and the velocity differences between the left meter and center meter 
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und between the right meter and center meter were compared with corres- 
bonding differences indicated by the curve. The data are given in Table 2. 

_ Although the data in Table 2 show that relative position has a definite ef- 
leet on measured velocity, this phenomenon does not appreciably affect the 
lata of Table 1 because of the interchange of meters between successive 
measurements. 


TABLE 2.—COMPARISON OF VELOCITIES BASED ON METER POSITION 


Mean Velocity Differences in Feet per Second 
Indicated Meter Position Minus Center Position 
Obtained from 


Meter 
Position 


Twenty 5 min meterings Velocity curve 


Comparisons of individual meterings of the three makes of meters, taken in 
airs, were made. On Fig. 4 these comparisons are plotted at panel point 4) 
sing the first and second minutes and third and fourth minutes of each 5-min 
aetering as individual meterings for meter combinations Price-Ott, Price- 
leyrpic, and Neyrpic-Ott. The points so plotted are in every case well cen- 
sred around the line of equal velocity throughout the range of measured 


TABLE 3. COMPARISON OF AVERAGE MEASURED FLOWS 


OBTAINED BY DIFFERENT METERS 
Fee! 2 ee ee ee eee 


Percentage Variation of Flow 
Obtained by Single Meter 
from Flow Obtained 
by All Meters 


Average Flow, 
in Cubic Feet per Second, for 
Measurements Computed Using 
Data from Indicated Meter Only 


Meter 


rice 


191,850 -9.13 
Itt 192,420 +0.17 
leyrpic 192,440 +0,18 
ll meters 192,100 


| a Some scatter of the points should be expected, but it is small, the 
oints lying almost entirely within a 5%, deviation from the equal velocity line. 

A further comparison was made in terms of the average measured flows | 
r all measurements using data from each make of meter exclusively, as 


in Table 3. 


HY 4 


April, 1960 


18 


S LNIOd TANVd ‘NOILLOUS VUVDVIN VITALS ‘SGOIMaAd 


puodsas sad jaa} ul ‘Ay D0;aA 
Se Gacy SE GUGy GE Ge 


LLO *SA OIdYAAN (9) ; DIdHAIN “SA JOIYd (9) 110 “SA 391Yd () 


Puodas Jad jaay us ‘Ajo0ja, 


HY 4 STREAM VELOCITY METERS 19 
CONCLUSIONS 


Based on the data presented herein, it is concluded that no greater ac- 
‘curacy can be achieved in Lake Survey flow measurements from a catamaran 
‘by using screw type meters instead of cup type current meters. 
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EFFECTS OF FLOOD FLOW ON CHANNEL BOUNDARIES? 


By D. A. Parsons! 


SYNOPSIS 
| Some general and particular things of interest, learned in studies of 
)stream bank stabilization methods in Western and Central New York State by 
ithe Agricultural Research Service, are given. The studies, for the most part, 
consist of attempts to relate observed effects of floods as evidenced by dam- 
ages to stream bank revetment, to the channel geometry and the qualities of 
the flood flows. 


| INTRODUCTION 


The Soil and Water Conservation Research Division, Agricultural Re- 
search Service, has been giving attention to stream bank stabilization prob- 
lems in Western and Central New York State, particularly to those in the 
: ee River watershed.2 Although the studies have far to go to meet the ob- 
jectives of a good understanding of flood flows and economical means for 


Note.—Discussion open until September 1, 1960. To extend the closing date one 
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is part of the copyrighted Journal of the Hydraulics Division, Proceedings of the Ameri- 
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_ 2 4streambank Stabilization,” by Leon F. Silberberger. Agricultural Engineering, 
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stabilization of channels, some general and some particular things have been 
learned that are of interest. Most of the observations have to do with patterns 
of damage to stabilized channels. Of interest, too, may be the measurements, 
the thinking, and the methods by which attempts are being made to relate the 
damages, or instability, of the channel bank to the flood flows that occurred. 

From the standpoint of evaluation and development of streambank-protection 
methods, it seems proper to separate the elements involved: the stream and 
its qualities that determine the potential for bank erosion, and the boundaries 
and their qualities of resistance to the attack of the stream. These two things 
are not completely independent, because in natural streams the flows create 
the channel and conversely the channel and flood plain geometry greatly influ- 
ence the flow qualities. Yet, how else can a rational design for channel sta- 
bility be achieved than by matching the resistance of the boundary against the 
intensity of the attack of the flow? This requires quantitative measures of 
both. 

Elements of the boundary are moved (eroded) because of the forces exerted 
upon them. It is common knowledge that these forces are primarily dependent 
on the velocity of the fluid in the immediate vicinity of the protruding portions 
of the elements. It seems then that to evaluate the potential of the flood flows 
to erode the boundary or destroy the bank lining, it is necessary to measure 
the water speeds near the boundary. This, of course, can be done, but the 
question arises, how close should the measurement be made? The speed of 
the water far removed from the boundary is not of direct concern. 

Perhaps the fluid shear stress at the boundary is a preferable measure. 
Its dimension is force per unit area and, although shear stress does not seem 
to be very descriptive of the actual forces that act upon the elements of the 
channel boundary, its dimension seems proper. But here, too, there are dif- 
ficulties. It has not yet been found easy to measure the boundary stress from 
spot to spot. It seems, too, that the boundary drag is a good measure of the 
intensity of attack of the stream only to the extent that it reflects the forces 
that would act upon a protruding element. 

Briefly, the work in New York is directed, among other things, toward the 
determination of the potential of flood flows to destroy the channel boundaries 
from spot to spot for many geometrical situations of channel alinement and 
bank height. Quantitative measures of some type are required. These may be 
the boundary shear stress, but thus far, reliance has been on comparisons of 
the flood-flow qualities of mean velocity head and mean tractive force for the 
channel with the measured sizes of riprap stones for the condition of incipient 
movement. The use of the measurable, mean qualities of the channel flow is 
not only helpful in investigations; it is quite essential to application in design 
of any relationships that are developed. 

Coarse materials transported by the flood flows are factors in the process. 
They may be capable of materially altering the flood flows and flood stages, as 
with ice jams, trees blocking bridge openings, etc., but, for the most part, 
their role in the bank-erosion process appears to be principally as agents o1 
helpers of the flow in the attack on the bank. Y 

In the case of ice, it is quite usual to find the motion of the pack ice, fol. 
lowing an ice’ jam, confined solely to the channel with vertical shear planes a 
about the toe of the bank. The ice along the bank thus acts as a buffer for the 
bank against the potentially destructive floes within the channel. However, i 
has been observed that this buffer ice is often insufficient to withstand the fol. 
lowing onslaught of the stream in those places along the boundary normall: 
subject to the heaviest action of the flood flows. 
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Heavy objects like rooted trees, stumps, and boulders may drag or catch 
on the channel bottom, causing higher water speeds along the banks, especially 
for the smaller streams, than would normally prevail. The materials carried 
nigher in the flow appear to be true agents of the flow, augmenting the ability 
of the water in its action. 

_ The manner of the erosive action of transported materials is varied. The 
processes include impact of moving objects with the bank lining, an abrasion, 
and sometimes a lodgment against a protruding element of the lining, thereby 
increasing the force on the element tending to move it out of place. The ef- 
fectiveness of these processes is dependent on the magnitudes of the water 
speeds in the immediate vicinity of the boundary and the transported object. 


: STREAM CONDITIONS 
: The usual stream conditions in Buffalo Creek, N. Y., range from a width of 
about 65 ft, slope 0.006, probable annual peak flow 2,500 cfs above Java 
Village to a width of about 120 ft, slope 0.0025, and probable annual peak flow 
ia cfs at the lower end. The bed material is sand, gravel, and cobbles, 

redominantly gravel. The stream has a recent history of degradation in some 
places, particularly in the upper reaches, and flows on nearly smooth bed 
rock in some places, especially toward its lower end. Floods are numerous, 
pecurring mainly in the winter and early spring. Ice conditions are some- 
times severe in the winter floods. 


: 
: GENERAL OBSERVATIONS 


Most streams are continuously changing their channels by aggradation or 
egradation of the beds and by the building up and tearing down of the banks. 
is is inherent in natural streams. Therefore, prior to initiation of any 
ajor channel work or change in flow and bed material transport regimen, 
f is need to determine the current trends of channel change, especially 
any progressive changes in channel-bed elevation. It is also essential to pre- 
erraine the effects of the installed measures on future stream-bed eleva- 
tions, in addition to the lateral migrations of the channel, before the plan for 
work can be truly called a design for improvements. Both benefits and dam- 
ages are partly dependent on the future elevation of the channel bed; so also, 
.s the exact design of most measures for stream-bank protection. Some highly 
undesirable channel conditions have been observed that were due to the neglect 
. one or both of these work-prerequisites. 

It should not be necessary to prove the obvious, that is, that the potential 
nf flood flow to erode the channel boundary varies from spot to spot. It is de- 
irable though to restate some elemental facts that tend to show this. They 
annot be over-emphasized, either from the standpoint of channel investiga- 
ions and evaluations of protective measures, or from the standpoint of actual 
ank-stabilization work. They have to do with meandering channels. Pending 
he time when studies of secondary currents, helical flows and variations in 
ractive forces over the wetted portions of the channel provide the means for 
etermining exceptions to these simple rules, they may not be ignored with 


mpunity. 


e 
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C. R. Allen presented a convincing description of the mystic tendency fc 
the meander pattern to move down valley. He said, in part: 


“In 1892 the writer meandered forty miles of the Des Moines River and 
plotted the same in relation to the U. S. Survey made in 1847. I found, 
which I presume is not new, but a well-known law, that whenever the 
reaches ran parallel with the general course of the river there was lit- 
tle change at the shore line by erosion; that wherever the reaches ran 
transversely to the general course... the banks have changed by ero- 
sion downstream.” 


This tendency for the down-valley migration of the meander pattern ha 
been observed many times and in many places since then, and when considere 
with the erosion of the outside bank in bends, obvious to all, it must be con 
cluded that the potential of the stream to erode its banks not only varies fro! 
one spot to another, but is persistently high in some places. For stabilit 
these places along the bank must be provided with a more resistant materi: 
than usually occurs in nature. It is just as important, too, to realize that ther 
is no need for protective work along the opposite shore. The opposite shor 
tends to build up by deposition. 

The potential of the flood waters to erode the bank also varies up and dow 
the bank slope, being less in the shallow water near the water’s edge tha 
deeper in the flow. A better knowledge of its variation from the toe of the ban 
to the water surface at stream edge would be of aid in revetment design. Muc 
of the present knowledge of this variation for straight channels has been sum 
marized by E. W. Lane.4 

The intensity of attack of the flood waters at a spot along the boundar 
generally increases with the rarity of the flood, and varies with time fo 
steady flows because of the turbulence and large-scale eddies of the floo 
waters. The chief concern in design is therefore with that greatest flood flo 
associated with the allowable risk of exceedance within the design life of th 
stabilization work. However, lesser flows may produce more severe condi 
tions at some places along the bank and should be considered. 

It has been found that the pertinent mean-flow qualities of the flood water 
within and adjacent to the channel are somewhat different for the larger flood 
than might be deduced by extrapolation of values from smaller flows com 
pletely confined to the channel. The topographical features of the valley floc 
and walls that tend to restrict the flow, including highway fills and artificie 
structures, become important. Sharp channel bends, immediately downstreaz 
from the bank areas under consideration, also become important in the large 
floods. High energy losses in the larger flows through these bends tend to re 
duce the energy gradient immediately upstream. 

Some particular and generally localized effects of floods have been ob 
served. For the most part, these observed effects are the complete or parti 
failure of established protective measures. These spots are almost entire! 
GAs nod Mies fec/Md factiak emne al nie oe 

3 “Protection of River Banks at Ottumwa, Iowa,” by C. R. Allen. Iowa Civil Engi 


neers and Surveyors Society, Proceedings of the Seventh Annual Convention, January 
1895, pp. 39-46. a 


eke of Stable Channels,” by E. W. Lane, Trahsactions, ASCE, Vol. 120, 1955 
p. ; 
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confined to the down-valley bank and in places where the boundary currents 
are strong. The observations indicate the following: 


1. There is a portion of the bend wherein the water is deepest and the at- 
tack of the stream greater than normal. 
_ 2. Great difficulty is experienced in holding the protective works along the 
down-valley bank in those places where the banks are low. High-speed water 
leaving the channel creates a condition that is more severe than usual. It has 
been found helpful in some places to build the bank up above the flood stage 
before revetment. 
: 3. Out-of-bank flood waters, returning in a concentrated fashion over the 
bank to the main channel, are often able to damage the protective works if it is 
mot somehow re-enforced. 
4. Inadequate distance between revetted banks, constituting a restriction in 
the channel width, increases the severity of the action of the flood flows on 
both banks. In a meandering channel this could ordinarily occur only near the 
downstream end of the revetment on one bank and the beginning of the revet- 
ment on the opposite bank. 
5. Protrusions into the flow are subject to high forces. 
6. Irregularities in the alinement of the bank induce severe action by the 
stream on those exposed portions. 
7. Failure of the bank-stabilization work is almost always in the nature of 
an erosion, that is, a little bit ata time. The damage to, or loss of, one small 
element of the lining may occur almost instantaneously; but several units of 
ime,a whole flood, or a succession of floods may be required to disintegrate 
e lining to the extent that the bank stability is lost. In fact, this process of 
sintegration is tacitly recognized in one method of stabilization. Periodic 
eeding of semi-resistant materials to an eroding area can effectively control 
e migration of the shore or bank. 
8. An effect opposite to the preceding is the continued progression and 
wild up in the down-valley direction (Fig. 1) of the shore after stabilization - 
of the opposite eroding bank of a rapidly migrating stream. A particular point 
of interest in this is the extent in plan of this deposition at the beginning of 
bends. 
i 9. The fluctuating boundary currents with associated pressure variations, 
aided. by surges, waves, seiches, and sometimes ground-water flow to the 
stream, tend to cause loss of the finer bank materials from beneath the pro- 
tective veneer along the face of the bank. Failure in this fashion of the sup- 
porting base of the lining has frequently caused complete failure of otherwise 
excellent revetment. 
' 


PATTERNS OF FLOOD EFFECTS IN BENDS 


Fig. 1 illustrates a typical stream-channel situation and a manner of 

tudy of stream behavior. A reference line, A-B, is drawn or sighted along 
ie eroding down-valley bank. It is not surprising to find that the alinement of 

iis bank is the key to the bed geometry and to the erosive forces in the chan- 
el bend immediately downstream. Instability of the bank, represented by the 
ine A-B, with consequent down-valley migration, immediately alters the sit- 
lation in the bend below. The design alinement and stabilization of this bank 


ao 


should therefore precede those of the one downstream. 
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Flood waters moving along parallel to line A-B tend to continue in the same 
direction, but as they approach the down-stream bend a force is exerted upon 
he flow, tending to change its direction. The slowly moving water near the 
bed and the banks is more readily turned. The fast water near the surface, 
-elatively free of suspended bed materials, crosses over to be in close prox- 
imity with the outside bank, somewhat as indicated by the arrows. 

Some of the observations of partial or complete revetment failure in bends 
are portrayed in Fig. 2. Also shown are locations of several deposition points 
n the upstream parts of bends following bank stabilization work, and the lo- 
ations of deep points in bends of unstabilized channels. Horizontal distances 
ave been divided by stream width to give a dimensionless representation. It 
is not expected to find true similitude in nature; yet, considering the order of 

nitude of the dispersion of the phenomena with the possible error inherent 
n the assumption, along with the great utility indimensionless representation, 
polarity of streams of different size is assumed. 

Even stream width is a quality of large dispersion, especially for migrat- 
ng streams. For rapidly migrating streams there is no well defined up-valley 
bank, and one needs to search for sites to measure and average. 

_ The position of the deposition point C in Fig. 1 is an important considera- 
ion since it is logically associated with the beginning point of need for a 
fugged type of revetment. A rapidly migrating stream would leave the bank 
m this area in a raw condition. This superficially indicates the need for 
strong revetment much further upstream than is truly the case. Common mis- 
udgments in stream-bank-protection work are to revet this bank too far up- 
itream and fail to go far enough downstream into the bend on the opposite 
ide. This is a result of the otherwise good, simple rule of reveting the bank 
hat is eroding. 

The true significance of the observed location of the downstream deposition 
imit at the beginning of bends has not-been established. Referring to Fig. 1, 
e condition seems to be that the stream flows freely in a curved path as it 
asses from one restraining bank to the other on the opposite side. Deposition 
rogresses down valley to conform with the stream path. Its point of ending, 
1 presumably would depend upon the exact location of the previously 
roding bank in addition to the curved alinement that the stream adopts. The 
xtent of the revetment on the opposite side is also a factor. The distance 
esignated Ar between the point C and arc with radius r was close to and 
veraged 0.04 W for 10 out of 12 measurements. 

The inference of the plotted representation of damages in bends is that the 
rosive forces of the stream begin to become severe at the point B, reach a 
maximum in severity at a point on the order of one stream width away from 
ne reference line, depending somewhat on the bend radius and approach con- 
itions, and then moderate further along in the bend. The variation from the 
verage picture is considerable. There are several reasons why this should 
eso. The tacit assumption of uniform revetment strength throughout each of 
1e bends is undoubtedly untrue. The approach conditions varied from bend to 
nd. Also, the adopted position of the reference line may not have been the 
fective one in every case. oe 

The picture of damages to protective works in Fig. 2 may be of aid by in- 
tet in a rough way, the places in bends that are particularly difficult to 


stect. But it leaves much to be desired, if one is faced with the need to 
yecify precisely the protective works to be installed. 


: 
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Whereas data are not adequate to indicate with confidence the strength of 
the required revetments, knowledge is not completely lacking. The question 
is whether to wait for additional measurements to augment the very limited 
information that is available, or use it as a basis for a first guess. Since the 
required data are difficultly and slowly obtained, and since the need for quanti- 
tative estimates are current, the guess will be made. The explanation of the 
estimate will also show the methods that seem to be required to derive the 
needed knowledge. Actually, in making these studies it seems necessary to 
set up trial or tentative relationships such as Eq. 3 (to be presented subse- 
quently) and make some preliminary assumptions in regard to the pertinent 
flow parameters and relative resistances. If the observations do not conform 
with the assumptions, new schemes are tried. 

There is need first to explain the origin of coordinates of Fig. 2. Each of 
the deep points and damage areas were located on a plan of the stream chan- 
nel, similar to that of Fig. 1. A line was then drawn tangent to the eroding 
bank at a point one stream width below the reference line. The angle between 
this tangent line and the reference line is a9. It was then assumed that the 


path of the stream as it entered the bend was that of a circular curve such that 


ee ee ee ee (1) 


r 
w  1-cos a 
The origin of coordinates for Fig. 2 corresponds to the point on Fig. ] 
marked P.C. If the tangent line is extended to intersect a line one stream 
width up valley from A-B, the distance from the intersection to the point P.C., 
the origin of coordinates, is (r/w) tan (a,/2). The Y and Y/W direction in 
Fig. 2 is the same as that of the reference line A-B in Fig. 1. There un- 
doubtedly is a better way of representing the observations. The angle ag may 
not be the best angle to select anda curve other than circular may be more 
representative. However, a better manner of representation of the benc 
geometry should not be greatly different for stream conditions like those ob- 
served. 
The straight line, 


Y: Xx + 
a ee ee (2 


in Fig. 2 is drawn to represent the positions of maximum intensity of attacl 
in bends. The requisite strength of revetment along this line, as determinec 
from field experience, is© 


“or 


Re diet ee 
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. 


ing of the ASCE and is believed to be much better for the very large and very sm 


‘ 
. . + : : . : 
© Eq. 3 is a revision of a simpler one given at the May, 1959 Cleveland, Ohio nes 
values of (¥/W)p. ’ 
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Distributions: 
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In Eqs. 3 and4 P'' and P are the minimum requisite strength, expressec 
as an equivalent median riprap stone diameter in inches or feet respectively, 
on a 1 or 2 bank slope, placed with only moderate care, and subjected tc 
stream conditions as found in western and central New York State; the sub- 
script p indicates the coordinate to the point of maximum attack in bends; anc 
Z is a flow parameter descriptive of the within-channel flood flow, either the 
velocity head (V2/2 g) in feet or the mean tractive force (To) in pounds pe} 
square foot. It happens to be that for the damaging flood flows that have bee! 
measured in Buffalo Creek, the average values of these two parameters aré 
about equal. However, this is not a general condition, and use of the velocity 
head is favored at the present time. Lines of equal required strength art 
presumed to be positioned about the line of the maximum, somewhat as is the 
envelope curve for the bulk of the damage data. Fig. 3 shows values of (r/w) 
(x/w)p, (Y/W)p and (P/Z) as they vary with ap. 

Eqs. 3 and 4 are strictly empirical and were chosen, because of their sim. 
plicity, from many others, that might have been used to represent the limitec 
information. The constant, 4 in Eq. 3 was chosen simply as a reasonable 
figure, but it gives results quite close, using the velocity head for Z, to th 
values for gravel as given® by S. Fortier and F. C. Scobey for straight chan 
nels. The line marked P/Z in Fig. 3 may be used in lieu of Eq. 4. 

Most of the damage data were for quarried stone revetment with a speci: 
fied size of approximately D = 34 p where D is the size in inches and p is thi 
proportion by weight that is smaller. Thus D5 = P''=17in. The flood flow: 
that caused the damage were measured in only a part of the cases. This wa 
done by means of post-flood surveys of peak water-surface elevations ani 
channel cross sections from which peak discharges, mean velocities, slopes 
and mean boundary stresses were obtained. The values of Tp and v2/ 2 ¢ thu 
determined were, of course, variable, but averaged about 3/2. The envelop 
curve for damages crosses the straight line at (Y/W)p = §.3. Using thes 
values, Eq. 3 gives 16 in. for the median stone size. 

A few other particular cases of the same and of different bank linings wer 
studied in the same fashion with reasonably good agreement. In one cas 
where (¥/W)p = 8.4, there was a vegetative lining of a combination of excel 
lent 3-yr growth of grass and woody vegetation including basket willow. It 
resistance was apparently exceeded for about 40 ft (0.6 W) along the outsid 
of a long-radius bend on Bennettsville Creek, N. Y. The stream survey an 
testimony of the land owner resulted in an estimate of 1lb per sq ft for the T 
and 0.55 for V2/2 g for the flood that caused the damage. Eq. 3 yields 9 i: 
and 5 in. respectively for the equivalent median stone size as a measure ¢ 
the maximum intensity of attack. For dumped, loose stone, these values woul 
be increased to 11 in. and 6 in., respectively. This compares with a maximul 
of 6 in. for the resistance of a good stand of long, green bermuda grass a 
determined from the data® of W. O. Ree and V. J. Palmer and the results « 
flume tests for critical tractive force. E 


x 
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5 “Permissible Canal Velocities,” by S. Fortier and F.C. Scobey. Transaction 
ASCE, Vol. 89, 1926, p. 940. 3 
6 “Flow of Water in Channels Protected by Vegetative Linings,” by W. O. Ree a 


Bee - Palmer. U.S. Department of Agriculture, Technical Bulletin No. 967, Februar 
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At the other end of the line it has been observed that stones approaching 
1 ton in weight (about 35 in. equivalent diameter) have been moved under ex- 
treme conditions along Buffalo Creek, N. Y., when partially isolated by loss 
of the finer stones or laid in an exposed manner. In a sharp bend with, say, 
(Y/W)p = 1.5 and To or v2/2 g = 1.8 representing an extreme condition in 
Buffalo Creek, Eq. 3 for median stones gives 33 in. There are no measure- 
ments to check on the validity of this estimate of the median size required for 
stability, but from experience with the creek it seems to be a reasonable 
value. 

A bend in Page Brook, N. Y., on the Verner Lewis farm, revetted with 
dumped riprap stones of 7-1/2 in. median diameter, lost a few stones ina 
length of about 0.2 W (7 ft). The (Y/W)p as defined was 4.54. Estimated 
values, based on the stream survey, for T, and v2/2 g were 1.00 and 0.49, 
respectively. Eq. 3 yields 11.7 in. using Ty and 5.7 in. using the velocity 
head. Here again the computed sizes are for stones placed with moderate 
care. One of the assumptions currently being made is that dumped stones need 
to be about 20% larger, and stones placed with much care to eliminate protru- 
sions from the face of the blanket about 20% smaller. Then the estimates for 
dumped stones in this bend in Page Brook for the flood experienced is 14 in. 


and 6.8 in., respectively. The bank slope in this bend was 1 on 1 ‘, a little 
steeper than is applicable with Eq. 3. 

A riprapped bend on the E. Foss farm along Buffalo Creek experienced two 
almost identical floods of about 25-yr expectancy. Post-flood surveys indi- 
cated flows of about 10,000 cfs and 0.5 ft per sec mean velocity. The (Y/W)p 
is 2.96, T,=1.5 and v2/2 g = 1.44. The specified median stone size was 17 in. 
Eq. 3 yields 21 in. and 20.3 in., respectively. Since this was contract work 
and much care was used in placement of the rock, the requisite size would be 
somewhat less. Values 20% less would be 16.8 in. and 16.3 in. The particular 
interest in this bend arose because the intensity of flood-water attack and the 
revetment strength appeared to be almost perfectly matched. Some of the 
smaller stones on the face of the bank were washed away and several of the 


large ones had been shifted slightly, but after two major flood experiences, 
the revetment was still in good condition. 


SOME ADDITIONAL ASPECTS OF REVETMENT DESIGN 


The stability of protective linings appears to be considerably affected by 
the extent to which the elemental parts jut into the flow. Although no quanti - 
tative measurements were made, there is evidence to suggest this. The 
helter-skelter arrangement of the stones of dumped riprap have appeared t 
be less effective than when the lining was placed by crane and clam shell with 
a moderate attempt at elimination of irregularities on the face of the Sacell 
Work done with a considerable effort to fit the stones together, chink the holes, 
and remove all large protrusions has generally stood up still better. Impres 
Sive, too, in this respect, is the comparative stability of closely fitted experi- 
mental concrete revetment blocks. These blocks were 24 in. by 16 in. by 4 in. 
each having twenty four 2 in. by 2 in. by 4 in. holes and weighing 83 lb. 
of about six hundred on the bank of a bend, in a spot that the formula sugges 
the need for 26 in. stones (800 lb) two blocks at exposed corners were lost t 
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r was lost, both upstream and downstream of the experiment. The relatively 
od showing of the concrete blocks in this respect is believed due to the fact 
at they provided little opportunity for the water and transported materials 
impinge upon their upstream edges. 

Irregularities in stabilized bank alinement that appear to induce a more in- 
use attack by the stream are of several types. A poorly defined notion exists 
at the stream tends to follow a smoothly curving path of decreasing curva- 
re throughout a bend, and that banks encroaching on this path are more 
Inerable than normal. Specifically, in several cases where the curvature of 
> downstream portion of a bend is greater than that of the apparent free-flow 
rvature of the stream as it enters the bend, the lower portion of the bank in 
2 bend is not only subjected to heavy action, but its presence seems to re- 
‘ve somewhat the intensity of the expected action immediately upstream. 
‘esumably, a large loss of energy in the downstream portion during large 
ods causes a reduction in the energy gradient upstream, with consequent 
duction in water speeds. 

The downstream shoulder of the bank that is created by a depression in 
> revetment for a stream crossing for farm vehicles or in the mouth of a 
1all tributary is particularly vulnerable. 

In the one case observed and shown on Fig. 2 where the upstream end of the 
vetment in a bend appears to have been extended out into the path.of the 
ream, severe action was experienced. If the downstream end of revetment 
es not curve down valley, somewhat in conformance with the path of the 
ww, it will in effect be a protrusion into the flow and need to be stronger than 
lerwise. 

The loss of fine materials from beneath a portion of the revetment has 
en experienced with the experimental concrete blocks. This was not evident 
til a flood eroded away the adjoining stone riprap revetment and the bank 
wterial along the edges of the experimental mat. Fine material from be- 
ath the blocks was lost adjacent to the exposed edges. Little, if any, was 
st from the small holes in the blocks. The mat appears to be stable because 
ich of it was laid on a layer of sand, gravel, and cobbles. The coarser ma- 
‘ial of the foundation was not moved out by the flood. The experience with 
S one flood, however, appears to demonstrate one mode of failure of mat 
es of lining. 


BOUNDARY SHEAR STRESS EVALUATIONS 


Measurements of the depths of scour of materials of known size and density 
ym some of the 2 in. by 2 in. holes in the experimental revetment have been 
de. This was an attempt to delineate the variation in the intensity of action 
the flood flows against the bank. Estimates of the boundary stresses have 
; made, but the technique of analysis needs verification. Also, the scour- 
data were very erratic for reasons only partly understood. For these 
Ens, they will not be reported here, except to say the depths of scour were 
initely greater in the larger floods; the scour was greater deep in the flow 
n where the water was shallow, and the mean scour differed for different 
ns along the bank. Good results are anticipated along this line from the 
‘rent studies of flows in laboratory channels at the Hydrodynamics Labora- 
Massachusetts Institute of Technology, Cambridge, Mass. A comparison 
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of the boundary shear stress variation throughout a 60° bend as reported? b 
R. E. Nece, C. A. Givler and P. A. Drinker with the pattern of flood effects i 
natural stream channels is encouraging. 


MAINTENANCE CONSIDERATION 


A discussion of stream bank stabilization work and flood effects would b 
incomplete without some mention of maintenance. A provision for some kin 
of maintenance effort is an absolute must for stream bank stabilization wor 
to remain effective. There are so many uncertainties about design and chanc 
happenings that may affect the stability that constant watching and prompt ac 
tion are a necessity. 

Actually, the expected quality and timeliness of the maintenance effort is 
factor in design. One farmer was observed to be doing an excellent job « 
maintaining stable banks of a small stream. The riprap stones that he use 
were from his own fields and were too small for the toughest spots in th 
bends. He made the stone blanket thicker in these spots. If a large floo 
washed away some of the material he would promptly haul in another loac 
The strength, extent, and quality of initial stabilization work should be in 
creased as the level of expected attention to maintenance goes down. 


SOME OTHER FACTORS INVOLVED 


The principal direction of the work in New York on this problem has bee 
toward the flow element of the process, because it seems to be the most neg 
lected and toughest portion. Needless to say, there is much yet to be done 
For the most port, too, the emphasis has been on complete stabilizatio1 
rather than on methods for reduction of stream-bank erosion to tolerab] 
rates, although an understanding of the flow capabilities should help in this 
The various methods for complete stabilization have hardly been touche 
upon. Cost is paramount in this field, and each method has peculiar probler 
of its own. The highly important questions of changes in streambed elevatic 
have been given only minor attention. 

The use of vegetation in streambank stabilization has been given muc 
study inthe Buffalo Creek work. Real progress has been made in this fiel 
The observations and the conclusions in regard to the place of vegetation hav 
been ably expressed8 by Harry Porter and Leon Silberberger, M. ASCE. | 
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FRICTION LOSSES IN LINES WITH SERVICE CONNECTIONS 


By David L. Muss,! M. ASCE 


SYNOPSIS 


The determination of friction losses in water mains is complicated by ser- 
ice connections along the line which result in different flows in each section 
f the pipe. One method of determining these losses2 has been to assume that 
1¢ amount actually withdrawn along the line is withdrawn as a single quantity 
t the end of the line. This assumption may result in the computation of fric- 
on losses which may be as much as 2.85 times the correct value, depending 
nthe number of connections and the ratio of the ‘quantity withdrawn to the total 
let flow. European practice has been to assume that 45% of the amount with- 
rawn is removedat the beginning of the section and 55% at its end.3,4 Although 
is approximation results in a smallerror(-5.7% to +3.5%) for alarge number 
f connections, the value of the friction loss computed for a small number of 
mnnections may be as little as 35% of the correct value. 

Two alternate methods of approach are proposed. Both are simple and ac- 
irate. The first method, applicable only to dead-end systems, or other situa- 
ons where the flows are known, involves the selection of either a correction 
ctor, to be applied to the flow entering the line, or tothe friction losses com- 
ited on the basis of the inlet flow or to the length of pipe. Factors vary from 
568 to 1.000 for equivalent flow and from 0.351 to 1.000 for friction loss or 


ee 
Note.—Discussion open until September 1, 1960. To extend the closing date one 
onth, a written request must be filed with the Executive Secretary, ASCE. This paper 
part of the copyrighted Journal of the Hydraulics Division, Proceedings of the 
nerican Society of Civil Engineers, Vol. 86, No. HY 4, April, 1960. 

1 Cons. Engr., Asst. Prof, of Civ. Engrg., City College of New York, N.Y. 

2 “Water Supply and Purification,” by W. A. Hardenbergh, Internat. Textbook Co. 
ranton, Pa., 1952. 

3 “Idralica,” by Umberto Puppini, Nicola Zanichelli, Bologna, Italy, 1947. 

4 “Vodosnabjenie,” by N. N. Abramov, N. N, Geniev and V. I. Paplov (Water Supply 


gineering), Moscow, Russia, 1958. 
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equivalent length. This method is precise and there is no error in its applica 
tion. 

The second method is applicable to either the dead end or loop system, ar 
the flows need not be known in advance. A portion of the flow, of the amou: 
withdrawn, is assumed to be withdrawn at the beginning of the line and tk 
balance from the end of the line. This portion varies from 0 to 44% dependir 
on the number of connections. This method is approximate, but the range « 
error is small and is limited to from -3.5% to +5.5%. 


INTRODUCTION 


The determination of friction losses in water mains is complicated by tl 
service connections which remove water from the main along its path. Th 
results in a different value of flow from each service connection. A conservé 
tive simplification of this problem isto assume that the quantity removed alor 
a given length of main is instead removed as a single quantity at the endof tl 
main as shown in Fig. 1. 


FIG, 1.—CONSERVATIVE SIMPLIFICATION 


That this assumption is conservative is illustrated in Table 1, which ind 
cates that the friction losses so computed for a line may vary from 0% to : 
much as 185% greater than that which actually exists for the condition of se: 
vice connections uniformly spaced along the line from its inlet to its outlet ; 


TABLE 1,—PERCENTAGE ERROR RESULTING FROM ASSUMPTION THAT 
QUANTITY WITHDRAWN IS REMOVED FROM OUTLET END 


Ratio of Quantity Withdrawn (q) to Total Inlet Flow (Q) 


computed by the summation of losses between each connection by rigorous cor 
putation in the manner shown in Fig. 3 (to be presented|subsequently). In tt 
discussion it has been assumed that (a) the Hazen-Williams flow equation, a 
plies, wherein the friction loss is assumed to be a function of the 1.85 pow 
of the discharge; (b) the connections are uniformly spaced along the main; a 
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¢) the withdrawals at each service connection are equal. Turbulence losses 
t each of the connections have been neglected. 


EQUIVALENT DISCHARGE AND LENGTH 


While conservation is undoubtedly desirable in water Supply engineering, 
alues so obtained may be unrealistic. One alternate approach might be (a) 
he determination of an equivalent flow (Qe) which would give the same friction 
oss as that which is actually encountered; or (b) the determination of a cor- 
ection factor (K'), which might be applied to the friction loss computed in the 
onventional manner or to the length of the line to determine the true loss. In 
his latter case, the equivalent length of the pipe would be equal to K" L. 


INFINITE NUMBER OF SERVICE CONNECTIONS 


With an infinite number of connections, the flow at any point x in a line of 
sngth L is equal to Qt, the total incoming flow, less a proportion, x/L, of the 


FIG, 2.—INFINITE NUMBER OF CONNECTIONS 


iantity withdrawn q—or Q; - (x/L)q. This is illustrated in Fig. 2. If the fric- 
on loss, hg, is a function of Q1.85, then: 
L 
1.85 _ : 
he=KLQ, -f K dx (@, 
0 


a =K LK' afeae » +++ (1) 


‘mad ita 


There Qe is the equivalent discharge in the pipe necessary to develop the ac- 
al friction loss and K’ is the correction to be applied to the length or to the 
iction loss computed on the basis of the total incoming flow in order to de- 


‘rmine the true friction loss: 
1.85 
Q 
K-( ) Cae wks pine pek pear CGO) 


Q 
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Expanding by the binomial theorem and substituting in the original equation: 


L 
a 1.85 0.85 x (1.85)(0.85) -0.15/ x \2 
b= K{ ax (Q, + 1.85 Q, bap eee Q (#4 
0 


2 (1.85)(0.85)(-0.18) g7h15 ( s 3) Faden ete. 4 


and 


h.=K LQ2*89 [4 - 0.9250(% + 0.2621/2-7 
: : Q Q 


+ 0,009820(%-\° + 0.0002260{ $-\* + 0.000809 (2-\? +...] ... 
Q. Q, Q 


The limits of the equation may be developed as follows: 


and 


2.85 
Y 1.85( 1 q 
hp=KLQ, Gato )t-0>8) st (' 


or 


seem 3)**) | 
«= aate (1- (1-8 ieee 


The values of Qe and K' computed from Eq. 9 and illustrated in Table 2 app’ 
only to the condition in which there are essentially an infinite number of cor 
nections. Table 2 illustrates the variations of these values as a function: 


TABLE 2.—VARIATIONS OF Q¢/Q; and K' FOR INFINITE NUMBER ‘ 
OF CONNECTIONS r 


ros [as [os [aa oe oe 
1,000 | 0.951 | 0,903 | 0.854 | 0.807 | 0,762 | 0.718 | 0.676 | 0.636 | 0.598 0.51 
1,000 | 0,910 | 0.826 | 0.747 | 0.673 | 0.605 | 0.541 | 0.486 | 0.434 | 0.389 0.3! 
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j/Q;, the ratio of the quantity withdrawn to the quantity entering the section. 
These show that the values of Qe/Q; range from 0.568 to 1.000; and the values 
of K' range from 0.351 to 1.000. 


FINITE NUMBER OF CONNECTIONS 


Values of friction loss for finite values of the number of connections n may 
9e computed by expansion of a series assuming takeoffs uniformly spaced L/n 


FIG, 3.—FINITE NUMBER OF CONNECTIONS 


listance apart and withdrawing a flow q/n at each takeoff. The procedure is 
Illustrated in Fig. 3 and the equations which follow: 
If the relationship kg = K L Q,1-85 = K L Q,1-85 kK" is expanded as in Fig. 


3, then 
1.85 
oe Wed Pe a AALS {2-20 
«| E(a, a (, -3) EY (ak 2) 


1.85 1.85 
(1-24) ee) mee eee se (Ait) 


i=n-1l : 
See [+-3! ae eer (42) 


Ay cgu tigi lege HM a Ate i 
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and 


se A ges Moh A ‘ 
h,=KLQ [3- 2 (a) ya+2+...4(n 1) 


2 
Pe seal a Tinie Oo et 2 tbe) ) in oe ee 
3 Q 

n t 

It should be noted that Eq. 12 reduces to Eq. 9 as n approaches infinity. The 
values of Qe/Qt and K' are found to lie within the same limits of 0.568 to 1.00( 
and 0.351 to 1.000 respectively. These values are given in Tables 3 and 4 an 
are plotted in Figs. 5 and 6. 


FIG, 4.—BEGINNING AND END WITHDRAWALS 


TABLE 3,—VARIATION OF Q,/Qt VS. NUMBER OF CONNECTIONS 


Ratio of Quantity Withdrawn (q) to Total Inlet Flow (Q) 


1,000 | 1.000 
0.976 | 0.951 
0.963 | 0.927 
0.955 | 0.912 
0.952 | 0.906 
0.951 | 0.903 


TABLE 4,—VARIATION OF K' VS. NUMBER OF CONNECTIONS 
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FIG. 5A.—SIMPLIFICATION (NO ERROR) 


FIG. 3.—ACTUAL CONDITIONS 
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FIG, 5.—NUMBER OF SERVICE CONNECTIONS VS. INLET FLOW CORRECTIONS 
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APPLICATIONS 


If the system is of the dead-end type, it isthe usual procedure to determine 
the flow withdrawn from each section of the main q and, by accumulation, the 
flow entering each section of the main Q,. The ratio q/Q; and the number of 
services, n, may be used to determine the equivalent length of the line (K' L), 
the true head loss, or the equivalent discharge through the main (Qe) for the 
correct friction loss. 

Example.—If 2.0 mgd enter the line and 0.6 mgd are withdrawn along the 
line from ten connections; the pipe is 1,000 ft long and C = 120: Qt = 2.0, q= 
0.6; n=10; q/Q; = 0.3; L=1,000. From Fig. 5 Q,./Q; is determined as 0.869 
and Qe = 0.869 x 2.0 = 1.738 mgd for correct friction computation. Also, hg = 
4.16 ft/1,000 ft or 4.16 ft. The same result may be obtained from Fig. 6 from 
which K' = 0.771 and Le = K' L or 0.771 x 1,000 = 771 ft for correct friction 
computation; hg = 5.40 ft/1,000 ft x 0.771 or 4.16 ft. 

Although this method is useful and precise with regard to dead-end sys- 
tems, it has very definite drawbacks when applied to loop systems. The draw- 
back, of course, is that in using the Hardy Cross method or a computor, the 
flows at the beginning and end of each line are not known in advance; and since 
q/Qt is not known, the quantity removed at the various junctions cannot be 
established. 


RATIO REMOVED AT THE BEGINNING OF THE LINE 


An approximation used in European practice3»4 and to a limited extent in 
the United States, is to assume that 45% of the withdrawal takes plate at the be— 
ginning of the line and 55% at its end. The various values of K' for this ap- 
proximation may be computed as follows: 

Let r q be the amount of flow withdrawn at the beginning of the section and 
(1-r) q be the amount withdrawn at the end of the section. The flow in the pipe 
would, therefore, be equal to Qt - rq and 


- 1.85 _ 1.85 _ 1.85 jy 
(Oe) aoprox LEQ. = LKQ, - ray" HF LKQ (KR), rox: + (15) 
(b), ores 7S (ee + 1.85 ae (-r q) 
1.85)(0.85) . -0. 
+ BONO) SOS Or ti] ce (16) 
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FIG, 6A.—SIMPLIFICATION (NO ERROR) 


FIG, 3.—ACTUAL CONDITIONS 
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Values of (K') approx were computed from Eq. 18 and the values compared 
with the actual values of K' computed by Eq. 12. The percentage errors were 
then computed and presented in Table 5. It will be noted that the error ranges 
from -5.7%to +3.5% for an infinite number of connections to -66.9% to 0.0% for 
a single connection. Accordingly, this procedure is correct only when there 
are a large number of connections since the values obtained for smaller num- 
bers of connections are all on the low side. 


TABLE 5,—PERCENTAGE ERRORS RESULTING FROM APPROXIMATION 
THAT 45% OF QUANTITY IS WITHDRAWN FROM BEGINNING OF MAIN 
AND 55% FROM ITS END 


The conservative simplification in which r is assumed as zero (that is, all 
of the quantity withdrawn from a section is assumed tobe withdrawn at its end) 
results in positive errors of from 0 to 185%. In contrast, when r is assumed 
as 0.45 (that is, 45% of the quantity withdrawn from a section is assumed to be 
withdrawn at the beginning of the section and 55% at its end) results in errors 
varying from +3.5%to -66.9%. : 


EXACT VALUES OF r AS A FUNCTION OF q/Q; AND n 


An exact value of r, the ratio assumed withdrawn at the beginning of the sees 
tion, can be readily computed from the expression Q, = Q; - rq, where Q, is 
the discharge which will produce the correct friction loss: 


. 

| 
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The ratio Qe/Q¢ is the ratio of the equivalent discharge to the total dischaw 
entering the pipe and the ratio q/Q; is the ratio of the quantity withdrawn in 
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ection to the total flow entering the section. These terms were developed 
reviously. 

From the values of Qe/Q, computed from Eq. 12 and tabulated in Table 3, 
1e eS values of r have been tabulated in Table 6 for various values 
fq /Qt- 

From the data in Table 6 it is readily apparent that r is more a function of 
than of q/Qt. The variation of q/Qt for any value of n varies as a ratio of 


TABLE 6.—VALUES OF r AS A FUNCTION OF n AND q/Q 


ooAN KE 


gr 


15 to 1 for the ranges q/Qt = 0.1 to q/Qt = 1.0. Obviously, a universal value 
rr cannot be used if a high degree of accuracy is required since r varies 
om 0.215 to 0.495 for the values shown in Table 6. 


APPROXIMATE VALUES OF r AS A FUNCTION OF n ALONE 


A study of the data illustrating the narrow range in r as a function of q/Q; 
iggested the possibility of approximating r as a function of n alone. Accord- 
gly, values of (K')approx were selected such that the negative error in the 


TABLE 7.—SELECTION OF TRIAL VALUES OF r 


umber of 0.97K' 
nnections H Minimum 


mputation of friction loss would not exceed 3%. Corresponding values of 
1.85 


/Qt were computed from the expression K' ( a) and trial values of r 
t 


re computed from Eq. 19for eachof the values of n. Minor adjustments were 
sequently made to spread the error more uniformly over the positive and 
ative range of errors. The trial values of r, which would be considered ap- 
ximately independent of q/Qt, are shown in Table 7. 
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Using the trial values of r, the corresponding values of (K')approx were 
computed from Eq. 18 and these values compared with the true values of K 
computed from Eq. 12. Errors were then determined for the complete rang 
of values of q/Q;t from 0.0 to 1.00 and the range of n from 1 to infinity. Th 
results are summarized in Table No. 8. 

The results of the study indicate that the approximation method using th 
suggested values of r will yield results ranging from -3.5% to +5.5%, well with 
in the range of accuracy of the pipe friction equations. Fig. 7 graphically il 
lustrates the relationship between the number of connections, n, and the ratio 
r, for the minimum range of error. 


TABLE 8,—SUGGESTED VALUES OF r 


Maximum Range of Error in K 
(+) 


Number of 
Connections 


Suggested 
Values ofr 


(0 tS es 
coohnNre 


A least-squares correlation indicates that the empiricalequation of best fi 
is an exponential : 


0.58 
r= 0.44 i a pc Obes tool PRR TRS (2¢ 


Analysis of the deviations using the Student “t” test indicates a significanc 
of greater than 99.9%. However the equation is too complex for routine us¢ 
and graphical selections of values of r appear more useful. 


APPLICATIONS 


Since this procedure is independent of the ratio q/Qt, it is therefore paz 
ticularly suitable for use in network analysis, since the ratio need not be know 
in advance of the final network computations. However, adjustments may k 
necessary if the wrong direction of flow was assumed for the line. The prc 
cedure is simple to apply. 

Example.—If 0.6 mgd are withdrawn from ten connections along 1,000 ft 
12 in. pipe, C = 120, n = 10. From Fig. 7, r is determined as 0.40, and tl 
amount assumed to be withdrawn at the beginning of the line is 0.40 x 0.6 mgs 
or 0.24 mgd, and at the end of the line, (1 - 0.40) x 0.6 mgd or 0.36 med. 
calculations then proceed in the usual manner. If, after the calculations, 1 
incoming flow is determined as 2.0 mgd, the friction loss in the pipe is dete: 
mined on the basis of Q;-rq or 2.0 mgd- 0.24 mgd or 1.76 mgd. The fricti 
po therefore 4.25 ft, which compares within 2.2% of the correct value 
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RADIOACTIVE TRACERS IN HYDROMETEOROLOGY®)? 


By L. Machtal 


SYNOPSIS 


Atmospheric motion is studied with the aid of radioactive tracers. Radio- 
‘tive materials are used to improve the meteorologist’s ability to track the 
th of masses of water in the atmosphere. The most promising tracer seems 

be tritium but the large number of sources of tritium complicates its appli- 
ition to hydrological problems. 


INTRODUCTION 


It is most reasonable that the study of the movement of atmospheric water 
por be conducted with water vapor itself as the tracer. Almost all hydro- 
eteorological research follows such procedures. But, to the conventional 
nsing instrument one water vapor molecule is no different from another. 
racking a specific mass of water in the atmosphere simply by successive 
servations of atmospheric humidity is often quite difficult. Rather, the path 
a mass of water vapor is normally computed from the wind observations. 
it today other possibilities are available since water vapor molecules can be 
entified by special labels; the isotopes of oxygen and hydrogen. In addition, 
. air mass containing the moisture can be identified and tracked by radio- 
tive aerosols or gases. 


Note.—Discussion open until September 1, 1960, To extend the closing date one 
onth, a written request must be filed with the Executive Secretary, ASCE. This paper 
part of the copyrighted Journal of the Hydraulics Division, Proceedings of the Amer- 
in Society of Civil Engineers, Vol. 86, No. HY 4, April, 1960. 

4 Research supported by the Div. of Biology and Medicine of the U.S. Atomic En- 
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In this paper, we shall review some of the information dealing with tritium 
and other radioactive tracers in the atmosphere. 


RADIOACTIVE TRACERS FOR AIR MASSES 


Water vapor in the atmosphere spreads with the air so long as it remains 
uncondensed. Horizontal transport and mixing by the atmosphere are vital ele- 
ments in the understanding of the atmospheric link of the hydrologic cycle. 

Air masses have been tagged by fission products from weapon tests. The 
measurements of this radioactivity, as nuclear clouds move across the United 
States from the Nevada Test Site, provide both the trajectories and diffusion 
rates. We have been able to verify trajectories2 for distances of 1,000 miles 
or more for thefirst time at altitudes below 30,000 ft. It turns out that because 
of better wind observations at lower altitudes our ability to reconstruct the air 
paths below 30,000 ft is better than might have been guessed from the verifi- 
cation of balloon flights above 30,000 ft. 

But more important is the fact that one now has more knowledge of the rate 
of lateral diffusion for the atmosphere.2,3 Previously, there were two other 
sources of information on horizontal mixing for air trajectories to great dis- 
tances: the first was derived fromthe spreadof moisture itself4 and the second 
based on the spread of pairs of constant level balloons.° The three sets of ob- 
servations from moisture, from balloons, and from radioactive clouds make 
an interesting pattern. The separation of balloon pairs suggests a coefficient 
of horizontal mixing of about 107 cm2 sec-1;the spread of radioactive clouds, 
about 108 cm2 sec-1; and the spread of moisture, about 109 cm2 sec-1. An 
explanation of the reason for the difference in these diffusion rates may help 
us understand the nature of mixing in the atmosphere. 

Itis the writer’s view that much of nature’s horizontal mixing of water vapor 
or other gases depends on wind shear and vertical diffusion. This phenomenon 
is schematically illustrated in Fig. 1. Fig. 1(a) is a top view, looking down on 
a cylinder which, at t = 0, is vertical (for the sake of argument). The atmos- 
phere’s winds vary in direction and speed with altitude so that at a later time, 
t = tj, the cylinder is stretched out. The winds in the upper part of the cylin- 
der were from the southwest and near the bottom, lighter and from the west. 
Looking down at the top at the later time, we find that each horizontal section 
of the initial cylinder has grown. The initial size is the inner circle. The 
cylinder would still have this radius if there were no growth. The closely 
stippled area, also circular in shape, reflects the growth of the cloud due tc 
true diffusion. In horizontally isotropic turbulence, the spread is the same in 
all directions. This growth might be derived from the coefficient of horizontal] 
mixing of the order of 107 cm2 sec-! as derived from balloons. The outer 
cross-hatched area shows the growth due to shear and vertical turbulence. As 
the cylinder tilts, as seen in Fig. 1(b), it is evident that vertical mixing will 
bring tracer material upward from below and downward from above. This pro- 
duces a further growth in the horizontal dimensions of the cylinder. The ad- 
ditional enlargement takes place in the line of the wind-shear vector as seel 


2 Machta, L., Hamilton, H. L., Jr., Hubert, L.F., List, R. J, and Nagler, K. M. 1957 
J. Meteor., 14, 165. ; 


3 Wilkins, E., 1958, Trans. Am. Geophys. Union, 39, 58. " 
4 Miller, J. E., 1948, Met. Pap. New York Univ. 1. 
Moore, C. B., Smith, J. R. and Gallswyk, A., 1954, J. Meteor., ped LG a 
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n Fig. 1(a). Thus, at t = ty, the cylinder has grown to the dimensions of the 
uter heavy line. From ordinary measurements of concentration one would be 
mable to discriminate which part of the growth is due to horizontal diffusion 
Jone and which part is due to the effects of shear and vertical diffusion. 

The measurements of the spread of radioactive debris suggests that this ex- 
lanation is correct. We find that the north-south extent of the cloud at say, 
0,000 ft is greater when the trajectories just below and just above 20,000 ft 


t=O t=t, 


‘Ss ORIGINAL SIZE 


GROWTH DUE TO NORMAL 
HORIZONTAL MIXING 


GROWTH DUE TO 
WIND SHEAR AND 
(a) TOP VIEW VERTICAL MIXING 


(b) SIDE VIEW 


FIG, 1.—VERTICAL COLUMN OF AIR 


2 to thenorth and south of the 20,000-ft path as they cross the erty samp- 
line. The coefficient of horizontal diffusion, 108 cm2 sec~!, noted pre- 
ously, is derived from those cases of very little directional shear. It might 
noted that these cases were in the minority. Further, all cases have wind- 
eed shear; that is, the wind speed changes with altitude. The growth due to 
ear and vertical diffusion, in the case of speed shear, is along the trajectory 
ther than laterally and is not detected by the kind of aircraft-sampling oper - 
ic h there is information. 

Be ewe diffusion coefficient for moisture appears to be one to 
ders of magnitude greater than for balloons or radioactive debris. The coef- 


ra 


s 
| a 


52 April, 1960 HY 


ficient obtained from the lateral spreadof radioactive clouds is higher than tk 
spread of pairs of balloons partly because the effect of shear and vertical dif 
fusion is not entirely eliminated. This indicates that the shear and vertic: 
mixing phenomena are of equal or greater importance than ordinary horizont: 
diffusion. Thus, water vapor clouds will spread horizontally fastest in regior 
where there are the greatest vertical wind shears. 

On a larger scale, we have followed particulate debris from the Pacifi 
atomic tests over almost the entire world. For example,§ in Fig. 2(b), we ar 
shown isolines of deposited radioactivity in millicuries per 100 sq miles : 
ground level during the 30-day period following the March 1, 1954, hydroge 
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FIG, 2.—HISTORY OF BRAVO CLOUD 
bomb test at Bikini. The hatched area of Fig. 2(b) is the intertropical cor 
vergence zone—the meteorological separation between northern and southe1 
hemisphere air. Fig. 2(a) shows the paths of various parts of the nuclear clot 
for only the first few days. The arrows toward the west show the sense of ¢ 
bris transport in the trade winds, below about 20,000 ft, while the part of 

cloud which movedinitially toward North America lay in the 20,000 to 50,00 
ft altitude range. All of these paths reflect the east-west or west-east moti 
which the weatherman believes to be characteristic of the atmosphere’s mo 
ment. The numbers in Fig. 2(a) indicate the number of days between detonat: 


6 Machta, L., List, R. J. and Hubert, L. F., 1956, Science, 124, 474, 
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ind first ground observation of fission products. The isolines of deposited 
adioactivity in Fig. 2(b) confirm the more rapid east-west rather than north- 
jouth transport of tracer material. Pieces of the cloud are torn from its nor- 
nal west to east orientation, as in the United States where southerly winds 
head of a low pressurearea carried debris into the midwest in about 2 weeks. 
fhe slower, general northward spread did not reach the west coastof the United 
tates until weeks later. Only small amounts of radioactivity reached into the 
jouthern hemisphere. But equally small amounts are reported in polar lati- 
udes. The clouds from the Hardtack Pacific test series in 1958, the only 
jummer Pacific tests, apparently were carried into the southern hemisphere 
n greater amounts than spring or fall test clouds. 

The tracking of radioactive debris has, by and large, confirmed our notions 
bout the broad movements of air in the troposphere, and some of the details 
f transport in both time and space. 


RADON—A NATURAL RADIONUCLIDE 


Radon is a naturally occurring radioactive gas with a 3.7-day half-life which 
manates from the rocks and soil of the land. Thus, its source, broadly speak- 
ng, is the complement of water vapor in the atmosphere. Where one finds 
nore water vapor one should also find less radon and vice versa. Insofar as 
ceanic areas have high moisture and low radon while landareas have the re- 
erse, this is true. But the atmosphere has a limited capacity for moisture, 
yhereas this is not true of radon. 

The similarity between radon and water vapor, both evolving from the earth’s 
urface, permits one to deduce information about the vertical fluxof water vapor 
rom observations of radon. Fig. 3 is a plot of radon concentration’ as the 
rdinate in Fig. 3(a) and thermal stability in the lower atmosphere as the ordi- 
ate in Fig. 3(b). The more positive the number in Fig. 3(b), the more stable 
ne atmosphere. The abscissa is local time of the day. The observations were 
1ade at Harwell, England, and are averaged over a period of several months. 
‘he zero gradient line in Fig. 3(b) is the adiabatic lapse rate, negative values 
re superadiabatic lapse rates, and positive values, more stable than the adia- 
atic lapse rate. The 0.24° C/76 ft gradient value is the isothermal lapse rate 
f actual temperature and values greater than 0.24 correspond to temperature 
iversion conditions. 

One notes a diurnal cycle in both curves. Throughout the night, the rate of 
mission exceeds the rate of upward mixing since, as seen in Fig. 3(b), the air 
; stable and resists vertical turbulence. With the beginning of solar heating, 
1e air becomes unstable and the upward flux of radon exceeds the rate of 
missionfrom the ground. The radon concentration decreases. The assumption 
; made that the rate of emanation from the ground is approximately constant 
nd that the concentration reflects the change in vertical diffusion. This cycle 
f radon concentration at ground level is no surprise since the diurnal vari- 
tion of vertical mixing is well known. 4 

Harry Moses, M. ASCE, and his colleagues at Argonne National Laboratory 
2ar Chicago, Ill., have made measurements at the ground (within 1/8 in. of 
\e ground) and up to 131 ft on a meteorological tower. The observed diurnal 
yelefor a summer day with light winds and cléar skies is given in Fig.4. The 


7 Gale, H. J., Peaple, H. J., 1958, Intern, J. Air Pollution, 1, 103. 5 
Moses, H., Stehney, A. F. and Lucas, H. F., Jr., 1959, Radiological Physics Divi- 
on Rep., ANL 5967, Argonhe National Laboratory, 165. 
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uppermost curve refers to the 1/8-in. level and the lowermost curve the trace 
at 131 ft. Something akin to the normal diurnal cycle is observed up to 18.75 ft 
but not at the 131-ft level. The amplitude of the diurnal cycle diminishes with 
height and the peak concentration is delayed in time during the night. The wine 
speeds during the night were very light and there was extreme thermal stabili- 
ty. The almost complete lack of atmospheric turbulence prevents the radon 
from being transported up to the 131-ft level until 5:00 A.M. At about 5:00 A.M., 
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FIG, 3.—DIURNAL CYCLE OF RADON CONCENTRATION AND THERMAL STABILIT" 


- the sun heats the ground and abruptly increases vertical mixing so that, at 131 ft 

the radon suddenly rises. Then, with deeper layers of the atmosphere involve 
in the mixing, the upward flux at 131 ft due to turbulence exceeds the growt 
due to ground emission and the concentrations at all observed levels on th 
tower diminish after 7:00 A.M. If we had observations at greater altitudes, w 
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could obtain the rate of growth of the ground layer “blanket” of air during the 
daytime also. It is clearly greater than 131 ft, of course. 

The vertical profile of radon or water vapor emitted from the ground or 
water can be converted to a measure of vertical turbulence, say, the coefficient 
of vertical diffusion, if the flux is known. It is here that radon may have an 
advantage over water vapor, for the radon emission rate is largely independent 
oi atmospheric conditions. The emanation rate can be estimated from these 
observations and has been given in the literature. The diffusion coefficient 
derived from the lowermost points at 1/8 in. and 3.17 ft in the middle of the 
night turns out to be about one order of magnitude larger than molecular dif- 
fusion, 1 cem2 sec-l, one of the lowest values of eddy diffusion in the litera- 
ture. In the afternoon, the coefficient inthis same layer is computed to be about 
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FIG. 4.—RADON CONCENTRATION FOR JULY 16-17, 1958. 


00 cm? sec-1. The night values in the 3.17 ft to 131 ft layer are also about . 
00 cm2 sec-! while daytime values are 1,000 (or more) cm? sec~!., 

In Fig. 5, similar results for the summer season are given for a cloudy 
eriod with less incoming or outgoing radiation. Wind speeds were somewhat 
igher than in the previous experiment and the night time thermal stability was 
uch less marked. For example, in the previous July case the temperature at 
5 ft on the tower was as much as 2° Celsius warmer than at 5.5 ft, but in the 
resent case, the temperature was never more than 0.5° Celsius warmer. The 
/8-in. concentration in the afternoon is about the same as in the last case but 
- night the maximum is less than half that of the July night with its strong 
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thermal inversion. The coefficient of vertical diffusion during this night is 
about five times greater than the previous case. The vertical mixing in the 
3.17 ft to 131 ft layer was also greater than before and this is reflected in the 
small diurnal amplitude in the radon concentration at all levels. It is interest- 
ing to note that the 1/8-in. concentration decreases fairly steadily from befor: 
midnight while the upper levels show an increase during the same interval 
These data confirm the differences in measurements of atmospheric proper: 
ties right above the interface from which the emission occurs as compared witl 
conventional anemometer -height observations. 

One must, of course, be cautious in extrapolating any conclusions found fron 
radon measurements to those of evaporating water surfaces: the ground has ; 
different roughness than a water surface, the ground is a better radiator an 
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FIG, 5.—RADON CONCENTRATION FOR AUGUST 27-28, 1958. 


absorber of heat than the water and has a larger diurnal cycle of stability, an 
the radon emanation rate may be too variable to attribute all the changes it 
concentration to atmospheric mixing. The emission from soil will not var’ 
greatly if not too windy days without snow or water cover only are studied 
Despite these and other cautions, it is likely that micrometeorology has mucl 
to learn from a study of the histories of radon concentration profiles. 


ee Ee SL ih 


NATURAL TRITIUM 


Tritium is a radioactive isotope of hydrogen with a halflife of 12.5 vntordil 
by the action of cosmic rays. It is produced in greatest quantities at abou 
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50,000 ft in the atmosphere. Thus, the main source of tritiated water vapor is 
the upper troposphere or stratosphere rather than the surface, as is the case 
for ordinary water. 

Tritium is normally expressed as T units, that is the number of tritium 
atoms per 1018 hydrogen atoms inthe sample. The limit of detectability of 
tritium is 1 T unit or slightly less. Very recently-disclosed observations of 
tritium in the stratosphere? give tritium concentration as about 105 T units, 
give or take an order of magnitude. Early 1953 rains at Chicago contained 
about 1 to 10 Tunits.19 Oceanic rains at the same time had about 0.5 to 3 T 
units .10 

The high stratospheric concentrations reflect the greater production of tri- 
tum in the stratosphere. But an additional reason for high concentrations is 
the very small amounts of moisture in the stratosphere. The stratospheric 
relative humidity is between about 1% and 15% and the temperatures between 
-50° and -80° Celsius. If the moisture content were higher, the number of 
Tritium atoms would not change, of course, but the number of hydrogen atoms 
would increase and the tritium-to-hydrogen ratio would decrease. When the 
1igh tritium concentrations are carried into the troposphere, they are diluted 
with the large amounts of tropospheric moisture and quickly reduce their con- 
-entrations to tens or hundreds rather than thousands of T units as in the 
stratosphere. 

Fig. 6 shows the time history of the tritium content of rains at three North 
American stations: (a) Chicago, (b) New York, N. Y., and (c) Ottawa, Cana- 
ja.11,12,13 he abscissa in Fig. 6 is time, increasing to the right, and the 
yrdinate, on a logarithmic scale, is the precipitation tritium concentration in 
[T units. Each collection point is shown with a separate symbol. Although the 
-anadian laboratory has intercalibrated with the Chicago laboratory toits own 
satisfaction, these and other data suggest that there is still a systematic dif- 
erence between the United States and Canadian laboratories with the Canadians 
rreater by a factor of about 2. Fig. 6 shows that the rains after March, 1954 
ire contaminated by man-made tritium. 

It is of some interest to explain the observations of natural cosmic-ray tri- 
ium found in rains. For simplicity, consider two reservoirs of water vapor: 
irst, the moisture in the upper troposphere with many tens or hundreds of T 
nits of tritium concentration. This moisture originated from the oceans so 
ong ago that it now canbe called “non-recent oceanic tritium.” The second is 
yater in the lower atmosphere which has recently evaporated from the oceans. 
"rom 1953 measurements of oceanic concentrations of tritium, we know that 
hese water bodies contain 1 T unit or less. All water vapor evaporating from 
he (oceans possesses about the same tritium concentration. Rains at coastal 
tations or island locations are made up of water which has almost entirely 
vaporated from the oceans recently and reflects tritium concentrations not 
uch greater than oceanic water. Thus, Hawaiian rains had an average tritium 
oncentration of lessthan 1 T unit in 1953. But continental rains are mixtures 
f the two sources of moisture. The bulk of the water is from the relatively 
ecent oceanic source but the greater part of thetritium comes from the upper 
“oposphere moisture reservoir. The day-to-day differences in tritium con- 
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9 Hagemann, F., Gray, J., Machta, L. and Turkevich, A., 1959, Science, 130, 542. 
10 Von Buttlar, H. and Libby, W. F., 1955, J. Inorg. Nucl. Chem., 1, 75. 

11 Brown, R. M. and Grummitt, W. E., 1956, Can. J. Chem., 34, 220, 
12 Giletti, B. J., Bazan, F. and Kulp, J. L., 1958, Trans, Am, Geophys. Union, 39, 807, 

13 Begemann, F. and Libby, W. F., 1957, Geochem. et. Cosmich. Acta., 12, 277. 
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centrations in rains over continents can be attributed in large part to differ - 
ences in the fractions of water from the two sources. This explanation, first 
offered by W. F. Libby, accounts for the greater amount of tritium in conti- 
nental than in oceanic rains even before there were bomb sources of tritium 
of an appreciable magnitude. 

The tritium concentration of oceanic water vapor is known. If one can esti- 
mate the tritium concentration of the non-recent oceanic reservoir, it becomes 
possible to compute the fraction of moisture derived from each source. Asa 
first and rather naive approximation, one can assume the tritium -to-hydroger 
ratio in the middle and upper troposphere, the non-recent oceanic source, to be 
inversely proportional to the amount of moisture at this level relative to the 
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stratosphere or about 100 T units. Simple arithmetic gives the following re 
sults: for a rain with 5 T units, about 5% of the water is from this non-recen 
oceanic source, for 10 T units about 10%, and for 40 T units, the non-recen 
oceanic source constitutes over 60% of the liquid water. Inspection of Fig. | 
before March 1, 1954, indicates the frequency of occurrence of the variou 
fractions of oceanic and non-recent oceanic moisture. The first half of 195 
may be a more reliable period than the second half of 1953 because of the lesse 
likelihood of artificial tritium. During this period, 5% to 10% of the rainwate 
came from non-recent oceanic water. Among the possible shortcomings ¢ 
this calculation is the assumption that the tritium concentration of a fallin 
raindrop is constant as it passes through water vapor with a different tritiur 
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mcentration. B. Bolin has recently calculated the exchange rate and found it 
) be surprisingly fast.14 

Some day, many tens of years from now, one may once again be able to use 
utural tritium when the bomb tritium decays away and becomes thoroughly 
lutedin the oceans. Today, one must deal with tritiumfrom man-made sourc- 
3 as well. Actually, for many geophysical purposes, the injection of large 
nounts of tritium into the atmosphere permits studies which would otherwise 
t be possible. Some of these will be discussed herewith. 


THE USE OF BOMB-TRITIUM IN HYDROMETEOROLOGICAL PROBLEMS 


The data in Fig. 6 show a systematic decrease in concentration after the 
ing of 1954 which can be approximated by a straight line with a slope corre- 
nding to a half-time of about 40 days. That is, the tritium concentration de- 
eases to half its value each 40 days. Since the radioactive decay of tritium 

12.5 yr, the decrease is not due to decay. Rather, we are certain it reflects 
eloss of tritium-rich water vapor from the atmosphere through precipitation. 
1e water vapor which evaporates from the oceans to replenishthe water con- 
ins less than 1 T unit, or a negligible concentration of tritium. The transfer 
the bomb tritium is thus one-way, from atmosphere to ground or oceans via 
ecipitation. 

The 40-day half residence time deduced from the graph must betreated with 
ution. It does not agree with much shorter estimates made by more con- 
ntional methods suchas the “7” day half residence time given by H. Lettau.15 
| explanation of this difference is necessary. First of all, the tritium con- 
ntration should be growing in North America at this time as the tritium dif- 
ses northward from the Marshall Islands. If this diffusive growth of tritium 
ncentration is subtracted fromthe decrease in Fig. 6, the half residence time 
comes less than 40 days. This fact cannot explain the discrepancy. 

There is a difference in the residence time for atmospheric aerosols de- 
nding on their altitudeor origin. Thus, dust which originates at ground level, 
d which is tagged, for example, with radon daughter products which are 
dioactive, suggest a half-residence time of about 1 day! while C. Junge and 
Gustafson17 estimate about a 3-day half-residence due tothe rainout of sea 
it particles. On the other hand, the dust which is injected into the upper 
posphere (say 20,000 ft to 40,000 ft) by nuclear explosions has been observed 


9€ removed with a half time of from about 15 to 30 days. 18 P.S.Goel et al.19 
ind a 30-day half residence in the tropics for dust in the upper troposphere 
ged with cosmis ray radionuclides. 
Removal by precipitation is probably important for both low and high alti- 
le dust. But precipitation scavenging is usually limited to levels below 
000 ft to 15,000 ft, on the average. Dust or water vapor originating from 
ground can be removed more quickly than dust or water vapor originating 
ll above the rain bearing layers which must mix downward before being 


14 Bolin, B., 1958, Proceed, Sec. Intern. Conf. on the Peaceful Uses of Atomic Energy, 
| a’ 1954, Arch. Meteor. Geophys. Biokl., A, i, U33. 

L6 T.ehmann, L. and Sittkus, A., 1959, Naturewissenschaft, 1, 9. 

8 Bee NG. Oomoad: BG; De resist z ni Bisby, E. M. R., 1957, A.E.R.E. 
9 Beal: A aaa ah N., Probhakara, C., Thor, R. and Zutshi, P. K., 1959, 


lus, 11, 91. 
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washed out of the air. This probably accounts for the discrepancy between the 
longer residence time derived from the bomb tritium and other estimates. This 
tritiated moisture was added mainly to the upper troposphere while ordinary 
water vapor is added from the earth’s surface. The same longer tropospheric 
residence time would apply to natural tritium feeding down from the strato- 
sphere. 

The initial impulse of tritium in the rains at New York and Ottawa exceedec 
1,000 T units. This means that water vapor from the Marshall Island areé 
contributed measurable moisture to rains in eastern North America. The same 
kind of arithmetic used to obtain thefraction of non-recent oceanic rain can be 
applied to the post-March 1954 period. To make such a computation we neet 
the estimate of the tritium concentration of the air contaminated by the boml 
tritium. The estimate of the amount of tritium made and available after th 
atomic cloud stabilizes in the troposphere and the degree of dilution by the time 
it reaches eastern North America, is the main uncertainty in the calculation. 

The fraction of the rain water containing tritium from the Marshall Island: 
is roughly 5% using such estimates as the writer could make for the bomb: 
tritium concentration. This number is much higher than would have been ex: 
pected. Is it possible that moisture from this limited equatorial region con 
tributes so much water to the rain in eastern North America? 

The arrival of bomb tritium in eastern North America is also faster thai 
might have been expected but here there is less uncertainty. The first rail 
with anomolous tritium occurred 10 days after the explosion although if it raine 
a day or two earlier, we might have had evidence of even earlier arrival. Fas 
though this is, we have evidence from the rainout of fission products durin; 
the Greenhouse test series of Marshall Island air passing over the St. Lawrence 
River region in 5 days. 


CONC LUSIONS 


This paper has reported on only a few examples of the use of tracers in th 
study of atmospheric motions but it is felt that they illustrate the usefulnes 
of this new tool in hydrometeorology. Water can be tracked by both stable an 
radioisotopes of oxygen and hydrogen and, in fact, a considerable effort ha 
been and is being devoted to this activity. Of these tracers, tritium seems t 
be one of the more exciting. But what are the prospects for current and futur 
research using tritium? . 

The atmosphere’s tritium is not in equilibrium with the oceans as was th 
case before man’s tritium came on the scene. The large number of tritiuz 
sources complicates its application to hydrological problems. As long as th 
atmospheric high yield nuclear tests are suspended, we can have nonew surge 
of tritium for study. It thus appears that we have either of two prospects fo 
the future; first, make anintensive study of the distribution of tritium in natur 
so as to be able to constructively use the variability, or, second, add new tri 
tium for small-scale experiments. The writer emphasizes small scale be 
cause to “see” new injections of tritium in the atmosphere over large areas ¢ 
the earth requires the addition of more tritium than we probably are willing 1 
use. At this time (1960), the additionof tracer tritium seems to hold the greate 
promise. ; 
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SEDIMENT PROBLEMS OF THE LOWER COLORADO RIVER 


By Whitney M. Borland,! M. ASCE and Carl R. Miller,2 M. ASCE 


SYNOPSIS 


The controlof a river to serve the purposes of manalters the river regime 
thereby creates a number of varied problems. The closure of Hoover Dam 
he Lower Colorado River, and the subsequent construction of other major 
ictures downstream, instigated a series of river adjustments that have re- 
‘ed corrective or protective measures. Some of the major problems that 
e been encountered are pointed out, methods of rectification and design 
siderations are presented, and some of the results that have been obtained 
given. Some probable future river problems are also briefly examined. 


INTRODUCTION 


‘here are many interests to be considered in any river control work in the 
er Colorado River Basin. State and Federal fish and wildlife agencies, 
idary commissions, recreation advocates, adjoining property owners, and 
rshave a vital interestin the river and its use. The material in this paper 
be restricted to the technical considerations involved in developing river 
rol operations by the Bureau of Reclamation (USBR), and no attempt will 
1ade to discuss the considerations of the other interested agencies. In the 


2te.—Discussion open until September 1, 1960. To extend the closing date one 
h, a written request must be filed with the Executive Secretary, ASCE. This paper 
rt of the copyrighted Journal of the Hydraulics Division, Proceedings of the Amer- 
Society of Civil Engineers, Vol. 86, No. HY 4, April, 1960. 
Head, Sedimentation Sect., Hydrology Branch, Div. of Proj. Investigations, Bur. of 
ee hed Tech, Research Branch, Soil and Water Con 
upervisory Hydr. Engr., Watershed Tech. Research Branch, n- 
de ee reek Research Service, Beltsville, Md. (Formerly with the 
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planning of river control work, however, all factors and interests involved ar 
duly considered, and the plan, when developed, is one that includes the variou 
river uses and the best engineering planfor present and future river operation 

Engineers associated with the planning of the Boulder Canyon Project, o 
which Hoover Damis a major feature, realized that the closure of Hoover Dar 
would cause extensive changes of river regime in the more than 300 miles c 
river channel, downstream, which extends to the Gulf of California. However 
the details of these changes and the introduction of new problems could not hav 
been fully anticipated. Fig. 1 shows the river course and structural contre 
between Hoover Dam and the Gulf. 

Estimates based on available information and data collected prior to th 
closure of Hoover Dam indicated that the average sediment load being trans 
ported past the damsite was in excess of 160,000,000 tons or 103,000 acre-1 
annually. A resurvey of Lake Mead in 1948, and sediment measurements afte 
1948, at the Grand Canyon gaging station, have confirmed this estimate. Trap 
ping this amount of ‘material in Lake Mead and releasing reregulated clea 
water to the channel downstream was certain to mark the initiation of river 
channel adjustments. These adjustments, increased by the later completion ¢ 
the Davis (1950), Parker (1938), and Imperial (1943) Dams, have extended fro1 
Hoover Dam to the Gulf of California, as predicted (Fig. 1). Public Law Ne 
469, enacted by Congress in 1946, established the framework under which th 
study, planning, and control of the Lower Colorado could be accomplished. 

History.— Prior to construction of Hoover Dam, the river moveduncontrolle 
through deep canyons and wide, alluvial valleys. Rampaging floods occurré 
in the springtime as the mountain snowpack melted, but the river became 
relative trickle in the late fall months. The sands and silts in the canyon be 
and the alluvial fill in the flood plains were alternately shifted, washed ou 
and replenished as the millions of tons of material moved down to the Gulf. 
is significant to note that the bulk of the water and sediment carried by tl 
Colorado River is derived from the drainage area above Hoover Dam. Bela 
Hoover, the river flows through desert terrain where the only water and sed 
ment reaching the river is derived from erratic, high intensity, summer storm 
In the valley flood plains, the river meandered and changed course frequentl 
destroying a parcelof land in one place and creating new land in another. Du: 
ing low run-off years, the valley areas suffered from lack of water, and duri 
high run-off years, the suffering resulted from floods. Historical accounts 
the fickle nature of the untamed Colorado River are numerous.’ Fig. 4 illu: 
trates the control of the Lower Colorado River effected by Hoover Dam. 

Lake Mead Sedimentation.—As mentioned previously, the sediment load 
the Colorado River, as it enters Lake Mead, averages about 103,000 acre- 
annually. This extremely large sediment inflow has led many to believe th 
Lake Mead would lose its usefulness in a relatively short period of time. Th 
is certainly far from the truth. Although the sediment inflow is large, t 
reservoir capacity is also large and, therefore, the annual per cent capaci 
loss is reasonable. At the present rate of sediment inflow, the 31,250,000 acr 
ft original capacity of Lake Mead would be reduced to a negligible amount 
about 300 yr, not considering additional life which may be added due to cot 
_ paction of deposits and deposition above spillway crest elevation. Howeve 
with the closure of GlenCanyon Dam, now under construction upstream (Fig. | 
studies indicate the future sediment inflowto Lake Mead will be only one-fou1 


oe ee ee eee 
3 “Silt in the Colorado River and Its Relation to Irrigation,” by Samuel Fortier ¢ 
Harry F. Blaney, Tech. Bulletin 67, U. S. Dept. of Agric., February, 1928. 


iy 4 


% 


caine 


SEDIMENT PROBLEMS 


UTAH z 
Y St.George Kanab \ © 
re i Te ON 
M oes be A S Lees Ferry 9 
4 ) $f MARBLE 
s) CANYON GLEN 
<{DAM SITES¢ CANYON 
DAM 
Las V eX& 
ooe awe Z a 
Boulder, S z 
city ° \ HOOVER Grand Ay 
he DAM VR Canyon 
orado 
ela BRIDGE CANYON coc 
SS oe DAM SITE Rotnene 
SS ‘. f 
ehiees Bulheo a 
RQ | Z OFS Naas A 


Needles 3 
\ 


CALIF. oe 
See < o oF hregeie 


pau oiiarns Alamo Z 

SA HEADGATE 

¢ R 

a BAM 

~ ¥ PALO VERDE ope esaany 
Blythe of DIV. DAM SCALE OF MILES 


Palo Verde od 


! 


ay IMPERIAL DAM 
LAGUNA 
DAM 


N. 
ieee 


MORELOS ¢ f “Yuma 
DAM San Luis 


navatoe 


“FIG. 1.-LOWER COLORADO RIVER, LEES FERRY TO GULF OF CALIFORNIA 


64 April, 1960 HY 4 


of the present rate, and the life of Lake Mead will be extended beyond 1,000 yr. 
The life of Lake Mead can be expected to be prolonged even further with con- 
struction of additional dams and reservoirs between Glen Canyon Dam and 
Lake Mead, on the main stem of the Colorado River, andon the sediment- 
carrying tributaries such as the Little Colorado and Paria Rivers (Fig. 1). 


AVAILABLE DATA 


Data available, in the Lower Colorado River basin, that can be utilized ir 
the determination of existing conditions, and in planning of river rectification 
work include: 


TABLE 1.—SUMMARY OF RESULTS OF CROSS-SECTIONING 
LOWER COLORADO RIVER? 


River 
Miles 


River 
Reach 


Quantities 
(1000 cu yds) 


Period 


Hoover Dam 1935 0. =- 1222¢ - 10,402 
to to 12.2- 25.6 - 13,671 
Lake Havasu 1951 25.6- 42.7 - 25,716 
42.7- 64.2 - 48,125 
64.2- 82.9 - 39,515 
82.9- 91.8 - 14,401 
91.8- 98.0 + 8,281 
98.0-111.8 + 94,737 

111.8-119.4 - 


Davis Dam 
to 
Lake Havasu 


- $3,541 


Parker Dam - 190,100£ 
to + 148,970 
Imperial Dam 


4 Report of River Control and Investigations, Lower Colorado River, Nov wgmbes 
1957.. b Aggradation = +; degradation =-, © River miles below Hoover Dam. 

miles below Davis Dam. "e River miles below Parker Dam. f Determined from ia 
cross sections and total sediment transport computations. 


a. River cross sections.—These cross sections, extending from Hooyl 
Dam to the International Boundary, are spaced approximately 1.2 miles apar 
and have been resurveyed at regular intervals since Hoover Dam closure. Re 
sults from the resurveys are summarized in Table 1. Cross sections are no 

being surveyed by the use of echo-sounding equipment where feasible. z 
b. Suspended sediment sampling.—Suspended sampling was initiated : 
Yuma in 1911, and at Imperial Dam and Red Cloud (near Taylor’s Ferry) i 
1933 (Fig. 3). A detailed, periodic, total load sampling program was initiate 
in 1955, between Davis and Imperial Dams to aid in the design of river channel 
ization work. The total load sampling stations are shown on Figs. 2 and 
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Total sediment load is presently being determined by the Modified Einsteii 
Procedure.4 Sampling under the current program is being accomplished by 
the equal-transit-rate method (ETR). This type of sampling gives a bette 
evaluation of the sediment movement taking place. A sample is collected a 
each of 25 verticals, spaced across the section. The 25 samples are combinet 
or composited into one sample for analysis. The results of these samples an 


Palo Verde 
Canal Intake — \. Palo Verde 
: Div. Dam 
Ehrenberg 


Taylors Ferry 
Palo Verde Drain— 


* * 
Cibola Valley 
Adobe Ruin 


Approx. Location 
Red Cloud Cable 


~ Arroyo Sen 
° 


Picacho 


Pilot Knob Wasteway 
and Power ay, 


JIMPERIAL DAM 


) —LAGUNA DAM 
Laguna 


Riy 
Yuma 


\MORELOS DAM 


: 4 Gaging Station 
joan Luis @ Total Load Sampling Station 


Canal 


10 (0) 10 20 


y 
Mey wei : : 
To SCALE OF MILES 


FIG. 3.—-LOWER COLORADO RIVER, PALO VERDE DAM TO 
INTERNATIONAL BOUNDARY 


computations are shown on Table 2. Average-size analysis of suspended an 
bed materials for the total load stations are shown in Fig. 5. Information ob 
tained includes the sediment loadat Needles, Calif., the inflow to Topock Gorg 
and Lake Havasu, and the movement into and through the Cibola Valley. Th 
river cross sectioning gives information on general river aggradation an 
degradation, and the detailed sampling gives information on the movement pas 


es 
4 “Computations of Total Sediment Discharge, Niobrara River near C : 
ody, Nebraska, 

U. S. Geol. Survey Water Supply, Paper No. 1357, 1955, 4 4 
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a specific point, including values that can be used in river shape-determinatio1 
equations and relationships. 

c. Bed material samples.—These samples have been collected in conjunctio! 
with the river cross sectioning and with the total transport-sampling measure: 
ments. Average bed-material-size distribution curves are shown in Fig. 5. 

d. River discharge records.—These records are available at several stra 
tegic locations throughout the river reach involved. The gaging station lo 
cations are shown on Figs. 2 and 3. 


For details on the sediment sampling, bed material sampling, river cros 
sectioning, and river discharge, reference is made to publications such as th 


TABLE 2.-SUMMARY OF AVERAGE YEARLY TOTAL SEDIMENT LOAD 
SINCE CLOSURE OF DAVIS DAM&;b 


Average Yearly Average Year]; 


rad Station Total Sediment Water Discharg 
Load (Tons) (Acre Ft)¢ 
Davis Dam R. S. 33 7,234,000 10,319,000 
to Below Needles 
Parker Dam Bridge 6,804,000 10,319,000 
R. S, 41 4,503,000 10,319,000 
R. S. 43 4,791,000 10,319,000 
Parker Dam Below Palo 
to Verde Weir 2,867,000 8,882,000 
Imperial Dam Taylor’s Ferry 3,646,000 8,882,000 


@ Based on flow-duration, sediment rating curve analysis. b Table 15 of “Interi 
Report, Total/Sediment Transport Program, Lower Colorado River Basin,” Janua 
1958. © Davis to Parker reach based on flow-duration curve for U.S.G.S. station ne 
Topock, Arizona, for calendar years 1950-56. Parker to Imperial reach based on flo 
duration curve for U.S.B.R. station estimated daily flows at Taylor’s Ferry for cale 
dar years 1951-56. 


ee 


U. S. Geological Survey Water supply papers, and reports on River Contr 
Work and Investigations.5,6 ; 


SCOPE OF CURRENT RIVER RECTIFICATION 


Until recent years, the river rectification work carried out was of an eme 
gency nature to alleviate the immediate problems. Examples are the channe 
ization below Needles in 1953, to lower the rising river level; channel cut-of 
in the Cibola Valley below Blythe, Calif., to lower river water surface at t 
Palo Verde Irrigation District outfall drain; and a cut-off above Yuma to mo 
the flow of the river away from an air-strip that was being damaged by ba 
erosion. Current and proposed river channelization and control work is nt 
more comprehensive in nature. Conditions between major channel structur 

Ss 3 


5 “Report of River Control Work and Investigations, Lower Colorado River Basi 
U. S. Bur, of Reclamation, Boulder City, Nev., November, 1957. a 

6 “Interim Report, Total Sediment Transport Program, Lower Colorado River Basi 
Sedimentation Sect., Hydrol. Branch, Bur. of Reclamation, Denver, Colo., January, 19 
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are being analyzed. The work plan between these structures includes an evalu- 
ation of all current and future conditions in that river reach. There are thret 
basic reaches involved in this comprehensive study approach: 


a. Davis Dam to Parker Dam (Lake Havasu) 
b. Parker Dam to Imperial Dam 
c. Imperial Dam to the International Boundary 


Note that the reach between Hoover Dam and Davis Dam (Lake Mohave) i: 
not included, since Davis Dam backs water to Hoover Dam and the whole reac! 
is inacanyon. For purpose of economic evaluation, the Parker to Imperia 
Dams reach has been broken in two subreaches (1) Parker to Palo Verde Dam 
(Fig. 2), and (2) Palo Verde to Imperial Dam (Fig. 3). 


AGGRADATION AND DEGRADATION 


The clear water released from the Hoover, Davis, and Parker Dams ha 
resulted in scour of the riverbed and bank materials in the channel below th 
dams, with resultant aggradation in the backwater reaches of the next struc 
ture downstream. General degradation andaggradation are illustratedin Fig. 6 

This scour and fill of the river alluvium results in a lowering of water sur 
face in one location and an increase in water surface in another location. I 
intermediate reaches an unstable, braided-type channel develops which allow 
the river to attack the vulnerable banks. Degradation introduces a continua 
supply of sediment to the downstream aggradation areas and results in deterio 
rationof riverbanks and bed. Aggradation results in a rise in river-water sur 
face and ground-water surface, growth of phreatophytes, and wide alluvial de 
posit areas. 

As has been discussed previously, the bulk of the sediment movement in th 
Lower Colorado River (since closure of Hoover Dam) is derived from the river 
bed and banks. This is important in the planning of river rectification wor 
because if this attack on the bed and banks can be brought under control, 
stable condition will eventually develop. This is the aim of the present rive 
control and rectification work, and planning program being carried out by th 
USBR. 


RIVER CONTROL AND RECTIFICATION PROGRAM 


Probably the best way to discuss the program for river control and rectif: 
cation in the Lower Colorado River is to beginat Hoover Dam and work dowr 
stream, examining each problem area in turn, pointing out the studies ar 
analysis made, the design considerations, construction program, and resuli 
obtained. The discussionis confined tothe major problems and problem area 
Many localized-type problems, associated with the operation and the river, al 
continually developing within the major problem reaches. 

Davis Dam to Lake Havasu.—With closure of Hoover Dam in 1935, and 
sequent closure of Davis Dam in 1950, a degradation cycle began in the rive 
downstream. Much of the material scoured from the river channel deposite 
in the lower Mohave Valley near Needles (Fig. 2). This lower valley was hi: 
torically an aggrading reach, caused in part by the constriction as the riv 
entered the Topock Gorge. The aggrading conditions were aggravated by tl 
closure of Parker Dam in 1938, and subsequent filling of Lake Havasu, whit 
created a backwater extending through the Topock Gorge. The aggradation 
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the area above Topock raised the river-water surface elevation and ground 
water elevation to a point that endangered the town of Needles. In 1946, emer 
gency legislation by the Congress enabled the USBR to purchase a dredge fo: 
the purpose of constructing a channel that would lower the water surface in th 
Needles area. The dredge purchased was a 20-in. cutter-head type and has | 
design capacity of 500 cu yd per hr and can move 1,300 cu yd per hr unde 
favorable conditions. The emergency channel constructed was 200 ft to 400 i 
wide, and banks were riprapped as necessary. This channel was completed i 
1953, and an immediate lowering in water surface of 5 ft was realized in th 
Needles area. 

Although this emergency channel work alleviated the situation, it was re 
alized that the effect would be only temporary because the source of materials 
the broad Upper Mohave Valley, was still contributing material tothe Needles 
Topock reach at anunabated rate. Further, the supply of material in the Uppe 
Mohave Valley, estimated at 112,000,000 cu yd, was sufficient to continue th 
sediment inflow to the lower valley for many years. The permanent solutio 
became obvious: The river had to be confined in a controlled channel fror 
Davis Dam to the mouth of Topock Gorge, a distance of about 43 river miles 
Actually, because the river is confined in a well armored channel from Davi 
Dam to Big Bend, a distance of about 10 miles, the channelization length re 
quirement was reduced to 33 river miles. To arrive at a comprehensive pla 
for a stable controlled channel, some specific information and data on existin 
conditions had to be obtained. Deep bed-material sampling was accomplishe 
by use of a “sand-bucket” and/or auger which were driven 10 ft maximu1 
depth or to refusal. Samples were retained for various depths and the vari 
ation in material size with depth was established. The total sediment load 
above Needles, at Needles, and in the Topock Gorge were determined by us 
of the Modified Einstein Procedure.4 A dominant discharge of 15,000 cfs wa 
determined from an analysis of Hoover and Davis Dam releases and a floodpe 
tential of 50,000 cfs from Davis Dam to Piute Wash, and 70,000 cfs from Piut 
Wash to Topock was computed for use in levee design. Determination of th 
dominant discharge was complicated to some extent by the daily fluctuation | 
releases from Davis Dam, which is used for power-peaking purposes. Th 
relatively high stage of Lake Havasu and its effect on the channel in the lowe 
part of the reach was an additional factor considered in the design. 

Since the sediment load in this reach of the river would, with channel cor 
trol, decrease with time from the present amount of about 18,000 tons daily 1 
a negligible amount, this factor had to be considered in the channel design. © 
is interesting to note that although the sediment concentration at Needles : 
only about 400 ppm(mostly fine to medium sand), the tonnage carried is hig 
because of the high discharge. When the channel has been established at de 
sign width, and banks have been riprapped as necessary, vertical degradatic 
within the controlled channel can be anticipated. This degradation will caus 
a siftingor sorting of the bed material until a sufficient layer of armoring m: 
terial is established. 

With the available data and the above factors in mind, the most compatib 
channel shape and alinement were established. Several available methods we: 
used in the channel design including relationships and equations developed | 
Leopold and Maddock, ’ Schoklitsch,8.and Lane.? The results of these relatio1 

7 “The Hydraulic Geometry of Stream Channels and Some Physiographic Implications 
by Luna B. Leopold and Thomas Maddock, Jr., Geol. Survey Professional Paper 252, 195 

8 «The Scholklitsch Bedload Formula,” by Samuel Shulits, Engineering, Lond 
England, June 28, 1935. : 

“Design of Stable Channels,” by E. W. Lane, Transactions, ASCE, Vol. 120, 1955 
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hips and equations as applied to the available data are shown on Fig. 7. Com- 
uted roughness values for the channel vary from 0.025 to 0.030. The amount 
f degradation that could be expected in the design channel width was deter- 
ined by use of the Schoklitsch bedload equation, competent bottom velocity 
1ethod, and tractive force criteria as applied to the deep bed material samp- 
ing results. Through the process of sifting and sorting of the bed material, 
he transportable sizes are scoured and the general channel gradient is flat- 
sned. Ata certain combination of bed material size, channel gradient, and 
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FIG, 7.—COLORADO RIVER CHANNEL STUDIES 


her hydraulic factors, stability is reached. It was determined that bed sta- 
lity would occur when a 6-in. armor layer of 4 to 10 mm material was ob- 
ined. Scour of 7 to 10 ft was determined forthe upper end of the controlled 
ach reducing to a few feet near Needles (Fig. 8). The final channel design 
ape, profile, and alinement is shown in Fig. 9. The overall plan, mnekidity 
ork accomplished to date, involves the dredging of over 20,000,000 cu yd of 
aterial and placing of about 236,000 cu ydof riprap. The dredged channel 
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ngth is about 30 miles, and is on a sinuous alinement with an average bend 
dius of 10,000 ft. Radii limitations are 5,700 ft minimum and 16,000 ft maxi- 
um. 

Construction of the design channel has now progressed to the lower end of 
ig Bend. The dredge is working upstream on the established alinement with 
“edge material being used for levee building as required. A channel width of 


NOTE 

99.9 Based on 1/4 cu. yd. max. Size __ 
and 1/10 cu. ft min. size!mean 
diameter computed from 4/3 1r3 
Graduation fluctuates with quarry 
type and operation. 


h rae 


‘Max. size limit 
(Approx. 28 in) 


FINCH 


FOANUCNI 
ow 
oO 


5% Rock fragments 
and dust allowable 


-Min. size limit 
(Approx. 7 in.) 


| 10 100 
MATERIAL SIZE - INCHES 
FIG, 11.—APPROXIMATE GRADATION RIPRAP, BIG BEND TO NEEDLES 


to 400 ft is cut depending upon the amount of material needed for levee con- 
uction. One pass ie the rede will create a channel 170 ft wide. Additional 
ening to design width is accomplished by river action. Riprap is placed on 
outside of the bends only, except in reaches with large radii or where other . 
ditions warrant riprapping on both banks. Generally, the riprap is placed 
he bank at design width and it falls in place as the channel widens. Riprap 
dation is shown in Fig. 11. The dredge is moving material at the rate of 
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14,000 cu-yd daily, working three shifts. Channel completion is scheduled fo 
early 1960. Recently, the proposed termination point for the channelization wa 
moved downstream a short distance to accommodate objections of local proper 
ty owners in the vicinity of Bullhead City. 

Because the sediment inflow to the Needles area will continue for man 
years, though at a decreasing rate, it was deemed advisable to acquire a small 
er dredge to maintainthe Needles to Topock channel while the upstream chan 
nel was being completed, and to maintain the completed channel infuture years 
A 12-in. cutter-head dredge, with a capacity of 300 cu yd per hr, was pur 
chased, in 1957, for this purpose and is currently operating in the area abov 
Topock Gorge. This has resulted in some stabilization of the river stage in th 
gorge and the water surface is being maintained fairly constant at Needles. | 
greater length of time will be required to determine the full results of thi 
program. 

During 1958, above-normal, sustained, high flows (20,000 to 25,000 cfs) wer 
experienced in the Davis to Lake Havasu reach. The designed channel, wher 
completed, functioned very well under these high flows. It is anticipated thi 
the most troublesome factors, under future operation, will probably be grav 
and cobble bars extending into the channel because of flash floods on tributar 
arroyos, and bar building at cross-overs within the sinuous channel. 

A resurvey of the active capacity of Lake Havasu was carried out in 195! 
and results indicate 58,100 acre-ft of capacity has been lost to sediment be 
tween operating elevations of 430 and 450. This represents 15% of the origin 
capacity within these elevations. This loss is primarily due to delta develoy 
ment at the upper end of the reservoir. Water is pumped from Lake Havas 
into the Colorado River Aqueduct for the Metropolitan Water District. 

With completion of the Davis to Topock channelization, a rapid reductic 
in the sediment transport of the river will be realized, thereby decreasing tl 
sediment inflow to Lake Havasu to a negligible amount and creating a controll 
water surface in the river at Needles. In addition to these benefits from tl 
channelization work, a significant amount of water salvage is accomplished | 
reducing water loss to phreatophytes and evaporation, since the wide shallc 
aggradational reaches are eliminated. Flood protectionfrom levee constructit 
and enhancement of recreational facilities are additional benefits. 

Data on river conditions, in the Davis Dam to Lake Havasu reach, will co! 
tinue to be collected to aid in evaluation of the design plan. When the riv 
channel has become stabilized, these data will aidin evaluating the various d 
sign equations which were used. Total sediment, load measurements and di 
charge records will help to evaluate the effectiveness of the plan in reduci 
the sediment load and increasing water salvage. 

Parker Dam to Palo Verde Diversion Dam.—No critical channel probler 
have developed between Parker Damand Palo Verde Diversion Dam—a distan 
of 65 river miles (Fig. 2). However, river cross sectioning and total sedime 
transport determinations indicate that a considerable tonnage of sediment 
being derived from this reach. Preliminary investigations have been carri 
out to determine the type of surface and deep bed material, and general riv 
channel conditions. Indications are that no immediate channel rectificati 
work is necessary, but it is evident that because of the sediment contributed 
this stretch of the river to downstream reaches, future work will probably hb 
come necessary. Headgate Rock Diversion Dam, whichsupplies water to Ind! 
lands, is located about 15 miles below Parker Dam. 4 
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Thefirst major problem to develop below Parker Dam was at the headworks 
ua the Palo Verde irrigation works which is located about 10 miles upstream 
rom the city of Blythe, Calif. Riverbed scour at the canal headworks con- 
inued until diversions could not be satisfactorily made. This lowering of the 
vater surface was arrested in 1945, by the construction of a temporary rock 
veir. Maintaining this weir proved to be a continuous problem. New rock had 
@ be added to replace rock moved on downstream during high river stages. 
Approximately 190,000 cu yd of rock were placed in the temporary weir between 
945 and 1956. Authorization was made by Congress, in 1954, for construction 
i a permanent structure a short distance downstream from the temporary 
veir. This structure, the Palo Verde Diversion Dam, was completed and was 
n operation during the 1958 irrigation season. The temporary rock weir was 
nodified in accordance with a requirement contained in the Congressional 
uthorization for the dam. This requirement, for nullifying the effects of the 
emporary weir, involved some detailed study to determine what cross section- 
i shape should be established, and how muchrock removal was required to ob- 
ain this section. Removal of the rock also constituted a unique problem since 
t involved removing material of 1 to 3 cu yd and smaller while the river dis- 
harge was about 7,000 cfs. The work was accomplished satisfactorily by the 
ontractor, however. 

Associated with the scour problem at the Palo Verde Diversion Headworks 
jas a sediment in-flow problem. Sediment diversion to the Palo Verde Canal 
as been heavy over the years—even before construction of upstream storage 
tructures. A desilting basin and dredge have been maintained by the Irri- 
ation District just below the canal heading. Because the new diversion dam 
ates and canal headworks can now be manipulated to minimize the sediment 
ntake to the canal, the desilting operations can be curtailed. 

Studies made inconnection with the permanent diversion structure included 
vater-surface profiles, upstream and downstream, and expected degradation 
1 the downstream channel. Degradation was determined by extending the slope 
rend that was occurring to an ultimate river slope of 1.0ft per mile. It was 
stimated that the bed would stabilize at the slope of 1.0 ft per mile. The up- 
tream, water-surface profiles were required to determine levee elevation 
long the Indian lands on the east side of the river. The downstream profiles 
nd degradation studies were required for tail-water determinations and struc- 
ire design. Ultimate channel degradation of 8 ft, beyond that existing at time 
f construction of the new structure, is expected. 

Palo Verde Diversion Dam to Imperial Diversion Dam.—The river channel 
roblems in this reach of the Lower Colorado River are currently being in- 
sstigated and studied in detail by the USBR. Problems include bank erosion 
ad channel deterioration between Palo Verde Diversion Dam and Ehrenburg 
ridge, near Blythe; channel aggradation with accompanying drainage, water 
'88, and channel capacity problems in the Cibola Valley, which extends from 
slow Taylor’s Ferry to the headwaters of Imperial Reservoir, andaggradation 
1d sediment inflow to Imperial Diversion Dam and headworks. An approxi- 
ate balance point between aggradation and degradation occurs near the Palo 
erde outfall drain in the Cibola Valley. Between Ehrenburg Bridge and 
aylor’s Ferry, the river is in a relatively stable reach which holds against 
e east side bluff and has a surface gravel bed in several locations. Fig. 3 
iows the location of these features. 

Current planning for this stretch of the river ihvolves channel straighten- 
gand bank protection in the upper reachand construction of a new channel in 
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the Cibola Valley. This channel-rectification plan is shown on Fig. 10. In the 
upper reach, some of the local land owners are placing protective works to 
stop bank scour. The primary objective of rectification in the Cibola Valley is 
to create a stable, controlled, river channel that will require a minimum of 
future maintenance. Other benefits, such as losses to nonproductive phreato- 
phytes, lowering of the water surface at Palo Verde outfall drain, and recre- 
ational and fish and wildlife improvement may result incident to this work. 
Severe drainage problems were reported in the lower Palo Verde Irrigation 
District lands because of the aggrading type channel in the Cibola Valley. As 
was the case upstream at Topock, there is a natural constriction of the river 
at the lower end of the Cibola Valley, near Adobe Ruin (Fig. 10). In 1947, a 
pilot cut-off was made just upstream from the outfall drain to ease the drainage 


situation. This cut was 1 5 miles long and 40 ft wide, and was accomplished 


by dragline operation. Additional widening was by river action. The work re- 
sulted in a lowering of water surface of 1 ft at the mouth of the out-fall drain. 
However, poor drainage and channel deterioration continued. It was concluded 
that major rectification in the Cibola Valley was needed. The first plan for 
such rectification was prepared in 1951.10 Since then, data acquisition has 
been carried out to aid in correctly evaluating the existing conditions and to aic 
in formation of a comprehensive plan for rectification. 

One of the most interesting and challenging aspects of the river planning 
and design in the Cibola Valley, is the proposed cut-off of about 9 miles ir 
length from Station 1863+98 to Adobe Ruins (Fig. 10). A cut-off in the alluvial 
Cibola Valley deposits must be designed to remain stable with minimum scou1 
upstream and minimum aggradation downstream under existing and future sedi- 
ment inflow and dominant channel discharge. A dominant discharge of 10,00( 
cfs has beenestablished from analysis of past and probablefuture river oper- 
ation. The sediment inflow to the reach is 10,000 tons daily, as determine 
from Modified Einstein computations at Taylor’s Ferry. The gravel bed ele- 
' vation below Taylor’s Ferry and the estimated maximum bed aggradation ele: 
vation at Adobe Ruin at the lower end of Cibola Valley must be met. 

For the existing bed material and sediment load, it was determined that : 
channel gradient of 1.2 to 1.3 ft per mile should be used to obtain stability 
Since the straight-line slope between stations would be considerably more tha 
the 1.2 to 1.3 ft per mile limit, it was evident that a sinuous type of cut-of 
channel had to be constructed to meet the slope limitations. Previous experi: 
ence by the USBR, the Corps of Engineers,1!1 and others, indicates that the bent 
radii should not be less than 5,000 ft for the hydraulic conditions in the sinuou: 
channel. 

The design width for the new channel in the Cibola Valley has been tenta 
tively determined as 450 ft maximum, by applying the acquired data and infor: 
mationto the Maddock-Leopold’ and tractive force? relationships. The result: 
were checked against data from actual field measurements at selected stabl 
reaches in the valley and agreed very well. As in the Needles area, the out 
side of the bends of the constructed sinuous channel would be riprapped. It i 
interesting to note that unlike the channelization in the Topock to Davis Dar 


— —q@“ 
10 “Report on Colorado River Channelization, Cibola Valley, Arizona,” U. 8, Bur, ¢ 
Reclamation, Colo, River Front Work and Levee System, Report No. RC-3-3.1, Boulde 
City, Nev., March, 1951. 
11 “A Laboratory Study of the Meandering of Alluvial Rivers,” by J. F. Friedkin, 'U. § 
Waterways Experiment Sta., U.S. Corps of Engineers, Vicksburg, Miss., May, 194: 
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each, there will be a continuous sediment inflow to the Cibola Valley until 
hannelization is carriedout in the Parker Dam to Palo Verde Diversion Dam 
each, (Fig. 2) which presently is not planned. A report is currently in prepa- 
ation which will present a comprehensive plan, and show benefits to be derived, 
y river rectification in the river reach between Palo Verde Diversion Dam 
nd Imperial Reservoir. 

Imperial Dam and Reservoir.—One of the most widely known and referred 
) desilting operations is the one at the intake to the All-American Canal at 
nperial Diversion Dam (Fig. 3). At the right side of Imperial Dam, the ma- 
rity of the river flow is diverted into the All-American Canal which has an 
itial capacity of 15,155 cfs. Water is also diverted at the left side of thedam 
) the Gila Gravity Main Canal which has an initial capacity of 2,200 cfs. Water 
assing the dam is for purposes of river regulation, which is discussed subse- 
iently, and for Mexican Treaty requirements. At the time of construction of 
nperial Dam, it was realized that large tonnages of sediment, primarily sand, 
ould be taken into the All-American Canal. Detailed studies, including labo- 
atory models were undertaken to devise a desilting works that would remove 
e sediment and return it to the river below the dam. The final desilting work 
lan consisted of six settling basins, arranged in pairs, equipped with twelve 
tating vanes or scrapers per basin that flushed the deposited material, as it 
attled, into sludge pipes. The sludge pipes dumped into the Imperial Dam 
luiceway channel.12 To date, the desilting has been-completely satisfactory, 
s illustrated in Fig. 12. Sediment concentrations at Station 60 on the All- 
merican Canal are a very small fraction of the concentrations entering the 
silting basins. 

Thereservoir above Imperial Dam, at thetime of closure in 1938, had a ca- 
acity of 85,000 acre-ft. The reservoir filled with sediment at a rapid rate 
ereby reducing the reservoir trap efficiency and allowing more sediment to 
2ach the All-American Canal headworks. The present reservoir capacity is 
aly 1,000 acre-ft. 

The sediment inflow will be reduced to some extent by bank protective works . 
1 the Palo Verde area andriver rectification in the Cibola Valley but will still 
2 a problem at Imperial Dam. 

In addition to the desilting works, sluicing operations are carried out peri- 
jically at Imperial Dam to remove sediment deposits upstream from the All- 
merican Canal headworks structure. Detailed sediment sampling and dis- 
large measurements are made during these sluices to aid in evaluating the 
eration, and for use in planning downstream river regulation. A desilting 
id sluice arrangement is also provided at the Gila Gravity Main Canal head- 
orks, but because of the relatively small amount of water diverted, the sedi- 
ent problem is minor at this point. ; 

The sediment problem will continue at Imperial Dam for many years to 
me because it is theterminal point for material carried downthe river from 
stream sources. Data collection and study will be continued by the USBR, to 
‘ep a check on conditions and to assist in planning of any future rectification 
ork, should such work become necessary. 

Imperial Dam to International Boundary.—Below Imperial Dam, the major 
‘oblem, at present, is that of river regulation. River regulation, as referred 

here, is the maintaining of a suitable channel that will transport the sedi- 
ents introduced at Imperial Dam without creating channel deterioration. This 


a Ee, A 
12 «Imperial Dam and Desilting Works,” Bulletin No. 6, Boulder Canyon Proj. Final 
ports, Part IV, Design and Constr., U. S. Bur. of Reclamation, Denvex, Colo., 1949. 
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3 accomplished by the release to the channel downstream of a certain per- 
entage of the water flowing into Imperial Reservoir. Water is released at 
oover Dam in accordance with a weekly master schedule based on require- 
ents of all water users. This schedule can be fluctuated within set limits to 
eet varying demands throughout the year. Based on studies of sediment data 
aliected at Imperial Dam and at Yuma, some relationships have been jestab- 
shed for use in predicting the amount of water needed for river regulation. 
ypical relationships are shown in Fig. 13. These curves are periodically re- 
ised to reflect latest data. The water used for river regulation can be the 
ame water released from Imperial Dam to meet Mexican Treaty requirements. 
hen the Imperial Irrigation District so requests, the water not required for 
ver regulation is put into the All-American Canal and returned to the river 
_the Pilot Knob Powerplant (Fig. 3). 

With the initial closure of Imperial Dam, degradation occurred downstream. 
his scouring actionreverted to aggradation as Imperial Reservoir filled with 
sdiment. Generally, the river channel has not yet recovered to the condition 
hich existed upon the Imperial Dam closure. With continued demands on the 
slorado River water and development of the water uses above Hoover Dam, 
e water available for river regulation below Imperial Dam will become less, 
id maintaining a suitable river channel will become more difficult. For this 
ason, new methods of silt control may become necessary. This situation was 
edicted in the studies and analyses made of the situation during the planning 
ases for Imperial Dam and desilting works.12 

It is interesting to note that the total sediment load now passing the Yuma 
ging station is only a fraction of that which passed the same point prior to 
nstruction of Imperial and Hoover Dams (Fig. 14). 

About 5 miles below Imperial Dam, there is an old irrigation diversion 
ructure, Laguna Dam, which was built in 1909 (Fig. 3). The only purpose 
is structure serves, at present, is to provide reregulation of releases from 
iperial Dam. Considerable sediment deposition has taken place in the wide 
ver area between the Imperial and Laguna structures. 

Below Laguna Dam, the interests of the USBR and those of the International 
undary and Water Commission are closely allied. Therefore, the problems 
the reach from Laguna Dam to San Luis, Mexico, which is the lower boun- 
ry in the limitrophe reach, will not be discussed in any detail. The USBR’s 
‘thest downstream river work has been the rehabilitation of the existing levee 
stem above and below Yuma.!3 This levee rehabilitation was necessary to 
ovide adequate protection against probable floods and any increase in water 
rface resulting from construction of Morelos Diversion Dam, below Yuma, 
Mexico. Part of the expense of this levee rehabilitation was borne by Mexico. 
e work consisted of placing 1,700,000 cu yd of levee embankment, including 
rap, as necessary over a total lengthof 41 miles. The levees were designed 
‘a flood of 140,000 cfs. This flood value was determined considering Painted 
ck Dam would control much of the Gila River flood potential. Riprap having 
sual thickness of 5 ft was placed at critical locations. Riprap limitations 
re a maximum size of 1/2 cu yd and a minimum of 1/10 cu ft in volume. 
Two interesting problems near the city of Yuma have resulted in rectifi- 
ion work by the USBR. At Yuma, another natural constriction of the Colo- 
io River exists (Fig. 3). Here, the river narrows from a wide valley into a 
‘ge-type opening about 300 ft in width. The area above this constriction is 


pa 
(3 “Brief History of Colorado River Levees Near Yuma, Arizona,” U.S. Dept. of 
rior, U. S. Bur. of Reclamation, Boulder City, Nev., November, 1955. 
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a natural depository for sediment, particularly during high flows when a back: 
water effect is created. The river course through this depositional area i: 
subject to considerable shifting. The river channel in the early 1950’s movet 
in against the left bank, attacking existing low-lying airport lands immediatel, 
north of Yuma. To correct this condition, a pilot cut-off was constructed t 
change the river course and give it a direct approach tothe constricted section 
This cut-off involved excavating 200,000 cu yd of material, placing riprap a 
the upstream end, and constructing a structure across the old channel to ob 
struct flow and force the river into the pilot cut. Shortly after completion o 
the cut-off in 1954, the entire river was flowing through the new channel. Thi 
relocated section of channel is continuing to function satisfactorily. 

Just below the constriction near Yuma, part of the city water supply is pumpe 
from the left bank of the river. During low river flows, considerable difficult 
has been experienced in maintaining an open channel from which the wate 
could be pumped. This difficulty resulted from the manner in which the low 
flow channel entered the constricted reach and from bar buildup where th 
water treatment plant effluent is discharged into the river. To alleviate thi 
situation, some deflection fencing was placed by the Bureau in 1957, unde 
special appropriation by Congress. Immediate improvement was noted, an 
the severity of the problem was reduced during 1958, by the relatively hig 
flows passing the Yuma area. 

Requirements for river rectification and control work in the Imperial Da1 
to International Boundary reach are likely tocontinue for many years to come 
Collection and analysis of data will continue, and a constant vigil will be key 
to foresee the development of critical conditions and thereby provide ample op 
portunity to take necessary action. Exchange of data and cooperation in ana 
lyzing conditions will continue with the International Boundary and Water Com 
mission as has been done in the past. 


SUMMARY 


The construction of Hoover Dam and other structures on the Lower Colorad 
River has essentially solved the flood problems, created a firm irrigation wate 
supply, and provided a large block of hydroelectric power. However, thes 
structures did not eliminate the problems of river control, but only change 
their character. The channel aggradation and degradation that have occurre 
with the construction of water diversion and storage facilities have created 
series of varied river control problems that require rectification. These rive 
alterations are continuing and rectification needs will continue as long < 
sediment-contributing reaches of the river remain uncontrolled. 

Rectification work, and planning for such work, is in various stages of de 
velopment. In the river reach between Davis Dam and Topock, channel dredgii 
is underway tocreate a newfully controlled channel. Work in this reach is 88 
completed (as of 1960). The channel has been designed to be compatible wi 
existing and future sediment load, bed material, discharge, and hydraulic fa 
tors. Inthe Palo Verde Diversion Dam to Imperial Diversion Dam reac 
studies are currently underway and a report will be made soon presenting 
plan for river rectification. The plan will probably involve bank control, cha: 
nel straightening and realinement, and a substantial cut-off channel in the C 
bola Valley. Downstream from Imperial Dam, sediment and hydraulic data a 
being collected, and studies andanalyses are being carried out for the purpo 
of determining river regulation water requirements and to provide a consta 
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eck on river channel developments. Various items of river rectification 
wk are carried out as required within these major river reaches to correct 
localized channel problem. 

Sediment movement and hydraulic data are collected ona regular basis at 
rategic locations throughout the Lower Colorado River Basin. Results and 
alysis of these data are published by the USBR.°>§ Total sediment transport 
mpling stations are established or dropped as required by the detailed study 
ogram. 

Since the river channel aggradation and degradation results from a shuffling 
d reshuffling of riverbed and bank materials above and below major struc- 
res, planning and design for river rectification is aimedat obtaining eventual 
mplete control between these major control points. Since any work must be 
stified on the basis of benefits to be received, it is not possible to carry out a 
mplete construction program between major structures, but, rather, work 
ist be concentrated on reaches or sections where benefits will justify the 
ctification costs. This intermediate type work is planned as an integral part 
the overall river control program, however. 

River rectification work generally involves channel straightening and re- 
nement, utilizing cut-offs, riprap placement, training structures, and levees. 
ie method used, in each case, is selectedon the basis of existing and probable 
ure conditions. Improvements in methods of control are constantly being 
ide from the experience gained on work accomplished so far, and the analy- 
3 Of newly collected sediment, hydraulic, and hydrologic data. Present indi- 
tions are that future needs for river rectification will be in the Palo Verde 
version Dam to International Boundary reach. The most urgent need at 
esent is in the Cibola Valley. 

Looking into the future, it is conceivable that the Colorado River will even- 
uly become fully controlled from Hoover Dam tothe International Boundary. 
th increased population trends in this region of the United States, it can be 
sected that the water resource, including the recreational potential along 
‘se 300 miles of river, will be utilized at an ever-increasing rate. A fully 
itrolledriver will allow the most complete use of the available water for the 
mbined benefit of all interests in the Lower Colorado River Basin. 
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FRICTION FACTORS IN CORRUGATED METAL PIPE2 


Discussion by W. O. Ree 


W. O. REE,/ F. ASCE,—Data presented in this paper provided an oppor- 
ity for the writer to test a relationship he had used to estimate Manning’s 
jalues for paved invert corrugated pipe. In 1953 the writer tested a full- 
ited, paved-invert, corrugated pipe of 18-in. diameter. Tests were made 
various depths of channel flow as well as for the full-pipe condition, Fric- 
n factors were obtained for this pipe but the problem of estimating n-values 
other sizes of paved invert corrugated pipe still remained. It was beyond 
iilable resources to test other sizes so resort was made to an estimating 
ae adapted from a method devised for open channels by Robert E. Hor- 


Manning’s formula was used in this analysis, The author’s symbols will be 
ployed in the discussion plus a few additional ones. These are: 


P = wetted perimeter, in feet; 
p = subscript indicated paved; and 
¢ = subscript indicating unpaved (corrugated). 


nning’s formula is: 


ae 1486 R2/3 gl/2 


stituting = for R and rearranging 


passe no) 3/2 _ 3/2 P 
V 


\ pair of values Np and ng will be assumed to satisfy this equation, thus: 


1.486 s1/2\3/2 3/2 Pp 3/2 Pe 
an eaees = Np + Ne A 


Se tember, 1959, by Marvin J. Webster and Laurence R. Metcalf. : 

aN Supervisor, Hydr. Lab., Agric. Research Service, USDA, in cooperation with 
Okla. Agric. Experiment Sta., Stillwater, Okla. 

“Separate Roughness Coefficients for Channel Bottom and Sides,” by Robert E. 
on, Engineering News-Record, Vol. II, November 30, 1933, pp. 652-653. 
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equating the equivalents, 


3 
n 
n3/2 5 PP, 


solving for n, 


(*/? Py + ngi/2 Be) i? 
shies Vea aaaaiae sien am 


This is the prediction equation. The data from the writer’s experime 
were used to test it. The experiment will be described briefly and the pert 
nent data tabulated. 

The pipe was part of a full-scale pipe outlet spillway set up for testing 
the Stillwater Hydraulic Laboratory of the Agricultural Research Service a 
the Oklahoma Agricultural Experiment Station. The pipe was new, corrugate 
full coated, and had a paved invert covering 38% of the perimeter measured | 
the minimum diameter. The test section was 200 ft long and was installed | 
a 3.57% slope. Depths of flow were measured at 20-ft stations by open colur 
glass manometers connected to piezometer openings in the invert of the pig 
Twenty-six channel flows were run with depths ranging from 0.17 to 1.30. 
Velocities were high and depths were small so fluctuation and variation we 
experienced in the depth determinations. The depth measurements were ave 
aged for the ten stations and this smoothed the data. Even so, the sixteen te 
flows contained within the paved portion of the pipe yieldedvalues of Manning 
n for the paving ranging from 0.0085 to 0.0112. Casting out the three extre1 
values reduced the range from 0.0090 to 0.0099 with a mean value of 0.00 
The ten remaining flows, which were deeper, yielded Manning’s n-values whi 
showed a more consistent behavior. The hydraulic elements for these tes 
are given in Table 1. 

The Manning’s n-values in Table Rl are plotted against the correspondi 
depth values on Fig. R1(a). Onthis same figure is shown the calculated n-dey 
relationship. The agreement between the estimated and actual values of Ma 
ning’s n is fairly good, being very good at full pipe flow. This agreement ga 
some encouragement as to the usefulness of Eq. 1. However, the writer wv 
unwilling to give credence to the formula, because of boththe nature of itsoz 
gin and the quality of the data used inthe proof. One uncertainty factor, hithe 
to unmentioned, was the value of nc. No data were available on coated pipe 
the value of 0.025 for plain corrugated pipe reported by Straub? was used int 
estimate of the composite n. Another matter for concern was that the te; 
were in the super-critical flow range. In view of these uncertainties the ic 
was allowed to lie dormant until the data reported by Webster and Metcalf we 
used to test the prediction equation. ‘ 

One difficulty in the application of the Webster-Metcalf data was the lac’ 
an experimental determination of ny (Manning’s nfor the paving alone). A vz 
of np was calculated by applying Eq. 1 to the lowest flow in the 5-ft diam 
25% paved invert pipe. The result of this determination agreed exactly with 


3 “Hydraulic Tests on Corrugated Metal Culvert Pipes,” by L. G. Straub and H. 
Morris, St. Anthony Falls Hydr. Lab., Technical Paper No. 5, Series B, 1950. 
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authors’ assumed value of 0.011 forthe paving. This value was used in the su 
sequent calculations. 

The value of nc used in the calculations for the 5-ft diameter pipe was t 
average of the values obtained in the channel flow tests on the corrugated pi 
without paving. This value was 0.0244. For the 7-ft diameter pipe no channe 
flow data were reported so a mid-range value of nc = 0.0237 was used in tl 
case. 

The results of the comparison between the predicted values of n andthe a 
tual values of n for the 5-ft diameter pipes are shown on Fig. R1(b) and ( 
The agreement is good. 

For the 7-ft diameter pipe with the 25% paved invert the only comparison tk 
can be made is for the full pipe condition. A calculation of Manning’s n yield 


TABLE R1.—SUMMARY OF TEST RESULTS OPEN CHANNEL FLOW 


D R A Vv s 
(FT)4 (FT) (SQFT) | (FPS) (%) 


ooooocoooccocoo 
OMT DS NS SOS cE aE Pine ETA 
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oooooooococs 
ae er eT S 


iii Tae aa eee ae ae See 
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° 


4D is depth of flow. 
b Temperature estimated from reading taken on nearest day. 


an estimated value of 0.0209. This agrees well with the actual values o 
0.0206 to 0.0211) obtained inthe range of Reynold’s number'compatible with - 
selection of 0.0237 for ne. : 

The good agreement between the estimated and actual values for Mannin 
n indicates that the prediction equation has some merit. It should be useful 1 
making n-value estimates for other composite pipes within the range of c¢ 
ditions reported.herein. 

The writer wishes to express his appreciation to the authors for mak 
these pipe-flow data available to the profession and for giving him an opp 
tunity to test an idea. ; 
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EDDY DIFFUSION IN HOMOGENEOUS TURBULENCE? 
Discussion by Mikio Hino, Takashi Ichiye and Charles G. Gunnerson 


MIKIO HINO.!—The author’s method of determining the turbulent charac- 
istics in channel flows is very simple and skillful, and leads to the success- 
verification of the Kolmogoroff similarity principle. 

dowever, this relationship (proved by experiments) can also be derived by 
theoretical consideration utilizing the well-known statistical characteris- 
| of the two-dimensional channel flows. Moreover, the constant in Eq. 21 
be shown to be dependent on the roughness of the channel bed. 

first of all, the mean sizes of eddies, L, for the fully developed pipe or 
nnel flows, are determined only by the ratio y/h, where y and h represent 
pectively the distance from the wall and the radious of pipe or the depth of 
channel, The reasons are as follows: As the experiment by J. Laufer2 
proved, the production rate of turbulence and the dissipation of energy are 
rly balanced in the range of approximately y/h <1/2, and the turbulent in- 
sities referred to the friction velocity U,, Vu2/U, etc, are functions of y/h. 
The production rate is given by 


dUigs Sake ( eres 3 
Production = T =— seer 1 | UR arr rar vee (24) 
sidering the fact that /v2 and /W2 arealmost proportional to U2, and as- 
rae 

vL 
nt flows is also applicable approximately for the flows considered, 


= constant, well established for isotropic tur- 


ing the relationship, x 


3 
Rene ey: 


Dissipation = A'y ites a L 


refore, by equating Eqs. 24 and 25, we have 4 


Les 7) 


September, 1959, by G. T. Orlob. 
i Graduate Student, Dept. of Civ. Engrg., Faculty of Engrg., Univ. of Tokyo, 


= fy (E):: Separate) 


is Structure of Turbulence in Fully Developed Pipe Flow,” by J. Laufer, NACA, 
i. 1954. 
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The experiment by J. Laufer? shows that Eq. 26 holds for the whole fle 
range. The relationship given by Eq. 26 is also assumable from the expe! 
mental fact that the Prandtl mixing length divided by h, l/h, is expressed as 
function of y/h only and independent of the wall roughness and the Reynol 


number. 
Therefore, the rate of energy dissipation, Eq. 22, can be modified as 


E = Uo € Se = fo(£) U2 ty orl, (: 


‘where Uy expresses the mean velocity of the flow. 
00 


Secondly, since the decaying time, te Rédé, is proportional to L/Vu2 
0 


LAV Ww, the diffusion rate defined by Taylor, Dz, is represented by 


Dz = we f Rede =cw ath 
oO u2 


ng 
— 


- fg (Z) vy 


Thirdly, the Lagrangian eddy size La, is related to L as follows, 


where U represents the mean flow velocity at the point. 
Lastly, combining these equations and putting y=h, we have 
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where 


yp Te; 


Ux 


2? k = Karman constant 
By 


Cr = 


a PTE sy ae 


U,; =Uat y=h. 


The relationship now obtained shows the effect of the bottom roughn 
which is not great. However, as the author anticipated, the effect will be 
apparent when the wide ranges of the bottom roughness are tested. “4 


3 “Investigation of Turbulent Flow in a Two Dimensional Channel,” by J . 
NACA, TN 2123, 1950. 


(4 DISCUSSION 97 


To check the validity of Eqs. 26 and 29, one needs the experimental results 
Iwing the relationships between U, Se, h, and La. 


TAKASHI ICHIYE.4,b—In the introduction the author described various kinds 
lefinitions of eddy diffusivity from an historical point of view, giving the im- 
2ssion that they are based on different concepts on diffusion. Instead, these 
initions can be derived in a unified manner from a statistical theory of dis- 
rsion of a particle or particles which is not necessarily pertinent to the 
ory of turbulence but is essentially based on the theory of Brownian motion 
random walk. Such an approach was initiated by Batchelor,5 who derived the 
inition represented by the Eqs. 1, 2, and 4 from the theory of “one particle 
lysis” and the neighbor diffusivity in Eq. 6 from “two-particle analysis” of 
tchelor and Townsend.6 As for the latter the validity of Eq. 7 over a wide 
ige of scales of turbulence and the theoretical explanation have been already 
en elsewhere. 7,8 

The diffusion of one particle in the turbulent field has been discussed by 
hors of some standard text books differently from the present paper.6,9,10 
¢ details can be obtained from these books, but it is worthwhile to repro- 
-e the discussion for comparison with the present paper. 

The standard deviation of the z-coordinate (the direction of which is arbi- 
ry) oz or Z2 of the particle which is at the origin at the initial time is given 


which w(t) is the turbulent velocity. (The notations of the paper are pre- 
ved as far as possible, though they are sometimes different from the con- 
tional ones.) 

By introducing the correlation function R(€) (or R&) this equation is written 


pie 
aoes ad Bismeragiaty diye: nape (b) 


Florida State University, Oceanographic Institute Contribution No. XXX. 

- Asst. Prof., The Oceanographic Inst., Florida State Univ., Tallahassee, Fla. Last 
The Application of the Similarity Theory of Turbulence to Atmospheric Diffusion, 
3. K. Batchelor, Quarterly Journal, Royal Meteorological Soc., Vol. 76, 1959, pp. 
-146. ; 
'Turbulent Diffusion,” by G. K. Batchelor and A. A. Townsend, Surveys in Me- 
uics, Cambridge University Press, Cambridge, England, 1956, pp. 352 399. 


“Horizontal Diffusion,” by F.C. W. Olson and T. Ichiye, Science, Vol. 130, No. 


L, 1959, p. 1255. 
“On Neighbor Diffusivity in the Ocean,” by T. Ichiye and F.C. W. Olson, Deutsche 


ngeraphische Zeitschrift, (in press). ; 

Turbulent Diffusion: Mean Concentration Distribution in a Flow of Homogeneous 
sulence,” by F.N. Frenkiel, Advances in Applied Mechanics, Academic Press, N.Y., 
@, 1953... ; 
ee Flow Theory: II. Turbulent Flow,” by S-I. Pai and D. Van Nostrand, 
ceton, N.J., 1957. 
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In this derivation we have assumed the condition of homogeneity that the co: 
relation w(tt+é) w(t) is a function of only. After some manipulation of tl 
double integral this is reduced to 


in which we have assumed no more properties of R& than being an even funt 
tion about 

From general properties of the correlation function RE two special cast 
of (c) can be derived immediately. First, since R& should vanish for a lar; 
value of &, we can define Lagrangian time scale of turbulence Lt as 


when t >> Lt, Eq. c is written as 


m= 2wityt-2w fg Ry dé AoE, ( 


The last term is constant and thus only the first term in the right hand si 
becomes important for a large value of —. The comparison of Eq. e witht 
probability density of the particle computed from the solution of the Ficki 
equation gives the eddy diffusivity equal to wa Lt. Therefore the definition 
the eddy diffusivity compatible with the Fickian diffusion or Brownian moti 
is possible only for the case that t >> Lt. This definition leads to the ed 
diffusivity which is constant in time and space. In some problems of meteo 
ology and oceanography eddy diffusivity which is variable intime and/or spa 
has been introduced in order to explain distributions of some properties. Hor 
ever, such an approach will not give any information about the relationshi 
between the mode of diffusion and physical nature of the turbulence. Accor 
ingly, the discussion of the Eq. 10 through 13 seems to be meaningless, b 
cause they just express the Eq. c in different ways by introducing a ia 
definition of eddy diffusivity. 

When t is small compared with Lt, the correlation Function R€ can be wr: 
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in which a characteristic time A is defined by 
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ul Bes (a)"/ eis. (aR ¢/a8?) epee uta: (h) 


hen t <A, Eq. g may be written as 


The curves of oz-t (or oz - x as in the present experiment)can be obtained 
latively easily from the diffusion pattern of dye or discrete particles. Eq. e 
jicates that the eddy diffusivity D = Dz(o) = w2 Lt canbe determined from the 
ope of oz - t curve fort >> Lt. Eq. i may be used for estimating w2 from 
= same Curve in the neighborhood of t= 0. Then the value of Lo=U Lt can 

obtained from D and wé thus determined. 

This process of determining the eddy diffusivity does not need any informa- 
mn or assumption about a functional form of the correlation function R(é) as 
the case in the paper. Instead, some authors derived the functional form of 
—) from the curves of oz - t obtained by experiments,9 though direct mea— 
rements with hot wire technique give more accurate results. Actually the 
i¢tional form assumed in Eq. 15 is not valid near & = 0, because the gradient 
(€)/d becomes discontinuous at § = 0. The casual inspection of the picture 
fluorescein dye pattern in the paper reveals that the spread of the dye is not 
actly expressed by the Eq. 16 near the origin. 

The isotropy of the turbulence is not confirmed yet in a shear flow which is 
asic flow in the present experimént. Even when the turbulence is isotropic, 
> apparent coefficient of eddy diffusion determined by the procedure similar 
the one described above or used in the paper is not so. The longitudinal 
efficient was found to be much larger than the lateral one in the flow of a 
ye11 and in an open channel.12 It seems tobe worthwhile to compare the re- 
it of the present experiment with those of Elder,12 since both were obtained 
an open channel. 

Representation of Kolmogoroff’s similarity principle may be done in several 
ys. In a two-particle analysis the relation similar to Eq. 21 has a well- 
inded background for a representation of this principle, when neighbor dif- 
ivity F and initial separation lo are taken instead of Dz(o) and La. (See 
tchelor® and Ichiye and Olson).8 On the other hand, in a one-particle prob- 
a as in the present paper, Eq. 21 expresses only a dimensional relation and 
physical meaning is obscure, This kind of relation is valid only in a range 
large wave numbers where the dissipation of turbulence is caused by vis- 
is forces (Batchelor’s equilibrium range). Inthis range the similarity princi- 

is expressed by 


ei/3 4/3 AP Pine Ore eee nee ee epen td) 


which V is the molecular viscosity and ¢€ is dissipation rate of energy. The 
racteristic length 7 is calledKolmogoroff’s microscale which corresponds 
‘roximately to the minimum wave number of the equilibrium range. 


ll “The Dispersion of Matter in Turbulent Flow through a Pipe,” by G.I. Taylor, 


ee edings, the Royal Society, Vol. A223, 1954, pp. 446 468. 
*The Dispersion of Marked Fluid in Turbulent Shear Flow,” by J : W. Elder, 


rhal of Fluid Mechanics, Vol. 5, No. 4, 1959, pp. 544-560. 
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The other difficulty in the application of Eq. 21 is the arbitrary definition + 
E. In a shear flow this quantity should include frictional velocity or shearir 
stress at the wall of a channel which in turn includes the molecular viscosit 
These factors may be implicitely included in Se of Eq. 22. But then the cor 
stant of Eq. 21 loses its nature as a universal constant. Aside from thes 
theoretical difficulties, the scatter of experimental points of Dz against E? L; 
in Fig. 11 seems to be too large to ensure the relation of Eq. 21, though Fig. 
may represent the so called 4/3-power law well. 

From the energy spectrum theory we can construct another relation be 
tween the eddy diffusivity and the scale of turbulence, which cannot necessari 
be claimed to be unique. In an inertial subrangel3 where the dissipation o¢ 
curs by inertia terms, the energy spectrum E(n) (in which n is a wave numbe 
is given by 


EGGiielts -2/3 1/3 


(Batchelor’s13 equation [6.5.1], hereafter indicated by [ B6.5.1]). 


Put no the wave number of the eddies which contain most parts of energ 
Then at no, E(n) has a maximal value 


The energy spectrum for a range of small wave numbers may be different 
different flow system but since we are concerned only withthe integrated for 
of this spectrum the contribution from this range is small. Thus we may pu 


EJ= Ep, (n/n) for,” 0.) mong kee ee (r 


in which m is a positive constant. Under the isotropic condition 


we = ; “f E(n)*dtit sn. 3) See ( 
0 F 
Substituting the value of E(n) of Eqs. k and m into Eq. n we have 5 
ES ey £7 ene 25.949 1 
we =B€E = al ce 
B ng 1c Pag toc nt mal Sykc ce f 


The correlation function w(t+€) w(t) is the same as the second equation 
(B3.4.4). Then the combination of equations (B.3.4.9) and (B.3.4.21) gives 


one Rg dé = 7 (4u)-1 "i n-1 E(n) dri + isi; aeeuee 
0 


13 “The Theory of Homogeneous Turbulence 


” by G. K. Batchelor, Cambri 6 [ 
versity Press, Cambridge, England, 1953. ; - — 
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which U is the mean current introduced in order to change the original dis- 
wee integral into time integral of the present definition of eddy diffusivity. 
ain substituting Eq. k and m into Eq. p we have 


: ! 1 
D, (0) = y U7! €2/3 ng 5/3 y = 47a (m- +3) 
e scale of turbulence Ly can be obtained by 


Lig = UD x) fw2)7) = yg bots... (r) 


om Eqs. q and 3 the similarity principle is expressed by 
ieee eect ret Sar is hg) oe (s) 


In a shear flow the mean velocity U may be related with €, but it is impor- 
it that the eddy diffusivity defined from Eq. e is dependent on 5/3 power of 
» size of the energy-containing eddies inaflow with a constant mean velocity. 
The results of the experiments are quite important and interesting, though 
erpretation and arrangement of the data may be revised. 

Acknowledgment: The author is indebted to Mr. C. G. Gunnerson for his 
yice in preparation of the manuscripts. 


CHARLES G. GUNNERSON,!4 F. ASCE.—The author’s extension of the 
our-Thirds Law” shown on his Fig. 8 and 12 is an important contribution. 
3 statement concerning the significance of eddy diffusion in the problem of 
ste assimilation by the receiving water of a large coastal bay immediately 
sgests the discharge of Los Angeles’ Hyperion Treatment Plant effluent into 
ita Monica Bay.15 

There are several difficulties in applying the eddy diffusion equation: 


Rare TSA ge TR eee rete (1) 


ich is a restatement of the author’s Eq. 7 and of the trend equations shown 
Fig. 12. As the author indicates, there are a variety of definitions of the 
ile of the phenomena, L. The trend lines on Fig. 12 indicate a range of about 
) orders of magnitude for the correlation coefficient, e. Finally, it is ap- 
‘ent from the literature cited by the author and others16,17,18 that there is 
universally accepted characterization of horizontal diffusion phenomenaand, 


14 Civ. Engr., Bur. of Sanitation, City of Los Angeles, Calif. ; 
15 *Sewage Disposal in Santa Monica Bay,” by C. G. Gunnerson, Transactions, ASCE, 


. 124, 1959, pp. 823-851. 
L6 “Diffusion in Bikini Lagoon,” by W. H. Munk, G. C. Ewing, and R. R. Revelle, 


nsactions, Amer. Geophysical Union, Vol. 30, No. 1, 1949, pp. 59-66. 

i7 “Horizontal Diffusion in the Sea,” by J. Joseph and H. Sender, Deutsche Hydro- 
phische Zeitschrift, Vol. 2, No. 2, 1958, pp. 49-77, (translation by G. I. Roden). 

8 “Diffusion of Sewage Effluent in an Ocean Current,” by N. H. Brooks, Proceed- 
2 First Internl. Conf, on Waste Disposal in the Marine Environment, Univ. of Cali- 
iia, Berkeley, Calif. (in press). 
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therefore, the definition of the coefficient of eddy diffusivity, k. Addition: 
difficulties lie in making adequate observations. 

A logical method for predicting the behavior of a waste discharge involve 
the introduction of a suitable tracer intothe receiving water. The choice of tI 
tracer depends on local conditions.19 Fluorescein dye has been used in San 
Monica Bay and has been released either at a constant rate or intermittentl 
The latter provides a convenient time scale. Additional information has be¢ 


FIG, Gl1.—DYE DIFFUSION STUDY IN SANTA MONICA BAY. NOTE APPARENT LA 
OF RELATION BETWEEN DIFFUSION, CURRENTS AND WAVES, 


obtained from the trajectories and dispersion patterns of the more than 50 
drift cards which were released from various stations in 1955-56. The gene 


al agreement of the results of these diffusion studies with other work has be 
previously reported.20 ' 


19 “Tracer Methodology and Pollutional Analyses of Estuaries,” by E. A. Pears 
Proceedings, First Internl. Conf. on Waste Disposal in the Marine Environment, Un 
of A Berkeley, Calif. (in press). . 

20 “Studies on Eddy Diffusion in Santa Monica Bay, California,” by R. B. nibs 
C. G. Gunnerson, Transactions, Amer. Geophysical Union, Vol. 39, No. 3, 1958, p. 5 
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A comparison of the dye stream in the author’s Fig. 5 with photographs of a 


ch larger dye stream in the ocean reveals the range of variation which ap- 
es to diffusion phenomena. 


FIG, G2.—DYE DIFFUSION STUDY. EACH SLUG IS ABOUT 350 FT LONG, 


he photographs were taken from an airplane of a diffusion experiment in 
a Monica Bay.21 A 10% solution of fluorescein was discharged from a skiff 
iored about 5 miles off shore. The dye was fed at a 12 gal per hr rate and 


“Summary Report, Oceanographic Investigation of Santa Monica Bay,” Bur. of 
ation, Los Angeles, Calif., 1956. 
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released intermittently at a depth of about 4 ft on a schedule of 10 min-ON an 
5 min-OFF. 

The experiment was carried onfor about 7 hr. During this time, the curre 
underwent one rapid right-angle change in direction. However, the velocity re 


FIG, G3.—DYE DIFFUSION STUDY. THE DYE PERSISTED FOR OVER 3 HR, . 


mained essentially constant at about 1/3 knot. Accordingly, the dye “slu : 
were about 350 ft in length. The widths of the slugs were obtained from ve 
tical photographs which included 2 boats, the positions of which were kno ‘ 
No measurements were taken within 500 feet of the source. | 

Fig. G1 shows the dye pattern about 1 hr after the abrupt 90° change in cu 
rent direction had taken place. Neither the current shift nor the lateral dif! 
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mn of the dye appeared to be related to the wind waves or swell. 

Pig. G2 shows the variety of patterns which were found. Some of the varia- 
ms may be due to the boat’s swinging on its anchor, 

Fig. G3 is an oblique taken during the time when the dye stream was rela- 
ély stable and when the measurements of lateral diffusion were made. The 
ht-colored area just beyond the zone where the dye disappeared is the sea- 


k, cm’ sec 


Scale, L, centimeters 


_ FIG. G4.—RESULTS OF VARIOUS INVESTIGATORS’ DETERMINATION OF 
: CORRELATION IN EDDY DIFFUSION EQUATION, k = eL4/3, 


‘d edge of the inshore water mass associated with the discharge of Hyperion 
uent and of cooling water from nearby steam plants,15 

The widths of the stream determined from vertical photographs were 
90thed and used for the diffusion calculations and were assumed to include 
ut 95% of the stream. ' 

The detailed results of these diffusion studies have been previously report- 
5 and are included in the author’s Fig. 12. Taken by themselves, those re- 
§ are reasonably consistent. 
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It should be noted that Fig.12 is typical of a number that have appeare 
since Inoue22 in 1950 first assembled an array of published coefficients of edd 
diffusivity. He found the correlation coefficient in c.g.s. units for the writer’ 
Eq. 1, to be, e = 0.01, although the data showed a scatter of as much as thre 
orders of magnitude. 

Several subsequent investigators have fitted their results to Inoue’s curv 
and found reasonable agreement.23,24 A few, including the author and Olso 
and Ichiye,25 have determined those curves (more precisely, those values 
e) which best characterize selected data. This selection is important, for: 
determines the validity of the final formulation. 

A partial summary of the published curves showing the relation of the co 
efficient of eddy diffusivity to the scale of the phenomena is shown on Fig. G4 
The light dashed lines surrounding the curves of the various investigators in 
dicate the extent and scatter of their data. The coefficients (e) are in metri 
units, although conversion to English units may be made from the secondar 
scales. 

Some new data are included. The current cross pair value is derived fro! 
observations made during bacteriological testing of Hyperion effluent. Th 
field was tagged by both dye and current crosses whose center of gravity wa 
at about the 4-ft depth.21 The 1-mile outfall field value is based on the sprea 
of the field from the 1-mile outfall at Hyperios.15 The drift card values ar 
based upon analyzing drift card returns reported by Stevenson, Tibby, ar 
Gorsline,26 using data from groups of stations which showed shore recoverie 
essentially normal to assumed streamlines. The data on radioactivity in Bikil 
Lagoon!6 are included to show their similarity to drift card observations. 

All of these coefficients of eddy diffusivity were determined by means of tt 
formulation of Munk, Ewing, and Revelle16 


where k is the horizontal coefficient of eddy diffusivity in cm2/sec, and oj ar 
o2 are the valuesof the standard deviation ofthe concentration ofatracer alor 
a line perpendicular to the flow at ty and tg, respectively. 

The measured widths of the various streams (dye, current crosses, Hyperic 
effluent) were assumed to contain about 95%of the tracer. Thus, the measure 
width, u equals 4 0. 

Substituting, Eq. 2 becomes 


22 “The Application of Turbulence Theory to Oceanography,” by E. Inoue, (in Jap 
nese), Journal, Meteorological Soc. of Japan, Vol. 28, No. 11, 1950. 

23 "éOn the Eddy Diffusion of Punices Ejected from Myojin Reef in the Southern si 
of Japan,” by H. Hanzawa, Records of Oceanographic Works in Japan, Vol. 1, No. 
1953, pp. 18-22. | 

“Discussion, The Measurement and Calculation of Stream Reaeration Ratio,” F 
E. A. Pearson, Proceedings, Seminar on Oxygen Relationships in Streams, - ‘Taft Sa 
Engrg. Center, USPHS, Cincinnati, Ohio, 1958, pp. 43-45. . 


5 

25 “Horizontal Diffusion,” by F.C. W. Olson and T. Ichiye, Science, Vol. 130, 195 

p. 1255, 

26 “The Oceanography of Santa Monica Bay, California,” by R. E. Stevenson, RA 

Tibby, and D. S. Gorsline, Allan Hancock Foundation, Univ. of Southern California, L 

Angeles, Calif., 1956; Fig. 82, Stas. 3430, 3436, 3437; Fig. 89, Stas, 3512, 3513, 151 
3519; Fig. 95, Stas. 3940, 3945, 3946, 3951; Fig. 96, Stas. 4025- ‘4032, incl. 
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sand? (bgt) SE Oe ilies ane ire eh = iy (3) 


The scale of the phenomena is taken as: 


Wi, +W 
patil nese (port 620. ont 0 


It is significant that for lateral diffusion where no boundary restraints ex- 
ed, the curve k = 0.005 L4/3 gives a reasonable fit. Where the flow was 


10%, 


of: 


Serer 
*” jdneaas 


(it) 10 10° (1oonr-1odO* 


lat (t,-t,) = 
= | mile 


24 hours 


METI TerlS 


t, seconds 


FIG, G5.—SOLUTION OF EDDY DIFFUSION EQUATIONS, 


g,2- o 
k = eL44 = 1 


Pe EE TEXT. 
2 (tg - ty) © = 
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bounded on one side by the shoreline, the curve k = 0.0026 L4/3 more accu 
rately describes the data. 

Reference to Fig. 4 shows the wide variation in the reported values of th 
coefficient, e, of Eq. 1. As noted above, the lower values may reflect boundar 
conditions. The explanation of the higher values determined by the author ma 
lie in the difference in nature of the turbulence in his laboratory channel an 
that turbulence in the stratified ocean. Undoubtedly, some of the difference 
are due to the variety of approaches which have been employed. 

In any event, the importance of eddy diffusion in waste disposal problem 
lies in the feasibility of predicting the spread or dilution of a pollutant fro 
this phenomenon. It is appropriate to examine the effects of assuming variou 
values of e in such a prediction. 

From the writer’s Eq. 1, 2, and 4 


2 2 
: (ei + Wo 4/3 e Wo” - Wy F 
2 oe 32 (to-t1) a'@ £6 aC. ©) @ 2) oe eee ‘ 


Fig. G5 shows a family of curves for Eq. 5 with values of e of from 0.001 | 
0.050, somewhat less than the ranges indicated on Fig. G4. The effect of af 
plying various ranges of times tothe lateral dispersion (w2 -w1) was examine 
throughout a large range of values by means of an IBM 650 computer.27 Fro: 
this it was determined that the functions are adequately represented by tl 
curves shown. 

Fig. G5 may be used to indicate the width w2 when w] and the elapsed tin 
(tg=tj) are known. In the example shown, if w1 = 1.6 x 105 cm = 1 statute mil 
and a value for e of 0.005 is assumed, the width wg after 24 hr will be 7.0 
105 cm or 4.3 miles. 

The 1-mile width was chosen since it is a reasonable first approximation: 
the initial width of the fieldfromthe new 5-mile effluent outfall of the Hyperi 
Treatment Plant. It is expected that initial dilution in the rising column 
be about 100:1, and subsequent dilution will be by horizontal eddy ree 
Similarly, the elapsed-time figure of 24 hr is based upon current observatio 
made in Santa Monica Bay21,25 and represents a typical travel time to shor 
It is of interest to estimate the dispersion of the field as it travels to shor 


thus: 3 
Wy (ty - ty) e Wo : 
1 mile 24 hours 0.0010 1.5 miles 
0.0025 2.5 miles : 

0.0050 4.3 miles 


0.010 8.5 miles 
0.025 28 miles 
0.050 72 miles 


Weighing the indicated values of w2 against experience, itis possible to as 
probable limits to e. 


27 Personal communications from D. E. Madden, IBM Corp., Los Angeles, Ca 
1949, 
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The field from the existing 1-mile outfall widens from 1 mile in width 
ached after about 4 hr from the outfall to 2 miles in width after 24 hr (tg - 
= 20 hr). Diffusion is to the seaward side only. Thus the value of 0.0025 for 
S probably too low. 

The value of 0.010 for e gives a value for w2 of 8.5 miles. Currents in 
ita Monica Bay average about 0.3 knot15 indicating a movement of about 6.3 
tute miles in 24 hr. A field of this sort would occupy a sector with an angle 
about 60°. Observations of water mass movements in the Bay do not support 
sumption of a lateral dispersion this large. Thus e must be less than 0.010. 
Obviously, more observations are required to characterize the dispersion 
a sewage field. The 5-mile outfall will be in operation in the early part of 
30, and some operating data will then become available. However, it is prob- 
ethat the entire field will be submerged most of the time due to mixing with 
d bottom water and thus acquiring the same density as intermediate depth 
ters. There are difficulties in sampling sucha field, and a number of ex- 
Sive surveys at all seasons are required for definitive work. 

Meanwhile, a value for e of 0.005 in Eq. 5 can be expected to permit rea- 
able approximations of the lateral diffusion of the effluent field. Prelimi- 
‘y computations indicate that it is satisfactory for predicting the areal ex- 
t of bacterial pollution.28 

An evaluation of boundary effects is always required in computations of 
persion of wastes. The existence and nature of clearly distinct water masses 
he oceans29,30 suggests that hydrographic boundaries may be essentially 
effective as land masses in providing lateral restraint to an ocean current. 
On a scale which is important in waste disposal, say 0.1 to 10 miles, it is 
dent that observations of diffusion processes may profitably be extended. 
ropriate studies might include analyses of the effects of currents, density 
erences between water masses, vertical stability, and other physical fac- 
s upon eddy diffusion. The author’s laboratory channel may give some in- 
mation. What happens, for example, when particles of a slightly different 
2 Or specific. gravity are used? How would two distinctly different streams 
‘ract in the channel? 


8 “Predicting Bacterial Pollution in Sea Water,” by C. G. Gunnerson, unpublished 
1958. ; 
D «The Oceans,” by H. U. Sverdrup, M. W. Johnson, and R. H. Fleming, Prentice— 


Englewood Cliffs, N.J., 1946. 
ea livemetne Meereskunde,” by G. Dietrich and K. Kalle, Borntraeger, Berlin, 
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THE HYDRAULIC DESIGN OF SPILLWAY BUCKETS2 


Discussion by M. B. McPherson and Ernest E. Brodbeck 


M. B. McPHERSON,! M. ASCE.—The principal purpose of this discussion 
present a comparison of performance between the slotted-type and solid- 
e roller bucket. A comprehensive solid bucket study covering a wide range 
parameter values has been reported for a radial exit at 45°, approach slopes 
l on 1 and 1 on 2, with flow uniformly distributed over the channel width.2 
The “sweepout” condition is identical with a “free-jet,” the term used in the 
id-bucket study. The authors report that bed arrangement was not critical 
she incidence of the “sweepout” condition. Limited tests with a solid bucket, 
ng both a movable bed at different elevations and a fixed floor, indicated no 
asureable affect on surface profiles or performance due to bed arrangement, 
With a solid bucket a “pulsating surge” was obtained for large values of 
‘R(in which hj is the total available head relative tothe bucket invert and R 
he bucket radius). A “pulsating surge” is a pronounced vertical, unsteady 
tion of the standing wave downstream from the bucket, and appears to have 
aé regime features in common with the “diving flow” condition for slotted 
kets. 
Nearly all of the original solid bucket tests were performed with a falling 
water. With a rising tailwater, cessation of a “free-jet” requires a higher 
water level than for the commencement ofa “free-jet” on a falling tailwater. 
pears that the ratio of the rising to falling tailwater limits might be less 
11.1 for hy/R greater than about3 3, A similar hysteresis effect is prob— 
7 associated with both the “pulsating surge” and the “diving flow” tailwater 
its. For the “sweepout” tests the tailwater was “lowered slowly” and for 
“diving flow” tests the tailwater was “raised slowly.” One would expect 
sistence of the “sweepout” condition at higher tailwater levels on a rising 
water andpersistence of the “diving flow” condition at lower tailwater levels 
| falling tailwater. Whether or not the 0.2 ft. added to the “sweepout” depths 
the 0.5 ft. subtracted from the “diving flow” depths of the slotted bucket 
lel data might include compensation for tailwater reversal effects is diffi- 
to say. However, if performance between the two limits can be assured, 
issue of rising or falling tailwater is not pertinent. 
“he use of a variable datum, the tailwater elevation, in nearly all graphical 
sentations by the authors, obviates the need for any trial and error steps 
Re a arene eas oe RSS Yt ne es ad 


October, 1959, by G. L. Beichley and A. J. Peterka. 

Prof. of Hydr. Engrg., Civ. Engrg. Dept., Univ. of Illinois, Urbana, II. 

“A Study of Bucket-Type Energy Dissipator Characteristics,” by M. B. McPherson 
f. H. Karr, Proceedings, ASCE, HY 3, June, 1957. (Corrections: HY 4, August, 
p. 57. Discussion by E.A. Elevatorski, HY 5, October, 1957, p. 33. Discussion by 
Dougherty, HY 1, February 1958, p. 43. 

Paper 1832, October, 1958, p. 41. 
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in the use of their data for design. For the purpose of comparison it is neces 
sary to use a fixed datum. Fig. Ml has been prepared using the same datu1 
as for the solid bucket study, the bucket invert. The flow parameter include 
the flow per unit width, q, in cubic feet per second per foot. The term hp i 
identical with B and hg with T. The slotted bucket model data were taken fro 
a USBR publication.4 The term Tin equals “sweepout” depth plus 0.2 ft. (es 
timated minimum tailwater level for roller action on a falling tailwater). Fa 
a solid bucket and a falling tailwater the approximate limit before a “free-jet 
is reached is at about h, = 0.2 ho. The average-fit Tip curve arbitraril 
drawn for the slotted bucket points indicates a generally higher tailwater re 
quirement than for the solid bucket, as noted by the authors, The term Tma 
equals “diving flow” depth minus 0.5 ft. All values of Tax, up to and incluc 
ing those for “maximum capacity” flows, have been plotted in Fig. M-1; curvy 
(c)-(c) is an approximate envelope of the largest “maximum capacity” (i. 
“design capacity”) flow points. The range between the Tmin curve and cury 
(c)-(c) is limited, but this fact does not reflect adversely on the recommend 
tion that slotted buckets may be particularly suited for lower ranges of tailwate 
depths. In the range between these curves the depth in the bucket, hb = B, ? 
obviously much greater for the slotted bucket. It could be demonstrated thi 
the surge height A is less for the slotted bucket. In using the model data4 } 
was taken as the difference between the headwater pool water level and tk 
bucket invert. (A plot similar to Fig. M-1, but for “sweepout” and “divi 
flow,” is available.) 

Fig. M2 is analogous to Fig. 15, The maximum hj/R for a solid bucket, | 
avoid a “pulsating surge” with a falling tailwater, is represented by line SI 
It appears that a smaller radius solid bucket may be used for q-paramete 
values exceeding about 0.045 (Fy at bucket invert as high as about 8, from tab 
in Fig. 5, reference 2). Smaller radii are indicated for slotted buckets with ¢ 
parameter values less than about 0.045 (Fj at bucket invert greater than abo 
8), representative of comparatively higher hj values. The solid bucket tes 
were directed towards applications withlow medium head spillways; the slott 
bucket appears to be more suited to high head installations, 

The flatter and more quiescent water surface immediately above and dow! 
stream from a slotted bucket would constitute a decided advantage in certa 
applications. Also, satisfactory performance has been demonstrated for tl 
slotted bucket with a movable bed only 0.05 R below the apron lip, which co: 
responds to an elevation about 0.1 R below a 45° solid radial bucket Lip. It 
questionable whether a solid bucket bed at 0.1 R below the lip would be sati: 
factory under all conditions, particularly at low tailwater levels; a bed elev: 
tion at the same level as the invert (0.29 R below the lip) has been recommen 
ed for general usage.2 : 

While itis acknowledged that “unsymmetrical spillway operation... (mak 
the solid) bucket undesirable in some installations,” is it to be inferred th 
this same limitation is not shared by the slotted bucket? Indeed, can eith 
type of bucket be used with other than uniformly opened crest gates? F 


the information given in the paper the reader has little basis for agreeing " 


4 4. “Slotted and Solid Buckets for High, Medium and Low Dam Spillways,” Prog 
Report II, Section 7, Bur. of Reclamation, Hydr. Lab. Report No. Hyd.-415, Den 
July 1, 1956. b. “Hydraulic Design of Stilling Basins and Bucket Energy Dissipator 
U.S.B.R. Engrg. Monographs, Denver, No. 25, September, 1958. (Both 2-a and 2 b 
tain tabulations of slotted bucket test data.) 
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PLOTTED POINTS ARE FROM USBR REPORT HYD.- 415, 1956 


(h, VS hg CURVES FOR SOLID BUCKET ARE FROM 
FIG.3, ASCE PROC. PAPER I266) 
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FIG. M1.—SLOTTED BUCKET VS, SOLID BUCKET CHARACTERISTICS, 
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sts showed the slotted bucket to be Superior to the solid bucket in all re- 
2cts” since the promise set forth in the Synopsis, that the slotted bucket “is 
npared with a solid-type bucket,” was not fulfilled. Specifically, what solid 
‘ket tests were performed, with what bed conditions, and where can the data 
acquired? Data from the Grand Coulee Dam model should not be used for 
nparison since the exit was not radial. The writer is impressed with the 
lity of the model tests reported, but cannot agree that the tests covered “a 
nplete range of bucket sizes and tailwater elevations, ” Comprehensive data 
flow-tailwater-roller relationships between the limits Tmin and Tmax are 
available. Large values of the q-parameter have been tested for only small 
‘R and small values of the q-parameter for large hy /R (or corresponding 
‘ameters). Fig. 21 purports to be rather universally applicable, but it ap- 
irs that only X = 5 ft was used in all the tests in which the surge height, A, 
3 measured, Since all the tests were performed with an approach slope of 
on 7, what is the permissable extrapolation of the results to other slopes? 
The curves in Fig. 22 appear to be identical with those in Fig. 15, pages 
3-3, reference 5, However, in the latter a few points for Shasta Dam and 
ind Coulee Dam (prototype) have been included, and the ordinate is Z (fall 
m reservoir headwater level to stilling basin floor), with the effective head 
the theoretical floor velocity equal to (Z - H/2) rather than the (H + h) in 
» 22, There appears to be a possible inconsistency between these figures, 
vever, there is good agreement between Fig. 1 of reference § and Fig. 22 
estimated depth and velocity on a spillway face. The curves of Fig. 15, 
erence 5, have been modified by Ven Te Chow, M. ASCE,’ so that the esti- 
ed floor velocity may be read directly. 
in Table 4 the value of F applies to conditions at the shifting tailwater ele- 
on. Unfortunately, the conjugate depth in Table 4 has been computed using 
'F, This is misleading. The values of Dj and Vj used to solve Fy should 
aken at the bucket invert for an equitable comparison of tailwater require - 
its, The equation at the bottom of page 34 cannot be applied otherwise. For 
le M-1 the average of Tyjp = 67 ft and Tmax = 71 ft, or 69 ft, at the design 
‘(901 cfs/ft in Table 3) was used to set the bucket invert at El- 3045. The 
e other values of T are the difference between the tailwater elevations 
n in Table 3 and El. 3045. The conjugate depths shown in Table M-1 should 
yide a better basis for comparison, 
teferring to Table M-1 and Fig. M-2 or Table 3 and Fig. 15, it should be 
‘d that the calculations for design flow conditions given for Angostura Dam 
stitute an extrapolation beyond the test data. In Table 3 a minimum R of 
tis estimated; a minimum R of about 50 ft is indicated in the extrapolated 
re of Fig. M-2 (hy/R = 2.7 maximum). How can a 40-ft radius for Angos- 
Dam be the actual size of the existing structure, the “recommended” ra- 
in Table 3 and then be changed arbitrarily to 47 ft in Table 4, for the de- 
discharge? The actual invert of Angostura Dam is at El. 3040 for which 
lesign discharge q-parameter would be 0.097 and hg/h1 would be 0.53, with 
74 ft, a conjugate depth of 67 ft and hj/R of 3.5 using the actual 40-ft 


“Hydraulic Design of Stilling Basins: High Dams, Earth Dams, and Large Canal 
stures (Basin II),” by J. N. Bradley and A. J. Peterka, Proceedings, ASCE, HY 5, 
ber, 1957, pp. 1402-7 and 1403-3. 

“Air Entrained in Fast Water Affects Design of Training Walls and Stilling Ba- 
” by D. B. Gumensky, Civil Engineering, Vol. 19, No, 12, December, 1949, p. 35. 
“Open Channel Hydraulics,” by Ven Te Chow, McGraw-Hill, N.Y., 1959, pp. 382- 
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bucket radius. As may be seen in Fig. M-1, there is no realistic basis fc 
comparing the invert 3040 and 3045 from the standpoint of maximum allowab 
tailwater, without gross extrapolation of curve (c)-(c) well beyond the graph. 


TABLE M1 (Disc. 220 


Revision of Tables 3 and 4, Paper 2200, for Angostura Dam 
with Crest at El, 3157.2 and Bucket Invert at El. 3045, 
Bucket Radius 40-feet.** 


q,cfs /ft. 901*** 657 365 146 
hg,ft. (=T) 69* 61 50 39 
Z,ft. (a) 153.1 146.0 136.5 125 
Vp ft./sec. 92.5 91.2 89.5 87, 
Vir (a) 0.97 0.96 0.95 0.9 
Vj). /sec. (c) 89.7 87.5 85.0 78. 
Dj, ft. (c) 10.05 7.51 4.29 1.8 
hy,ft. (b) 134.6 126.5 116.3 97. 
hy /ny 0.51 0.48 0.43 0.4 
hy/R 3.4 3.2 2.9 2.4 
a 0.10 0.081 0.051 0.0 
Ve m3? 
Fy (c) 4.98 5.62 7.23 10. 
Do, ft. (c) q 66 56 42 26 


(Teonj) 


Notes: (a) See Fig. 15, Paper 1402, October, 1957. j 
(b) Net, or effective, total head relative to bucket invert. : 
(c) At bucket invert. k 
*Ave. of Trin and Tmax, max. Q, Table 3, Paper 2200 , 

**Actual Bucket Radius, 40-ft.; Invert at 3040 

***Design Flow ’ 


It is hoped that the authors will present pertinent details and a summ 
results obtained with the slotted bucket models for Angostura Dam and > 
Superior-Courtland and Cambridge diversion dams,®as well as similar infc 
mation for whatever solid bucket tests have been performed. 


PROBABLE CORRECTIONS 


The word “not” should be stricken from the first line of page 205" 
In Figure 15 for R = 12-in, the flows should be 3.5 and 3.25 rather tha: 


3.00 and 2.0 given. The units of flow should be cfs/ft. rather than cfs. 
reference 2). . 


ee Energy Dissipators,” by E. A. Elevatorski, McGraw-Hill, N.Y., 1 
pp. 175-178. 
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The terms to the right of the 1/2 in the equation at the bottom of page 34 
Juld be in brackets or parentheses, 

On the last line of page 25 the relation beneath the square-root should in- 
ide (H + h) rather than (H + H). 


ERNEST E. BRODBECK,? F. ASCE.—The authors merit the appreciation of 
se concerned with the hydraulic design features of energy dissipators. The 
iter has always felt that a paper presenting hydraulic design criteria based 
theoretical analysis or model test data should be augmented by an example 
wing its practical application. This has been accomplished, in part, by the 
hors. 
In the design of hydraulic structures, particularly large dams, it is essen- 
| that the designing engineer envision operating conditions and evaluate their 
ential effect onthe performance of the completed structure. In this regard, 
authors have given cognizance to such conditions by stating, “Before adopt- 
a design, all factors which might affect the tailwater range should be in- 
tigated; that is, large or sudden increases in spillway discharge and effects 
lischarges from outlet works or power plants,” However, a more detailed 
cussion of “all factors” should have been included to assure proper recog- 
on of operational features and conditions as they affect bucket performance, 
With the above in mind, the writer offers the following comments with re- 
d to bucket performance vs, operating conditions. For example, consider 
ungated overfall spillway with a flood flow discharge. Under such condi- 
1S, the tailwater rise or fall is generally compatible with the spillway over- 
v. Therefore, the design data for minimum and maximum tailwater as pre- 
ted in the paper would be directly applicable. Howeve r, with a gated spill- 
structure, unusual operating conditions could be created. During a period 
reservoir rise it may be incumbent upon the project operator to open the 
Iway gates for flood flow release. Required sudden gate openings with in- 
ased discharges could result in a condition wherein the magnitude of the 
lway overflow could be considerably greater than the ability of the exis.ing 
water to maintain or confine spillway overflows to roller action in the 
set. Asa result of the high discharge and deficient tailwater depth, sweep 
would occur as shown in Fig. 7A. 
similarly, a sudden decrease indischarge by closing of spillway gates, when 
ailwater is atahigh stage, could result in a diving jet similar to that shown 
"ig. 3 and 7C. Such a condition could prevail as a result of a sudden de- 
ise in spillway flow being released upon a tailwater depth established by a 
-eding greater flow release. 
.S both sweep out and diving flow represent undesirable flow conditions 
a potential for channel scour, it is incumbent upon the hydraulic engineer 
1alyze carefully the spillway gate operating schedule so as to maintain over- 
S and tailwater depths within the operating range required toassure proper 
sr bucket action. In order to accomplish this, a knowledge of the rate of 
and fall of tailwater depth vs. discharge is required. Also, the problem 
egradation or downstream channel changes which may affect tailwater 
2s should not be overlooked. 
here remains another design feature which is often not given too much 
ideration. This feature concerns itself with the establishment of the height 


-Hydr. Engr., North Atlantic Div., Corps of Engrs., New York, N.Y. 
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radius of slot C apparently does not prove that the corresponding incipie 
cavitation parameter would be little affected. Since “field experience sho 
cavitation below a slot with such a 12:1 convergence” and since the lowe 
boundary pressure was found below the slot it is implied that damage dow 
stream from the slot is necessarily associated with the minimum bounda 
pressure, However, it is not clear whether sharp or rounded upstream a 
downstream corners or some combination are at issue. With the W.E. 
rounded upstream corner tests, the initiation of cavitation occurred within t 
void of the slot. If cavitation is also found to commence within the slot, w; 
the sharp upstream corner of the Type C slot, the minimum boundary pressu 
occurring at the end of the convergence may not be a Significant criteria f 
anticipating incipient cavitation. As to field damage below the slot, perha 
vapor produced in the slot implodes upon reaching the low pressure zone 
the P.I. It is not clear whether or not the author has conducted any cavitati 
tests on the Type C slot. What are the dimensions of the slot with “the ov 
wardly offset, 1-inch-rounded corner with a 12:1 convergence”? Does t 
term “cavitation pressures” used in the same sentence mean actual inciden 


TABLE M2 
Velocity Offset 
greater height, y, 

than (f.p.s.) (in.) 
45 1/16 

45 1/8 

30 1/4 

20 1/2 


of cavitation in the tests? If so, at the point of minimum observed boundaj 
pressure? ; 


Into -the-Flow Offsets.—The writer is particularly interested in Fig. : 
n 


Using 31 piezometer taps, Williams5 made extensive pressure distrib i 
measurements for conditions including abrupt into-the-flow offsets in aclos 
conduit and an open channel. The closed conduit tests covered 13 values o: 
ranging between 0.06-in. and 1/4-in. for an approach conduit height of 0.500 
0.833 ft, with a maximum ratio of offset to conduit height of 0.04. The approa 
section, between a transition piece leading from a large entrance tank and | 
offset, was 12-in. long; the conduit was 6-in, wide. Available flow limited : 


5 “Study of Misalignments in an Open Channel and a Closed Conduit,” by J. C. | 
liams, Jr., M.S. Thesis, Lehigh Univ. Library, Bethlehem, Pa., May, 1952. 
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roach velocities to a maximum of 24 fps. His results can be expressed in 
‘rms of the parameter given in the following table, in whichH is the approach 
lezometric head, hy is the minimum piezometric head (3/64-in. downstream 
om face of offset) and Vo is the average approach velocity. Values of the 
rameter approached and reached a constant asthe approach velocity was in- 
reased. Sample asymptotic values from his Fig. 12 are given in Table M3. 
V (velocity above offset, as in Table M2 was used 

Stead of the approach velocity Vg, the above pa- TABLE M3 

meter values would be increased by a maximum 


9% (for the 1/4-in. offset and 0.50-ft. approach 

ight). Referring to Table M2, it is to be noted . Offset H - hy 
at the incipient cavitation parameter values are NCR if V.2/2 
eater than the pressure parameter values by ; geic® 
illiams, typical of abrupt flow-disturbances. One 1/16 0.66 
ason the pressure parameter values are lower is 1/8 1.06 
cause thetrue minimum pressure cannot be mea- 1/4 1.12 


red, The minimum pressure occurs at the offset 
rner, the initial point of flow separation. The pressure gradient immediate- 
downstream from the offset is too steep for satisfactory extrapolation of 
otted boundary pressure distribution data to the offset corner. 

Holi6 has reported on cavitation tests in water tunnels with triangular and 
rcular-arc offsets. The triangular elements were 90° isosceles with one leg 
rpendicular to the floor and the other facing downstream. The triangular ele- 
snts were mounted on a flat plate with an airfoil shaped leading edge. The 
ute was installed in the central plane of the circular tunnel test section. For 
hase 1” the test elements were placed 11-in. from the leading edge of the 
ite (“Forward Station”) and at 37-in. (“Center Station”). “One would expect 
it the pressure field would be dependent upon the location of the element in 
> boundary layer,” and “Also, the shape of the boundary-layer velocity pro- 
2 would influence the pressure field.” Boundary-layer velocity profile mea- 
rements were made for all triangular element tests. In Table M4 are ap- 
oximate incipient cavitation parameter values for one series of tests(Phase 
with velocities of 30 to 60 fps, where H is the approach piezometric head, 
is the vapor pressure and Vo is the approach velocity. If one assumes that 
streamline near the point of separation of the triangular elements is simi- 
‘in curvature to that for a rectangular-shaped offset, a direct comparison 
in order. The general magnitude of parameter values in Tables M2 and M4 
> in reasonable agreement. As the boundary layer becomes small compared 
h the offset height, the cavitation parameter in Table M4 increases. It may 
seen that a given ratio of offset height to boundary layer thickness (1.1) is 
sociated with a different cavitation parameter value at the two stations. The 
io of boundary layer displacement thickness to momentum thickness is 1.40 
the Forward Station and 1.33 for the Center Station, a boundary layer shape 
ameter. “For given values of velocity and boundary-layer shape parame- 
, the incipient-cavitation number increases with relative height of rough- 
s. Furthermore, the relative height of roughness appears to be the most 
nificant variable for describing the incipient-cavitation number.” 
Williams® reasoned that the local velocity approaching an offset wasa con- 
(ling factor, with the boundary-layer velocity profile thereby governing the 


3 “The Hiception of Cavitation on Isolated Surface Irregularities,” by J, W. Holl, 
ME, No. 59-Hyd-12 (to be published in the ASME Journal of Basic Engineering). 
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of training walls. The model tests presented in the paper were accomplishec 
in a flume which voided the condition of training wall overflow from tailwate1 
back-up. As the training walls at abucket are seldom constructed to the maxi- 
mum height of tailwater, some overflow at this location could prevail. If suck 
overflow were of a significant magnitude, a diving jet adjacent to a training 
wall could occur with consequent damaging scouring action concentrated at the 
toe of the wall. Such bucket action could be avoided by raising the wall to sucl 
a height so that training wall overflows would be insignificant insofar as the 
stability of bucket roller action is concerned. Proper operation sequence 0; 
spillway gates adjacent to the wall could also alleviate such undesirable floy 
conditions. 

The foregoing is offered indiscussion only to emphasize that the final bucket 
design should be thoroughly investigated and anlayzed in the light of conditions 
that may prevail during normal and emergency spillway operation. 
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HYDRAULIC CHARACTERISTICS OF GATE SLOTS2 
Discussion by A. Thiruvengadam, M. B. McPherson and Warren H. Kohler 


A. THIRUVENGADAM.!—Fig. 20 shows the head-velocity relationship for 
incipient cavitation for abrupt offsets into the flow. This datum would be much 
more useful if the offset is considered as a single isolated roughness protrud- 
ing into the flow. ‘Then the roughness Reynolds Number would be V z/y where 
V is the velocity, z is the offset in feet, and vy is the kinematic viscosity. The 
relative roughness would become Z/H where H is the depth of flow beyond the 
offset. The critical cavitation parameter 


where Hy is the vapor pressure and g is the acceleration due to gravity. Then 
Fig. 20 would provide data for the plot of the equation 


The writer has plotted this in Fig. Tl. If such a relationship could be used for 
shecking the design of stilling basins, the designers would be much benefitted. 

- The writer records his thanks to N. S. Govinda Rao, F. ASCE, for his sug- 
sestions in preparing this discussion. 


M. B. McPHERSON, 2 M. ASCE.—The major objective of this discussion is 
© solicit further clarification and amplification from the author. The broad 
scope of his excellent paper has obviously required considerable condensation 
f certain selections. 

Gate Slots.—Presentations of pressure distribution profiles are given in 
erms of the “reference pressure immediately downstream from gate leaf,” 
n Fig, 8 through 15 and in Fig. 17. It appears that the cross section “imme- 
liately downstream from the gate leaf” must coincide with the “downstream 
dge of the slot.” The figures in question are presumed to be for two- 
limensional tests without a gate leaf. Therefore, the reference point used is 
;omewhat vague. Is the pressure at the reference point representative of an 


erage or normal gradient, away from the gate slot and unaffected by both 


ee 


2 October, 1959, by J. W. Ball. 
1 Tech. Asst., Civ.& Hydr. Engrg. Sect., Indian Inst. of Science, Bangalore 12, India. 


2 prof. of Hydr. Engrg., Civ. Engrg. Dept., Univ. of Illinois, Urbana, tll. 
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CRITICAL CAVITATION PARAMETER K.- 


02 04 06 08 10 12 14 16 18 20 22 


“5 
OFFSET REYNOLDS NUMBER op x10 


FIG, Tl.-RELATION BETWEEN OFFSET REYNOLDS NUMBER AND CRITICAL 
CAVITATION PARAMETER WITH RELATIVE OFFSET AS THE THIRD 
PARAMETER, 
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né leaf wall-aperture and the seat? Was the conduit cross-sectional area 
ken upstream from the gate slot, at the gate slot, or at some point down- 
tream from the slot, to compute the “velocity head at reference point”? 

In the original data given by Kindsvater3 for “Slot B” of Fig. 11, the refer- 
nce point was taken at “an arbitrary upstream control pievometer,” and the 
normal friction grade line” was extended from this upstream control point 
ast the gate slot position. From Fig. 68 of the reference it appears that the 
ighest smoothed-curve pressure-parameter value would be about +0.20 and 
1¢ lowest about -0.063 (lowest plotted point = -0.070), using the “normal grade 
ne” at the downstream edge of the gate slot and “the computed average ve- 
city head.” What procedure was employed by the author to obtain the +0.19 
nd -0.12 extremes given in Fig. 11 for Kindsvater’s data? The minimum 
ressure for Slot B is shown at a distance of about 6.1-in. downstream from 
1e gate slot in Fig. 11 whereas it is given as about 13.5-in. in Fig. 68 of the 
eference (that is, at the P.I.). 

Incipient cavitation parameter values have been presented by Brown4 fora 
2:1 tapered slot similar to Type B, but for W/D=1.55 Table M1 is presented. 


TABLE M1 


Radius at Taper length, Taper 
ownstream edge of slot offset, 
ige of slot touPsds total** 


1.0425 W 0.250 W 0.0625 W -0.36 
),0425 W 0.500 W 0.0850 W -0.29 
0775 W 0.750 W 0.140 W -0.28 
).1050 W 1.000 W 0.187 W 

).114 W)* (1.228 W)* (0.205 W)* 


*Slot B, Figure 11, W/D = 1.8. 
**12:1 Taper plus downstream corner curve. 


he last column in Table M1 has been expressed in terms of the author’s pa- 
umeter since the reference head and velocity were “at the slot;” the equivalent 
hy was taken asthe vapor pressure head by Brown. “The data presented in- 
cate about the same degree of cavitation for each slot.” In the W.E.S. tests, 
‘avitation flashes first appeared in the eddies within the slot proper.” 

In Fig. 11 and 12 it is shown that with W/D = 1.2, a 12;1 convergence and a 
1arp upstream slot corner, the lowest pressure parameter value was decreased 
ightly from -0.16 to about -0.14 with “no adverse effects” when the down- 
ream from the P.I. (where the lowest pressure occurs) with a 12:1 conver- 
‘nce and a sharp upstream slot corner is about -0,12, as opposed to an in- 
pient cavitation parameter of about -0.4 indicated by gross extrapolation of 
e table inthe above paragraph. For a rounded upstream slot corner incipient 
vitation reference heads are greater the more the taper offset is reduced. 
1¢ fact that the pressure parameter was little changed by omitting the 0,056 


3 Discussion of “Cavitation in Hydraulic Structures—A Symposium,” by C.E. Kinds- 
ter, Transactions, ASCE, Vol. 112, 1947, pp. 101-105. 

4 «Hydraulic Models as an Aid to the Development of Design Criteria,” by Frederick 
Brown, Corps of Engrs., Waterways Experiment Sta., Bulletin No. 37, Vicksburg, 


88., 1951, pp. 21-23. 
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magnitude of the pressure parameter. Unfortunately Williams did not mal 
velocity traverse measurements and his analysis was purely qualitative. EF 
argued that a given sized offset located near the entrance of a prototype cor 
duit, where the boundary layer-is thin, would have a higher pressure parame 
ter value than if it were located downstream in the zone of established flo\ 
The trends in Table M4 appear to lend some credence to this argument, 

For truly representative and generalized results, it will probably be neces 
sary to restrict future offset studies to fully developed boundary layers, © 
make complete velocity profile measurements, and to express the results 
terms of boundary layer characteristics. 

It is hoped that the author will provide the details on the tests represente 
by Fig. 20, in his closure to the discussions, Also, it is suggested that he de 
fine the criteria used for delineating incipient cavitation. 

An interesting comparison might be made between the circular-arc analys 
and accompanying incipient cavitation results presented by Holl and tl 
circular-arc pressure profile data given by Hickox. 7? 


TABLE M4 


Offset Ratio of y H - hy Ratio of y H - h 

height, to mean 5) to mean | 
2 

y; bound. layer Vo" /28 bound, layer Vo /2 

(in.) thickness (approx.) thickness (appro? 


1/16 0.90 0.138 

(1.3)* 
1/4 1.55 0.552 1.42 
1/2 2.2 1.10 1.8 


*By graphical interpolation. 


oe ————OOOOOOO—O—oror—rCoOoOoOnnnanqom eee: 


WARREN H. KOHLER. 8—The information and data given in this paper prt 
vide a sound and authoritative background as well as a progressive approa 
and analysis of the gate slot problem. The accuracy and reliability of the tes 
have been proved by the success of prototype gates which were designed ont 
basis of data secured by the model tests, This success was not achieved, hov 
ever, on the basis of random application of available data to designs. To a co! 
siderable degree the success of prototype gates is due to the close coordin 
tion and cooperation which exists between Laboratory and Design Division 
the Bureau of Reclamation. The designing and building of prototype gate 
based on small scale models, necessarily involves some changes and compr 
mises between the ideal conditions in the model and what is practical and ° 
tainable on the prototype from the design, fabrication, and cost standpoint 
This adaptation requires care, ingenuity, and continued consultation with t 

: = 


T« 
Cavitation in Spillway Tunnel,” by George H. Hickox, Fontana Proj. Hydr Mor 
Studies, T.V.A. Tech. Monograph No. 68, Knoxville, Tenn., Chapter 4, ee : 


ie Head, Large Gates and Valves Sect., Mech. Branch, Bur. of Reclamation, Denve 
olo. . . 
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aboratory to avoid deviations from the test models which could seriously af- 
ect operation of a prototype gate. This projection of the data from a 6-in. 
aodel to a gate which may be 6, 8, or 10 ft in size does pose a number of de- 
ign problems. The following comments are primarily intended to point out 
ome of the design problems involved intranslating the model to the prototype. 

While not mentioned in the paper, the problem of providing an economical 
ate design with sufficiently smooth fluidways to avoid local areas of cavita- 
ion under high-velocity flow posed a design problem. Completely machined 
luidways would have been hydraulically satisfactory, but the machine work 
nd extra joints required for machining would have been very expensive. The 
esign problem was solved by changing gate bodies from the conventional cast 
esign to welded design. The change avoided the necessity of machining fluid- 
jays, as the rolled surface of the plates forming the fluidway were smooth 
nough without requiring additional machining. Some problems in welding gate 
odies to avoid undue distortion were encountered, but numerous bodies with a 
luidway plane tolerance of 1/16-in. in 3-ft have been successfully and economi- 
ally fabricated. The use of welded construction makes curvedor other shapes 
ownstream from the gate slots economically feasible, as bending the plate to 
uit atemplate curve is relatively simpie and cheap when compared to the cost 
f machining such a surface. 

As a matter of fact, one of the basic reasons for considering a small offset 
n the downstream side of gate slots was that it would be economically prohibi- 
ive to try to match the openings of the upstream and downstream bodies of 
ates exactly. This led to the investigations which showed that a small offset 
ctually was ahydraulic improvement at the slot. The faring of this offset back 
) the initial fluidway width did, however, pose some problems hydraulically 
hich the model tests solved. 

The elimination of the radius on the downstream edge of the gate slot was 
articularly helpful from a design as well as a hydraulic standpoint. Offhand, 
‘would not seem that reduction in the span width of a gate which is 6 or 7 ft 
ide by 2 in. would be material, but in the narrow-slot type of gate where only 
tongue on the side of the gate extends into the slot, the reduction in bending 
,oment amounted to nearly 20% on the critical tongue extensions. 

The author speaks of the use of slot fillers, but from the standpoint of pro- 
type design, these pose serious problems. Ifa slot filler does not aline al- 
ost perfectly withthe gate slot sothat no edges protrude into the high velocity 
ow, it canbe an even worse problem from a cavitation standpoint than an open 
lot, Unless the designer is prepared to pay the penalty involved in design and 
linement costs for providing almost perfect alinement, slot fillers should be 
voided. 

The author describes the improvement which can be achieved by moving the 
ute slot upstream from the spring point onthe bottom of a gate leaf to provide 

Cavitation-free slot for both free discharge and back pressure conditions. 
he problem in this design is the seal; and while several designs which possi- 
y might have been successful were investigated, the basic seal problem is 
ta simple one. The complication lies chiefly in developing a seal which is 
ifficiently rugged for high heads and has enough flexibility to seal in three 
anes involving right-angle corners. The gates for which this design was 
rst proposed had to withstand heads in excess of 400 feet. For this reason, 
tried and questionable seals were not acceptable, and the conventional ar- 
ingement and seals were used, It may be that some new thinking will result 
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ina satisfactory seal which can be tried on a lower-head gate. For the tim 
being, however, the design has fallen victim to the gate designer’s axiom, 4 
you can’t seal it, forget it.” 
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CAVITATION DAMAGE OF ROUGHENED CONCRETE SURFACES@ 


Discussion by A. Thiruvengadam 


A. THIRUVENGADAM.!—This paper deals primarily with the inception of 
avitation on rough surfaces. The author deserves to be congratulated for his 
ttempt to correlate the inception pressure with shear velocity which is defi- 
itely a new approach. 

The writer wishes to clarify one point. The title of the paper is “Cavitation 
amage on Roughned Concrete Surfaces.” But the paper deals only withthe in- 
eption of cavitation. The cavitation inception and cavitation damage are two 
ifferent aspects of the cavitation phenomenon obeying different physical laws. 


VARIATION OF RATE OF DAMAGE W.R.T. TIME 


Mix.No.7. Comp. strength: 47 50 bs/Sq 


Time in Hrs 
EIGSTIF 


The experiments conducted by the writer at the Indian Institute of Science 
aboratories have clarified certain points on the effect of surface roughness 
1 cavitation damage. The experiments were conductedina venturi-type cavi- 
tion tunnel similar to the author’s tunnel. Fig. T1 shows that the intensity of 
image (rate of weight loss per hour of test) is maximum at the beginning of 


2 November, 1959, by Donald Colgate. ‘ 
1 Tech. Asst., Civ. & Hydr. Engrg. Sect., Indian Inst. of Science, Bangalore 12, India. 
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RATE OF DAMAGE ON HIGHLY POLISHED STONE W.R.T. TIME 


St Per 3 Hours — In Ounze 


Time in Hrs. 


PIGS T2s 
EFFECT OF SURFACE ROUGHNESSON DAMAGE 


ROUGHNE SS 


FIG, T3. 
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e test on a natural concrete surface (compressive strength of the concrete is 
10wn in the figure itself). The rate of damage falls down and reaches a con- 
ant value after a particular time depending on the compressive strength. A 
ranite block was polished till the surface became glassy in appearance, and 
sted. The compressive strength of the granite tested was 16400 psi. Fig. T2 
1ows how the damage develops and falls down after prolonged testing. Then 
fferent surface roughnesses were produced on the same granite and tested 
x three hours each. Fig. T3 shows the variation of the rate of damage with 
e a roughness (see Table T1 for definition of surface roughness 
lopted). 


TABLE Tl 


Definition Description as per public works 
department code 


1¢ line dressed Chisel dressed so that no portion of the face dressed is 
more than 1/4 inch from the edge of a straight edge laid 


along the face of the stone. 


Chisel dressed so that no portion of the face dressed is 
more than 1/8 inch from the edge of a straight edge laid 
along the face of the stone. 


vo line dressed 


Chisel dressed so that no portion of the face dressed is 
more than 1/16 inch from the edge of a straight edge 
laid along the face of the stone. 


ree line dressed 


The surface of the stone is dressed until a straight edge 
laid along the face is in contact at every point. 


ur line dressed 


The surface is polished till the surface becomes glassy 
in appearance. 


nely polished 


From these data, two points are clear. One is that the intensity of cavita- 
mn damage on polished surfaces increases withtime and then decreases. The 
itial increase in damage is due to the progressive cracking on account of the 
peated collapse pressures of the bubble. The decreasing trend may be due 
the cushioning out of the shocks by the water filling the pits formed. There- 
re it is clear that polishing the surface is only a temporary remedy. The 
cond point is that the roughness as such has no effect on cavitation damage 
yond a certain limit. 
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EARLY HISTORY OF HYDROMETRY IN THE UNITED STA TES2 
Discussion by J. W. Johnson 


J. W. JOHNSON,°9 M. ASCE.—During the 1860’s and 1870’s when hydraulic 
ning activity was at its height in California, several historically important 
gineering developments occurred. Some of these projects were the first of 
sir kind in the west, and a few of them were the first of their kind in the 
rid, all taking place in a relatively small area of the Mother Lode region. © 

y name some of the more important, it was here that the first large-scale 
velopment of reservoirs and canals was made in the United States for min- 
y, irrigation, and power purposes; the first use in the country of iron pipe 
der heads of several hundred feet; it was the birthplace of the Pelton wheel, 
d here occurred the invention of the needle nozzle; in this area was one of 
2 first plants in the west for the generation of electric power under high 
ads; the first long-distance, high voltage, power-transmission line in the 
rld; and the first successful long-distance telephone line in the world. 
any of these historic engineering works are still in operation. 

An important contributor to the developments during this era was the well- 
own hydraulic engineer, Hamilton Smith, Jr., who, in the writer’s opinion, 
ght well have been mentioned in the paper. Smith is best known for his 
ok61 “Hydraulics,” the background for this treatise being obtained while he 
$s connected with the North Bloomfield Gravel Mining Co. While at North 
oomfield, Smith became the recognized authority on all matters relating to 
draulic mining in California. Between August 1871, when he was appointed 
perintendent of the North Bloomfield company, and 1881, when he left the 
ited States for mining work in South America, 62 he not only designed and 
nstructed large dams, pipe lines, and tunnels, but he also conducted impor- 
it experiments on orifices, weirs, water wheels, and on the flow of water in 
ye lines. Smith’s paper, “Water Power with High Pressures and Wrought 
m Water Pipe,” was awarded the American Society of Civil Engineers’ 
omas Fitch Rowland Prize in 1884.63 The results of his experiments on 
fices and weirs still are quoted in present-day hydraulics textbooks. 64 His 
ta on flow in pipes were obtained by tests conducted on such pipe lines as 
» Texas Creek and Humbug inverted siphons of the North Bloomfield sy stem65 
| NS ee a a ee ee a 
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59 Prof. of Hydr. Engrg., Univ. of California, Berkeley, Calif. 
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64 “Hydraulics,” by G. E. Russell, New York, 1942, p. 125 
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and the Cherokee siphon in Plumas County. To provide more precise data 
friction losses than was possible with experiments on large pipes, tests we 
made by Smith on relatively small pipes at the famous New Almaden quicks 
ver mine in Santa Clara County. 

The orifice and weir experiments by Smith were conducted at Columbia k 
and North Bloomfield and were financed by the three large companies supply 
water to the mines in that region, The purpose of the experiments was to p1 
vide precise discharge coefficients for the type of orifice that was usedin n 
tering water to the various mining operations. A complete description of © 
details of the tests is given in Smith’s “Hydraulics.” 

Hamilton Smith also is credited with giving modern tunnel building its st 
when he built the mile-and-a-half Malakoff tunnel at North Bloomfield in 1! 
to explore and work gold-bearing gravel in an ancient stream channel. Pr 
to construction of the Malakoff tunnel, tunnels generally were excavated fri 
one end, sometimes from two ends. In the Malakoff tunnel, Smith worked fr 
both ends but, in addition, sand eight shafts and worked both ways from ea 
thus giving sixteen faces to work on simultaneously. Examination of Smit! 
personal papers, now in the University of New Hampshire Library, shows t 
much information on tunnel construction was obtained in 1871 by correspor 
ence with Edward Frost of Frost Brothers, consulting engineers of Bost 
who gave Smith considerable data on costs, rate of progress, and other d: 
obtained in constructing Hoosac tunnel on the Troy and Greenfield Railroad. 

With the mounting opposition to hydraulic mining, Smith left California 
1881 for a successful mining career in South America, Africa, and Lond 
Finally he formed a partnership in New York with Henry C. Perkins, who k 
worked with Smith at North Bloomfield. Smith died on J uly 4, 1900, at his ho: 
in Durham, New Hampshire. His will provided funds for the library of 1 
University of New Hampshire, which is known as the Hamilton Smith Libra: 
Unfortunately, only a small portion of his personal books and papers is ava 
able to yield further information on a man who found time to perform class 
engineering research at a time when the gold output of the mines was the i 
portant item to most people. 


NEW APPROACH TO LOCAL FLOOD PROBLEMS2 
Discussion by Bernard L. Golding 


BERNARD L. GOLDING,2 M. ASCE.—The author and the Tennessee Valley 
thority are to be congratulated for the action they have taken in calling the 
ention of both the public and Congress to this “positive method” of flood con- 
1, that is, the inexpensive method of flood prevention by the establishment of 
€am encroachment lines or flood plain zoning rather than the expensive 
thod of flood control dams, levees, etc., after a damaging flood has occurred. 
sy are also to be congratulated on the action they are taking in the Tennes- 
» Valley watershed area in supplying the local communities with the neces- 
‘y hydrologic and hydraulic information they need to setup a system of flood 
in zoning. 
The communities of the Tennessee Valley watershed area are fortunate that 
overnment agency is available intheir areato supply them with the data ne- 
sary for the establishment of such zoning. As mentioned by the author, in 
ny states there is no such agency, either Federal, State or local, for which 
amunities may go to for help, even though they may recognize the need to 
e action. However, the lack of such an agency should not deter acommunity 
a county from the establishment of such flood plain zoning or regulation. 
2 actual hydraulic and hydrologic problems involved in determining the ex- 
t of the areas of the flood plain that will be inundated by a particular design 
)d are not as difficult as was indicated by the author, such that a major 
ncy must be available to solve them. Enoughtechnical information is avail- 
e from various sources, so that local engineering forces or local engineer- 
forces supplemented by consulting engineers can do the job effectively. 
st communities already have the necessary legal assistance to draft the 
essary legislation once the technical problems have been solved. Counties 
has Fairfax County in Virginia and Westchester County in New York have 
d their engineering forces supplemented by consulting engineers to do ef- 
tive flood plain zoning work in their areas. 
Technical information necessary for the solution of hydraulic and hydro- 
ic problems is readily available to the engineer from many sources. Dis- 
rge information is available from the U. S. Geological Survey who have 
lished a series of papers on regional flood frequency, in which discharges 
various frequencies have been correlated with the drainage basin charac- 
stics for rivers and streams in a particular area. Discharge information 
Iso available from the Bureau of Public Roads which has made similar re- 
ial correlations. Methods for the computation of water surface curves 
‘kwater curves) have been presented in many technical publications readily 
| January, 1960, by Herbert D. Vogel. Le 
' Head, Hydr. Dept., Howard, Needles, Tammen & Bergendoff, New York, N.Y. 
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available to any engineer. An excellent, simple, easy-to-understand met! 
for computing water surface profiles has been presented by Joe M. Lara, 
ASCE, and K. B. Schroeder.3 An entire manual for the computation of wa 
surface profiles written by the Corps of Engineers is also available.4 Paps 
on headwater losses through bridges are also available from boththe Geolo 
cal Survey and the Bureau of Public Roads. 

The actual cost of a flood plain zoning program to a community is qu 
small in comparison to the cost to the community that would result if eve: 
small housing development or several industrial plants constructed on the flc 
plain were damaged by a flood. Besides the cost of engineering services, | 
only additional expense is the cost of doing the necessary survey work wh: 
generally consists of cross sections necessary for the hydraulic computatic 
and any necessary boundary or real estate surveys necessary to establish 1 
actual zones or encroachment lines. Most communities or counties genera 
have full time survey parties which can do the necessary survey work dur: 
the times when they are not busy on other work. The cost of aerial sury 
work is not behind the financial capabilities of even small communities anc 
may be used to obtain cross sections and show the extent of zones, 

The necessary hydraulic and hydrologic computations done in conjuncti 
with a flood plain zoning program in small communities willin many instance 
provide other indirect benefits. Any existing bridge waterway areas which a 
too small so that headwater losses through them could cause extensive upstrez 
damage during flood periods will readily become apparent. The results of su 
computations will also enable local engineers to design adequate bridge wate 
way areas for new local or county bridges, which in many instances are bac 
underdesigned in this regard. Often the cost of replacing a single county brid 
owing to a washout will be only slightly less than the cost of an entire progr 
of flood plain zoning for a small community. 
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3 “Two Methods to Compute Water Surface Profiles,” April, 1959, 
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